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Das Verfahren von Theodorsen?) ist’in seiner ursprünglichen Form auf die konforme Abbildung eines 
kreisnahen Gebietes auf den Einheitskreis zugeschnitten und nur auf gemischtem graphisch-analytischen 
Weg durchführbar, wenn die Rechenarbeit in erträglichen Grenzen bleiben soll. Dagegen läßt sich für die 
umgekehrte Abbildungsaufgabe (Einheitskreis > kreisnahes Gebiet) der Prozeß auf rein rechnerischem 
Wege durchführen und so bei Bedarf eine beliebige Genauigkeitssteigerung erzielen. Für diesen Fall werden 
von einem neuen Gesichtspunkt aus die zugehörigen Formeln abgeleitet und Überlegungen über die Kon- 
vergenz des Verfahrens angestellt. 


The method of T’heodorsen?) gives, in its genuine form, the conformal mapping of a nearly circular 
region onto the unit circle. It requires a combined graphical-numerical process, if the computation work 
shall remain endurable. It is possible, however, to perform the reciprocal mapping (unit circle — nearly 
circular region) in a purely numerical way and so to achieve an arbilrary increase of accuracy if required. 
In this case, ihe neccesary formulae are derived from a new point of view, and the convergence of the method 
is investigated. 

La methode de T'heodorsen, dans sa forme originale, aboutit & la reprösentation conforme d’une domaine 
a peu pres circulaire au cercle unitaire et se ne peut executer que par un procede graphique-analytique 
combine, d moins que la calculation n’exige pas trop d’effort. La contre pour la repr£sentation r&ciproque 
(cercle unitaire > domaine a peu pres circulaire) on peut determiner une methode purement analytique; 
ainsi on alleint, suivant le besoin, une pre£cision quelconque. Pour ce cas on d£rive les formules necessaires 
d’un nouvel point de vue; la convergence de ceite methode est examinee. 


Merox TeopopcenHa?) ımpumeHnm B CBeii IIePBOHayalıbBHOÜ PopMe IA KOHPOPMHOTO OTOÖPa- 
RT %KeHHA KPYTOOOPasHBIX OÖNacTei Ha eNCHHYHEIH KPyT HU MOIKeT ÖBITB OCYIHeCTBIeH TOABKO 
CMEIHAHHBIM TPAhuYecKH-AHAJNUTHYeCKHUM IIYTeM, eCHIH BEIKIANKH He MOIKHBI ÖBITb CJIMIIIKOM 
CJIO®KHEIMH. HaıporTuB, oöpaTHaa sayaya IpeoßpasoBaHua (eNeHMYHBIÄ KPYyT > KpyrooöpasHaa 
OÖNACTB) MOLKeT ÖBITB PeImeHa YHCTO AHANNUTUYeCKHM HYTEeM, 4YTO IMO3BONAeT IOANYYHTB 
HCKOMBIÄ pesysIBTaT c I00K Tpeöyemoä CTeNeHL® TOYHOCTH. 1A 3Toro cHyyaa BEIBONATCH 
N HOBBIM IIYTeM He00XonuMbIe POPMyJIbI U HCCHEeNyeTcA CXONUMOCTb HAHHOTO MEeTopa. 


ft 


I. Das Verfahren 


Bi _ 1. Aufstellung der Grundgleichungen 
" Für die Näherungsfunktion der Abbildung wird in Analogie zu Theodorsen der 
- Ansatz 
N 
wfnk) zit me Hain". ..2.2.2.22.2. (ll) 
m=0 


gemacht; f®(z) soll das Innere des Einheitskreises der z-Ebene näherungsweise auf das Innere 
des kreisnahen Gebietes @ der w-Ebene konform abbilden®). 

Dabei wird k®(z) dadurch bestimmt, daß durch f(z) möglichst viele Punkte 2* = e'? 
der Peripherie des Einheitskreises exakt auf den Rand R von @ abgebildet werden sollen. Es 
zeigt sich, daß dies für 2n Randpunkte möglich ist, wenn man 2n Punkte z£(k =1,2,...2n) be- 
liebig fest auswählt und nur verlangt, daß ihre Bilder w; auf Rliegen sollen, ohne also deren Lage 
genauer festzulegen?). 


1) Unter Benutzung der ekrtation des Verfassers (eingereicht bei der Pädagogischen Fakultät der 
TH Dresden, Referent: Prof. Dr. Fr. A. Willers, Korreferent: Prof. Dr. O.-H. Keller. Tag der münd- 
lichen Prüfung: 27.7. 1949). 

"#.2)T.Theodorsen, NACA-Report Nr. 411, me) und T.Theodorsen.und J.E. Garrick, 
©: teport Nr. 452, 1933. 

ist für die hier ‚angestellten Überlegungen gl eichgültig, ob es sich um die Abbildung des Inneren 
e oder des Äußeren auf das Äußere hande lt. Die Rechnung verläuft in analoger Weise, nur daß 
j reiten | Falle ein Polynom‘ in 1/z ist. Hier wird deshalb nur der erste Fall behandelt. _ 

> re diese Mehr: kann man sich ein klares Bild vom Verhalten der Näherungsfunktion machen, 
sher üblichen Ableitung, bei der die pa für die exakte Abbildung ergebende Potenz- 
en ( Be ERROR Sn wird. 

en 


ip 


Ey 
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Obwohl danach die Auswahl der xt völlig frei ist, soll in der Folge nur der Fall betrachtet. 
werden, daß die 2; äquidistant auf dem Kreisumfang verteilt sind, daß also für die AB BE 
von 2t = ek gilt | a ü 
. Dh AD, ‚n-Ao = P% 9 = Re (1.2). D 

Das kreisnahe Gebiet @ ist ein Sterngebiet, dessen Rand R in Polarkoordinaten durch 
die eindeutige stetige Funktion ey) gegeben ist, so daß also die Randpunkte die Form 
w* = o(y)e'Yhaben. Die Ableitung 0’(y) sei stückweise stetig; weitere Einschränkungen für 0’(y) 
E:- werden weiter unten ausführlich besprochen. Der üblichen eindeutigen Festlegung der Ab- 
| bildung dient zunächst einmal die Form des Ansatzes (1.1), derzufolge der Nullpunkt der z-Ebene 
in den Nullpunkt der w-Ebene übergeht; außerdem werde vorgeschrieben, daß f(”’(0) > 0, also 
positiv reellist. Mit c„= a, — ib, ist f'®’(0) = e®—i%, und die Vorschrift wird durch dieFest- 


De ne 1 


na 


setzung b, = 0 erfüllt. 
Gemäß der obigen Vorschrift sollen die Punkte zZ: in Punkte wf übergehen, die an Rliegen; 
& es gilt also 
N i i (n) (z \ 
j { 0 (yr) e Ad = ek . % @ et eco Ile 0 or enge (13er 
R | Division durch e”% und Logarithmieren liefert 1 
: In o@yr) + pe — 9) = Re). | 4 
; Be Bezeichnet man den Real- bzw. Imaginärteil von h'"(z£) mit A, bzw. &,°), so ergibt sich nach } j 
Trennung in Real- und Imaginärteil ; i Yo4 
In 0 (y) = Ar EN EEE er, de (1.4a), h 
1 17) ie 27 5) DE ET I ARTS Po © (1.4b). # 
Um die Beziehung zwischen A; und eg; abzuleiten, wird k®(zf) hingeschrieben, 
N 
Oz) = (am — ib) (COS MpL + isin myy), # 
\ m=0 
und in Real- und Imaginärteil zerlegt, i ee 
n & \- =. 
—- N „(Gm cos MPr + 7, sin MQr) NR REN we DE 
if. | ar (k=1,2,...,2n) RR, 
B; = P2 ‚(m sin moL — bj cos MY) N a Me a RE RR ES Ba 
Dabei ist außer d, =0 auch noch b„ = 0 gesetzt, da sich dieser Koeffizient nicht aus den Io | 


benutzten Gleichungen "ergibt und über ihn frei verfügt werden kann. Nun werden die durch 
harmonische Analyse aus (1.5a) gewonnenen Koeffizienten _ 


1 an ‚ 1 m | ‘ : Wr 
ur. an = SAY ui =0, Be 
i=1 i=1 Re 
De j an \ eg IRRE (1.6) 


1 ; 
Om = Ih cos mon Bm = Ih sin mg 
l=1 


ee ee 
in (1.5b) eingesetzt, 


& -1 In: Bin mn) = Sn er Nee. re R 


Für die Koeffizienten f;, ergibt sich unter Verwendung einer bekannten su 
ann, bei Beachtung der Beziehungen und Bezeichäiiiil von (1.2) 
i-e9r TB 
Pr = ee, eig et Pr=0. 


‘ 


rn 


N 5) Der Index n wird hier fortgelassen, da jetzt der & 
wir 


°) &, hat eine anschauliche Bedeutung, es ist die , W 
den der erde 2; gedreht werden muß, um in den Vektor 
2)—cos[(?r+1)a/2] 
?) 3 sin ma = „re elBrtneB, SR 
2 Nas ei, AR x 


m=1 


5 je! ” EEG, "Min 2077 ir ZU‘ w sp 
a ba a N ES ae En U ER a 
Tan-\ kg EAN nn n 


/ 
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Für gerade Werte von k—1 verschwinden also die fz,, so daß in der Summe (1.7) nur 
n Summanden auftreten. Schreibt man (1.4a) mit (1.4b) in der Form 


4. = In 0(9 + &) Me 2 2 RN EI IE 
< so ergibt sich mit (1.7) 
an 

& = Pr 3". 2 12000 Van ob 


eine Form, die für ein Iterationsverfahren geeignet ist, das zur Berechnung von {A} und {e} 
dienen kann®). Mit einer Anfangsnäherung {e}, beginnend, berechnet man aus (1.9a) {A},, damit 
nach (1.9b) eine neue Näherung {e}, usf., bis die Iteration innerhalb der gewünschten Genauig- 
keit zum Stillstand kommt. Eine hinreichende Bedingung für die Konvergenz des Iterations- 
| verfahrens wird im II. Abschnitt abgeleitet. (1.6) liefert dann die Koeffizienten von h® (z) 
ö und damit die Näherungsfunktion. 


2. Praktische Ergänzung und weitere Formeln. 


Wie man sieht, kann man mit der Durchführung der Rechnung eine geübte, aber mathema- 
tisch wenig ausgebildete Hilfskraft betrauen; bei gegebener Umrißfunktion o(y) ist nur die 
Kenntnis des Formelpaares (1.9) nötig. In rationeller Weise läßt sich die Arbeit nur mit einer 

; Rechenmaschine durchführen, die für (1.9b) die Aufsummierung von Produkten unter Berück- 

En sichtigung ihrer Vorzeichen gestattet; die Schreibarbeit und die damit verbundenen Fehler- 

möglichkeiten werden damit auf ein Minimum reduziert. 

Da die Hälfte aller f,, verschwindet, zerfällt das System (1.9b), für k =1 bis 2n hin- 
geschrieben, in zwei Teilsysteme, denn zur Berechnung der eg; mit geradzahligem Index k werden 
nur die A, mit ungeradem Index / benutzt und umgekehrt. Man behandelt deshalb die beiden 
Teilsysteme getrennt und wendet am besten die Schiebestreifenmethode an. In 
jeder Gleichungszeile treten nämlich die gleichen Koeffizienten f;; auf, da diese nur von der 
Differenz k — 1 abhängen, nur jeweils zyklisch verschoben. Man schreibt die Folge der A, mit 
ungeradem Index vonZ=1bis! =2n — 1 und mit geradem Index von I = 2bis I = 2n getrennt 
untereinander und fertigt sich einen ß-Streifen von doppelter Länge an, auf dem die Werte ß;; 
vond=k—l=2n—1bis d= —(2n — 1)für alle ungeraden Werte von d (also nach fallen- 
den Werten von d) aufgetragen werden. Zur Berechnung von g; legt man den f-Streifen so an 
die Spalte der A, heran, daß für nebeneinanderstehende A; und fz; die Summe von / und d gleich 
dem gewünschten % ist. Dann wird nach (1.9b) die Produktsumme gebildet, ohne daß ein 
Zwischenergebnis festgehalten zu werden braucht. Da je zwei ß;,-Werte sich nur durch das 
Vorzeichen unterscheiden, können die beiden‘ zugehörigen Produkte mit einer Einstellung 
berechnet werden. Steht für die Berechnung von 4; nur eine Tafel der dekadischen Logarithmen 
zur Verfügung oder benötigt man eg; für die weitere Rechnung ausschließlich in Gradmaß, so 
kann man die zugehörigen Umrechnungsfaktoren zu den £ß;; schlagen. 

Infolge der Linearität des Systems (1.9b) braucht man jeweils nur die Verbesserungen 
fd — {etm-ı aus {A} — {Alm—ı zu berechnen und zu {e}„—ı zu addieren, um {e},, zu erhalten. 
So vermeidet man vielstellige Multiplikationen, muß allerdings die Rechnung von vornherein 
mit der vollen benötigten Stellenzahl durchführen und zur letzten Iteration wieder die Werte {2}, 
direkt benutzen, um Abrundungsfehler auszugleichen und eine Rechenkontrolle zu haben. Als 

. ausgezeichnete Kontrolle eines jeden Iterationsschritts erweist sich die Gleichung (3.2), in der 

S U, und V; durch A, und ex zu ersetzen sind; diese Rechnung muß mit großer Stellenzahl durch- 
geführt werden, um die Fehler (absolut genommen) kleiner e;-Werte erkennen zu können. 

_ Wendet man das Verfahren häufiger an, so kann man unter Verzicht auf die bequeme 
Schiebesstreifenmethode durch Zerlegung in geraden und ungeraden Funktionsanteil den Glei- 
chungen (1.9b).eine Gestalt geben, bei der die Zahl der Multiplikationen auf die Hälfte reduziert 
ist?). Dazu ist die Aufstellung einer nur für eine bestimmte Punktzahl 2 geltenden Koeffizienten- 

_ matrix mit n?2/4 bzw. (n— 1)/4 nicht verschwindenden Elementen notwendig. 


Dazu bringt man zunächst (1.9b) auf die Form!®) 


n n—1 
= I Brrht I Ami Aen-ı eo Ro (2.1) 


F 8) Durch die geschwungenen Klammern {...} soll die Gesamtheit aller Werte des in den Klammern 
stehenden Buchstabens bezeichnet werden. C% 
0 9%)H. Wittich, Jb. dtsch. Luftfahrtforseh. 1941, S.I 52—57. 


» Fr Nr j 10) Bei diesen Reehnungen beachte man die für die ßgz und Ar geltenden Rechenregeln: 
gen h Periodizität, 
00.000 H. fu=Pk+2n, 1 = Pr, I+2n i 
DE Pre — Pur Vertauschung der Indizes, 
IWW. Bu=—Pß-u- ' Gleichzeitige Vertauschung der Vorzeichen der Indizes. 
AR 1 25* 
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| oder 
i 5 Bin/nt Pr, 2n An + 2 (Brı Ar+ Pr, —ı 3 n-1): 
Fr Unter Benutzung der Identität | Dr 
e Prı Ar+ Br, —ı Agn—ı = > (Brı+ Pr.) (Ati) + > (Br ı — Br,—ı) (Ar — Aen-ı) RT I 
3 | RE 
Br kann man mit den Abkürzungen NT. # 
= Atom di= h— honı (=1,2...n—1) .. (a0 09 
sp Agn ’ En SIR, Due AR en 
1 1 | OR 
Arı = (+ Br, —ı)» B.ı = 5 (B.ı — Br, —ı) 0,1. ne _ (2,2b}.2 f\ }- 
: ; x { ER 3 
° die Summe in folgender Form schreiben: ke R7 » 
n—i L> - &:2 ben R ‘2 N 
& = S,+ D; = B3 Arı s+ er B; ıd, (k = 0, T, ... n) .” uRenR (2.3a). 2 RN A 
; Die übrigen e;-Werte findet man wegen PR 
Ri; Re Apn-nı= —Arı, Bonn Bar 
Be; aus der Beziehung / ER; vs ri 
En: en =— + Di SF E RI) 223), 
a "ER, $; und D; müssen also getrennt berechnet werden. Wegen Bj; =— Ay braücht nur die ER h 
2: Matrix der A;, berechnet zu werden, die noch dazu wegen Ay-,n-ı =— 4, zum Zentrum 


anlisymmetrisch ist, so daß die oben angegebene Zahl von Koeffizienten übrig bleibt. Zur 
I; Kontrolle kann man die Zeilenprobe | Ren; 


n—1 - > 7 NEN 


At Ant 2 Au =0 (k=12%..,n—1) a 


benutzen. 
Weiter sollen noch einige Formeln mitgeteilt werden, die es gestatten, Real- und Imaginär- 


teil von 5h® (z*), 1(9) und e(p), an beliebigen Punkten p = 8 =mr0 (öl < en an 


TPN unmittelbar aus den EP. zu berechnen, also ohne Kenntnis der Koeffizienten c,„. Die Methode a 
La ist die gleiche, die zur Darstellung der e, führte. In die N für A(9x,5) bzw. &(9,5) — a 
N in (1.5) ist nur @, durch @,,5 zu ersetzen — werden die Fourierkoeffizienten (1.6) eingesetzt. 
Nach Umformung der trigonometrischen Summen erhält man die folgenden Formeln: 


- 


| 4(9r,2) 2 %1,8 x A rn = % E 8, ER 2, zii 
i ae - Be 
mit 3 , BT 


r i : 38.3 
pr Se DE EN ale Bar a 
0 - sin nd clgZ P-n0; RR EHE RR U 30 


 &(P4,0) = 


mit 


Durch die Beziehungen BT “ 


} 


Ak, —ö = 0, R 


Bi die numerische Berechnung erleichtert, un 


NE: Re u Dr 
können zur. Kontrolle dienen. | 
Par: die euer von IR 


Y > an, N . 2 ; 4 ? ” . in 
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2n 
Ale) ERBEN en re (2.6) 
mit ® 
0 tur .k = |, 
= En N eig2 pri für k |; 
en 
€ (Pr) = Bhhı N RER, ZUR SDR (2.7) 
mit 
0 für k—l gerade, k #1, 
FR für k=]1, 
R 2 
Prı = N 


— für k— 1 ungerade. 
2n- sin?> P-1 


Setzt man A(p) = const., so ergibt sich die Kontrollformel 


an 


Ser 0. 


le1 
(2.6) und (2. 7) kann man eine der Gleichung (2.3) analoge Form geben (siehe H. Wittich, 
a.a.O.). Auf die unter Umständen schlechte Annäherung der Ableitung einer in den Funktions- 
werten relativ gut angenäherten Funktion sei nur las 


x ® 
\ II. Die Konvergenz des Verfahrens 
| 3. Ein Hilfssatz 
i Grundlage der Konvergenzbeweise ist folgender Hilfssatz: (3.A) 
; Besteht zwischen zwei N von je 2n Zahlen {U} und {V} die Beziehung 
N Ar RN NS N, 
n wobei die ß;ı durch (1.8) definiert nu so Ex 
BR 2n 3 


Dabei bedeutet > und 3, daß über alle BA bzw. es Werte von k=1,2,...,2n 
(u) (0) 
summiert werden soll. Zum Beweis stellt man durch harmonische Analyse das Ei ienieteiiele 


; 2 Polynom n-ter An auf, das an den Stellen @; die Werte U, annimmt, 
= 2 (An cos mp + Bm sin m pr), EEE EEE AR (3.32). 


Dann genügt {in folgender Beziehung, in der ee gleichen Koeffizienten An Bm wie in (3.3a) 
BREITER 


=>; (An sin mr — Em COSmMp) 2 2 22er +@836). 


- m 


Man sieht das’ sofort, wenn man (3.3) mit: d. 5) vergleicht: 
{U} und {V} stehen im gleichen Verhältnis zueinander wie {A} und {e}. Da nun (1.9b) 
2 aus (1. 5) folgt, folgt auch (3.1) aus (3.3), und wegen der Eindeutigkeit von (3.1) kann es kein 
ee System {V} geben als das nach (3. 3b)  darstellbare. 

Für {U} und {V} und die Koeffizienten Am B m ihrer trigonometrischen Polynome gelten 
i hun die Ba aengen?) a... iu 


AR 


i , 2. Bu Berlin 1950, $. 222, 


al N v2 IE 


% gestatten. Die Rechnungen verlaufen wie der in der oben beschriebenen Weise!!) und liefern: 


erden als die Ableitungen anderer Summen ge- 


rn 
i 
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$, Es ist also zunächst | s 

Bi , ZN= 8 U—2n(di+ A) 

Z h Dem letzten Term gibt man > andere Gestalt,- 

} 4 | 2(A5+ Ar) = (As — An)? +(A,+ An)? 

F & nun ist entsprechend (1.6) | R 

Br. I m Aid: Art nm 2 Ui. 


Setzt man diese Werte in (3.5) ein, so erhält man aus & 4) die behauptete Gleichung (3. 2), aus 
der sofort die Ungleichung 


2 an 2n 

j ee 

f k=1 k=1 

3 folgt, so daß der Hilfssatz kurz in folgender Form formuliert werden soll: Er 

ä Aus (3.1) folgt (3.6). | (3.B) | 


4. Die Konvergenzdes Iterationsverfahrens 
\ Als hinreichend für die Konvergenz des Iterationsverfahrens gemäß (1. 9) erweist sich die 
B: . Bedingung °?) 
A M = Max e v2 
. je) 
Anschaulich gesprochen bedeutet dies, daß die Normale des Randes R mit dem "Fahr- 
strahl vom Nullpunkt nur Winkel bilden darf, die kleiner als 45° sind. 


ERBE ERERN ER | 


N Für den Beweis wird zunächst (1.9b) für einen beliebigen m-ten Iterationsschritt hin- 
geschrieben!*) und (1,9a) en 
ak >. 
2n 5 1 
&,m = az Prı Am == Bu ! Ino(p+ &,m—ı) LEARN ER (4.2). 
ze 


Jetzt denkt man sich die gleiche Gleichung für den nächsten Iterationsschritt darüber geschrieben 
(m wird durch m +1 ersetzt) und bildet die Differenz 


he: ä 2n N PS $ en 
E | Eu mtı — &hım — 2 Pi Ine@+Eam)—meolpitema)]) -.... A 
Wegen der Stetigkeit von o(y) kann man weiter schreiben. - NS £ j 

2n 1’? 


Ehmt1 T ähm = Be Dahn — &m—ı) re ehe . (44), 


. wobei D, m eine Zahl ist, die zwischen dem größten und kleinsten Wert Mech den [0 We 
im Intervall (Pr Em» Pit &,m—ı) annimmt. Nach @ a ist nun 


> (enm+ı — Em)? <S Di,m Be en 
oder, wenn man D},„ durch den wegen (4.1) sicher nicht kleineren Wert Mm: erset 
2 (um — m)? < m (nm — em)? 


Bezeichnet man den In dinienkeöihlen Vektor {dm vorühergehend 


(4.5) in der Notation der Vektorrechnung Ph: =. 


23 


lem — em < Mm Y | 


‘ 12) Siche die Berähe Bedingung bei) G.Opitz, 
der Satz „Wegen .... &(p).“ unter GI. ee: „D: ? 
Maße Null gegen s(p).* „ 
1.77...) Die ge durch | K 
VER Ki ee) 


und daraus schließlich nach einigen Umformungen 


Mm 
[em+» — eml<; Tas Cr (4.7) ; 


für beliebiges ganzzahliges p > 0. Da die rechte Seite wegen (4.1) beliebig klein gemacht werden 
kann, indem man m genügend groB wählt, bedeutet dies aber gerade, daß e,„ gegen einen Grenz- 
vektor e konvergiert, 
une al50 DEE. ee (4.8). 
MX Mm © ® j 
Daß dieses Grenzsystem {e} tatsächlich ein Lösungssystem darstellt, sieht man, wenn man 
in (4.2) m gegen Unendlich konvergieren läßt. Rechts kann man wegen der Stetigkeit von o(y) 
Funktions- und Grenzwertbildung vertauschen und gelangt so mit (4.8) zu (1.9). Damit ist die 
Konvergenz des Iterationsverfahrens und die Existenz einer Lösung nachgewiesen, und zwar 
unabhängig von der Anfangsnäherung {e},. 


Die Eindeutigkeit der Lösung wird leicht indirekt bewiesen. Gäbe es zwei verschiedene 
Lösungssysteme {e}' und {e}YP),so müßte für beide (1.9) gelten. Durch Subtraktion erhält man 
ei — el — B: Ba Inolp + eP)— no tel?) »...,. (4.9), 


und behandelt man diese Gleichung wie (4.3), so ergibt sich entsprechend (4.5) 


; 


2n 


= (ei gu we <M: = (ec) [5 euUmy2, 


also, da die Quadratsumme wegen er vorausgesetzten Verschiedenheit beider Systeme nicht 
verschwindet, M 21 im Gegensatz zu (4.1). Es gibt demnach nur eine Lösung. i 


Bu 


5: Einige Eigenschaften der exakten Abbildungsfunktion 
Setzt man in Theodorsens Ansatz für die Abbildungsfunktion 


Er Bee en ee (5.1) i 

u, . m=0 2 

= h(2*) = (9) +3 &(p) (esist 2* = eiP), so gilt für die Winkeldrehung auf dem Rande die Integral- 

= gleichung (s. Theodorsen, a.a.0.) wi. 

we j PR ; ; 

} x me ae PREERN 17 

&(p) | ip): ci dp NER GREEN (5.2) 

0 x 4 

Br mit | Er „% }. 
Br. = | A(g) = In az ae ER (5.2b). 


Unter Berücksichtigung der Periodizität von A(p) läßt sich das uneigentliche Integral 
(5. 2a), von dem der Gauchysche Hauptwert zu nehmen ist, umformen in 


7 


1 
(=>; e Badener. (5,3), 
3 ; j 0 
BARON wobei gesetzt wurde 
E ER Ä ? u .0(p, 0) = Rp — 0) —Alp+ 0 )]- cig5 0 Re, (5.4). 


er 8; 0) ist für « = 0 nicht definiert; durch geeignete Definition läßt sich jedoch diese Un- 
ns ste igkeit beheben. Dazu müssen noch einige Zwischenüberlegungen angestellt werden. 


® Setzt man in (5.1) für z die Randpunkte 2* = e'? ein, so ergibt sich nach Differentiation 
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Dies wird in (5.5) eingesetzt, und man erhält 
%.p) = a N En RC 
epPriep) SCH 


In allen Randpunkiin, in deren Bildpunkten der Rand R glatt verläuft, also o’(y) stetig 
ist, ist dies auch N (z*) und damit A’(p). Besondere Untersuchung erfordern nur die Unstetig- 
keitsstellen von o’(y), denen Soken. im Rand R entsprechen. Ist 5 = f(zö)ein Eckpunkt, so läßt _ 
sich f(z) in der Form 


on 


w— uw, = (% — 2) VPl&— 28). 7. x VER (5.8) Ei 


darstellen; P(2) ist eine Potenzreihe in 2, deren Absolutglied P, nicht verschwindet, 0 270 ISE En A 

der Winkel, vom Innern von @ aus gerechnet, den der Rand im Punkt w5 bildet. Damit wird e 

Ä (3.7) zu 3 

Zn) Em ‚net — a) Per —ed)4 A): Pa) ER 
4 F Z : Det 1 er L 

I N e'y)+te(y) ws + (2* — 28)*: Pe* — 20) \ % 

F Man sieht, daß A’(p) für 2*— z$ über alle Grenzen wächst, wenn « < 1 ist, d.h. wenn @ ch 4 

Ei ausspringende Ecken besitzt. Ausspringende Ecken sollen deshalbin Zu- 

Br kunftausgeschlossenwerden. Für«a >], also einspringende Ecken, geht dagegen 

x 

h 


A X (p) für 2* > 28 stetig gegen Null, da die Unstetigkeit von e’(y) in diesem Fall unwirksam wird. B 
Fi Man hat so das Ergebnis: 5 
N, $ (p) ist im gesamten Intervall 0 <p <2n stetig... . 2.2.22... (5.A) 


Da ®’(p) gleichfalls in einem abgeschlossenen Intervall stetig ist, ist es. im gesamten 
Intervall 0 Sp < 2 beschränkt. 


N. Daraus folgt weiter sofort: 
En ©(p) genügt einer Lipschitz-Bedingung ..... re . (5.B) 


IR Nach Bereitstellung der nötigen Sätze wird die On im hu an 1.6. 4) fort- | 
Br. gesetzt. Wegen (5.A) kann man (5.4) dem Mittelwertsatz der Differentialrechnung zufolge De 
Be schreiben RN 


e j 99; oc) = —27 (p+ 70): acigs o, en le BEERE EL RR. (5.9). 


Ist nun C ein beliebiger (nicht notwendig analytischer) Weg in dem Bereich k Y f 

ar 0sps2n, 0<o<a, der in dem Randpunkt P,(g', 0) endet, so wird, wenn man a 
den Punkt P(9, o) längs C’ nach P, rücken läßt, wegen (5.A) ©(p, o) stets dem gleichen Grenz- ir 
af? wert zustreben, Er 


* Durch diese Definition der Randwerte © (9, 0) ist 9p, 0) im ganzen Gebiet 0sp = RE % 
0 <So<a stetig, und daraus folgt noch schärfer a 


- .°.0(9, 0) ist im abgeschlossenen Bereich 0 < 9S2n,0<So<a gleichmäßig stetig (5 .0). ER 


Für den folgenden Beweis soll weiter der Begriff des Stetigkeits- oder Schwaikuneeh a 
moduls°) eingeführt werden. Betrachtet man ©(9, o) bei: festgehaltenem 9 als Fauknope von Rh 
so ist der Stetigkeitsmodul &(p, 6) definiert durch I 

w(9, 6) = Max |9(9, 0,) — 9(9, 0,)| für alle |, — | <ö . Fr en 1. 


die im Intervall 0So,, ,<r liegen. Wegen (5 RN ist (9, gr eine mit wachsendem * ii 
ö nicht abnehmende Funktion von 6, für die silt 2 HINTER 


im 9 )= 0, AR) U RE RE 6.10). d 
| 
h 


gleichmäßig für alle o. 
o(9, 6) Aunvergiert mit 6 gleichmäßig für alle 9» gege 


6. Die Konvergenz der Nä hi 
Um die exakte Abbildungsfunktion mit der 


ö en auch an den diskreten Punkten Pr beta ghtet, 2 
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valle mit den Teilpunkten 99. 91, : - -, 9n undnähert das Integral (5.3) nach der Tangenten- 
trapezregel!*) an; das Restglied wird mit Hilfe des Stetigkeitsmoduls abgeschätzt. Unter Berück- 


sichtigung der Periodizität erhält man mit A(gx) = Ar. E(9) = & 
/ 3 RT TEN 1 ee 
€ = I Puhtzy (Pr): Iyr| el ee Na a Et Se (6.1). 


$ Da jetzt der Grad der Näherung veränderlich ist, wird er zu A, und e, als oberer Index 
e; zugefügt; damit wird (1.9a) 
9% lo +E). 


Nun wird die rechte Seite als Funktion von e{”) an der Stelle &, entwickelt, 
ie — no +&)+ EN ( — 5) = + FM. (& — &); 
; Fi ist ein Wert von o’(y)/o(y) im Intervall (+ &. 9 + &). Obige a wird nun in 
a. 9b) eingesetzt, 
an 


ei) = =2 ut BrıPE (ei — 5), 
und davon wird (6.1) abgezogen, 


F {n) lausr pin. Ne N LO NR 2 
ET. = Bu I lei —r 8) z% o(P,, ni 
Schreibt man das letzte Glied auf die-linke Seite und wendet (3.B) an, 
5 (n) = I 712 SR) (n) ZANZ 
Dale &)+ 5 Ye’ @(Pr -)E<23 [Fr (9 — &)]°; 
k=1 n k=1 
so kann man rechts F{”2 wegen (4.1) durch M? ersetzen und links nach Ausmultiplizieren den 
einen quadratischen Term fortlassen, ohne die Ungleichung zu verletzen. Man erhält so 
an. % an er Kr 
(1— M?) 2 (rl &)® = 2 KL: "© (9%, (er en (0:2) 
Da die kleinere Seite sicher positiv ist, kann die ganze Ungleichung quadriert werden. 
Dann wird auf die rechte Seite die Schwarzsche Ungleichung angewendet, 
2n 
W (1 — an): 3 (en. . &x)? = = YE' (Pi )- 
2 Nun ersetzt man rechts y; durch 1 und 'o(9;,r/n) durch den größten Wert la: 
R den »(g, /n) im Intervall O Sp <2r annimmt; die Existenz von w, ist durch (5.0) gesichert. 
i 1 en 2 2n j 
| \ oe Se De 1250. 
2 ' 


Hier läßt man nun n gegen Unendlich konvergieren und erhält mit (5.D) 


lim $, lim -, N UT N RE a (OA); 


Nn—x 


Bsi2ie mittlere quadratische Abweichung von ec) (9) von &(p), genommen an den Punkten 9;, 
_ konvergiert also gegen Null. Das Bemerkenswerte an diesem Grenzübergang ist, daß sowohl die 
Br Funktion wie auch die Punktzahl sich mit n ändert. Zum Nachweis der mittleren Konvergenz 
Re n yon e”)(p) gegen &(p) ist deshalb die Einführung einer Folge von Hilfsfunktionen &(p) not- 
vendig. Ba) ist das az tieche Polynom Arten Grades 


, zn) = Fun sin mp — Arm cos MY), on = mn =0. .. . . . (65), 


er; 
Pr} 


ir den 2m Punkten 9% übereinstimmt, 

FRER ‚a fi : & Em) < RN Et = En. 
ı \ icth. der prakt, Aal 24 or es Berlin. 1950, 8. 136. 
M AR 


a kon KR ARE Pac a ZA 
a N ee 2 


h 
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Damit nimmt (6.4) die Gestalt » $ Bu 7 # 
EN Fe 

lim S, = lim 1 Siam _ em) = NT a INC n i 

en no N ARE ER. 

BR 


an. Aus (1.5b)und (6.5) folgt, daß auch die Differenz &)(p) — E®(p) vom n-ten Grade ist, 


Er (9) — E® (9) = 3 (Amn Sin MP — Bmn COS MP), fon =&nn=0 . : . . (6.7). 
m=0 


’ 


Da e®(g) das konstante Glied und die n-te Harmonische nicht enthält, ist fon = Bon 
und f,n = Ban. Für (6.7) gilt nun (s. Willers, a.a. O.:2)). 


an n—1 N 


—_ -13 (ei) — Em)? — 2 (Bön+ Pant 2, (mn Ban) See (6.8) 


reg 


und andererseits!?) 


er z Ba! z 
Pr ee ee io 
ge Ba er ee ARTEN TE | 


PERS % 


Fr | ni 
= a [ Le) — EI dp = 2 nt Al + nt) 3 (09) 


i so daß man die in (6.6) auftretende Summe durch ein Integral ersetzen kann, "rs 
B | RENT HRERE NN ER) | 
Ei } 4 LI 
u Durch Aufspaltung wird aus J„ "ö 
Br 27 , X y : 
‚A In = [ {Le (@) — & (P) — E@) — E(p)}}?dp 
er‘ d ! ’ 


und durch Umformung mit Hilfe der Schwarzschen Ungleichung 


2 [Vier =unra- fern= anni] em. 
Mit (6.10) wird daraus Re ä | 
vr [Le (9) — ep)? dp — [tm (9) — Ep)? an) 


Für die einzelnen Teile dieser Ungleichung gilt mit wachsendem n: 
a) 8, geht nach (6.6) gegen Null, ar 

b) Ban = Ban geht wegen der Stetigkeit von &(Y) gegen Null!8). ni "1 
c) Da &®(@) wegen (5.B) gleichmäßig gegen €(p) Kante sb: der zweite e Term R ” 


Yr (Sn vn Ban) = 


en der rechten Seite gegen Null. \ 
00%. Daraus folgt, daß auch der erste Term der rhlen Seite gegen Null seht, 
A 1 [ten (9) — ep?’ dp = % RN Re Ye ‚12 
d.h. die Folge der mg) a im Mittel gegen 2(9) Ne > a A 


Damit ist der Beweis abgeschlossen. Es seien nochmals kurz die Bedingungen ange; 
unter denen er geführt wurde: Das Gebiet @ ist ein Sternbereich, dessen. and m 4 
nahme einspringender Ecken eine sich stetig drehende Tangente besitzt. Auß 
3. (4.1) auf „‚kreisnahe‘‘ Gebiete begrenzt. Die nz der Krümmung ® 
vorausgesetzt. SR ß 


1?) Siehe z. B. Fr. A. Willers, Meth. der pra 
18) W. Quade, Deutsche Mathematik 5 (1 
19) W, Quade, "Deutsche Mäthematik 5 


Kt _ Eingegangen am 4.10. 1999. 
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Dreidimensionale Probleme der Charakteristikentheorie 
partieller Differentialgleichungen 
Von R. Sauer in München 


Die Arbeit handelt von 3-dimensionalen Problemen der C'harakteristikentheorie partieller Differential- 
gleichungen. Im ersten Teil werden Probleme erörtert, die sich infolge gewisser Symmetrien auf 2- oder 
I-dimensionale Aufgaben entweder sireng oder in „halblinearer Näherung‘ zurückführen lassen. Im zweiten 
Teil wird eine Gitterkonstruktion entwickelt, welche dem bekannten Massauschen Verfahren für 2-dimen- 
sionale Aufgaben analog ist. 


The author investigates 3-dimensional problems in ihe theory of characteristics of partial differential 
equations. In the first part, problems are treated ihat because of certain symmetries can be reduced to 1- 
or 2-dimensional problems, either exactly or by a „‚semilinear approwimation‘‘. In the second part, a con- 
struction of lattices is developed that is analogous to ihe well-known Massaumethod for 2-dimensional 
problems. 


Dans la communication pre£sente, l’auteur Etudie des probl&mes a 3 dimensions de la th£orie des caracieri- 
stiques des Equations aux derivees partielle. La premiere partie comprend des probl&mes qui, en conse£- 
quence de cerlaines sym£iries, se reduisent a des probl&mes a 1 ou 2 dimensions, soit rigoureusement soi 
dans une „approximalion semi-lineaire‘‘. Dans la seconde partie, il est d&velopp£& une construction de 
grilles analogue a la methode de Massau bien connue pour les probl&mes ü 2 dimensions. 

B aannofi pa6ore uecnenyerTcH Tpex MepHas sanaya Teopun xapakTepucTuk AndepenH- 
IIHAJIBHEIX YPABHeHHÄ B YACTHBIX IPOHSBONHEIX. B epBoä yacTu paccMaTpuBamTcH 3anayn, 
KOTOPEI® BCJIEICTBHe HEKOTOPEIX CBOÜCTB CHMMETPHH MOTYT ÖBITB TOYHO HIIH B „‚HOMY-JIH- 
HeÄHOM‘‘ IIPHÖINKeHNÜ CBeNeHEI K 3anayaM ONHOIO HIN ABYX usmepennü. Bo BTOopo& yacru 
paspadaTsIBaeTcA MmeTon Pelle ToK, AHAJOTHYHUBIH H3BECTHOMYMETOAY MACCO MIA INIOCKUX 3aay. 


1. Einleitung 


Zur numerischen Lösung von Anfangswertproblemen partieller Differentialgleichungen, 
wie sie z.B. in der Gasdynamik und in der Strömungslehre seichter Gewässer auftreten, ist bei 
2-dimensionalen Problemen das Massausche Charakteristikenverfahren ein sehr brauchbares 
Hilfsmittel und wurde in zahlreichen Arbeiten [1] [2] [3] für derartige Aufgaben verwandt. 

Neuerdings ist das Bedürfnis nach geeigneten Charakteristikenverfahren für 3-dimensionale 
Probleme immer mehr in den Vordergrund getreten. Hierbei ergeben sich 2 Möglichkeiten, 
nämlich 1. die Reduktion 3-dimensionaler Aufgaben auf 2- oder 1-dimensionale und 2. die Er- 
weiterung des Massauschen Verfahrens auf 3 Dimensionen unter Vermeidung eines praktisch 
nicht mehr tragbaren Arbeitsaufwandes. Beide Möglichkeiten sind mehrfach versucht worden. 

Die Reduktion auf geringere Dimensionszahl durch Drehsymmetrie oder Kegelsymmetrie: 
führt bei Lösungen dreidimensionaler Probleme, die nicht allzu stark von dreh- oder kegelsymme- 
trischen Lösungen abweichen (z. B. Strömung um Drehkörper bei kleinem Anstellwinkel), zu 
„‚halblinearen‘‘ Näherungsverfahren. Bei diesen wird die reduzierte Aufgabe entweder mit dem 
üblichen 2-dimensionalen Charakteristikenverfahren behandelt oder sogar noch weiter auf eine 
1-dimensionale Aufgabe gewöhnlicher Differentialgleichungen reduziert [4] [5] [6]. Ziffer 2 gibt 
einen systematischen Überblick über diese Reduktionsverfahren. 

Die Erweiterung des 2-dimensionalen Massauschen Charakteristikenverfährens zu einer 
3-dimensionalen Gitterkonstruktion bildet den Hauptgegenstand der vorliegenden Arbeit und 
wird in Ziffer 3 ausführlich behandelt. Für spezielle Fälle wurden solche Gitterkonstruktionen 
schon mehrfach entwickelt, erstmals wohl von C. Ferrari [7] und neuerdings insbesondere in 
einem Buch von A. Ferri [8] über stationäre Überschallströmungen. Es erscheint jedoch nützlich 
diese Gitterkonstruktionen unter allgemeinem Gesichtspunkt systematisch darzustellen und zu 
ihrer praktischen Anwendung eine Gebrauchsanweisung anzugeben, wie dies in ähnlicher Weise 
für 2-dimensionale Charakteristikenverfahren bereits früher von K. Oswatitsch [9] und vom 
Verfasser dieser Arbeit [10] geschehen ist. 


2. Reduktion auf 2- und 1-dimensionale Probleme 
Wir-erläutern die Reduktion 3-dimensionaler Probleme auf 2- und 1-dimensionale am Bei- 
spiel der Überschallströmung idealer Gase (# = cy/c, = const = Verhältnis der spezifischen 
Wärmen bei festem Druck bzw. Volumen). Sie genügt den 4 quasilinearen Differentialgleichun- 
gen erster Ordnung 
(uw — a?)u, + uvu, + uwu, + vuv, + (0 — a?)vy + vww, + wuu, + uwvuy, 
+ (w — a)w, —=0, 
a? 
v(u, — de) + Ww (u, — W,) N2@&— 1) ,=0, 
w(%, — w,) + um, — u,) 1) =0; 


us, + vs, +u,=0 


(ı). 
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für die 4 Funktionen , v, w (Komponenten des Geschwindigkeitsvektors) und s ‚ (Entropie) SS 
3 cartesischen Koordinaten x, %, z. Die Schallgeschwindigkeit a ist durch 


(u +v? + w2), Ag = Consb N 

IN 2 
als Funktion des Geschwindigkeitsbetrags gegeben. Die Untersuchung beschränkt sich auf 
Bereiche mit Überschallgeschwindigkeit d.h. mit w+v? uw? >a?. _ 


“—1 
| 
?=m— 


2.1. Drehsymmetrische und kegelsymmetrische Lösungen 


Zur Ermittlung drehsymmetrischer Lösungen (z-Achse = Symmetrieachse) führt man 
Zylinderkoordinaten o,r,2 (Bild 1)ein und erhält dann aus den Gln. (1)ein System von Differen- 
tialgleichungen mit nur 2 unabhängigen Veränderlichen z, r, also ein 2-dimensionales Problem. 

Da die Charakteristiken dieses Problems im Überschallbereich reell sind, läßt sich das 
Massausche Charakteristikenverfahren verwenden [1]. 


Bildi. Zylinderkoordinaten o, r, 2 und Geschwindigkeits- Bild2. Polarkoordinaten R, #, ® und Geschwindigkeits- 
komponenten %, vi, © komponenten U, PV, 


Zur Ermittlung kegelsymmetrischer Lösungen (4, v, w, s— const auf jeder Halbgeraden 
durch den Symmetriepunkt <= y=z=0) führt man Polarkoordinaten R, d, ® (Bild 2) ein 
und erhält dann aus den Gln. (1)ein System von Differentialgleichungen mit nur 2 unabhängigen . 
Veränderlichen d, », also wieder ein 2-dimensionales Problem. MEN 


Die Charakteristiken dieses Problems sind jedoch nur in dem Bereich reell, in dem die zu 
den Geraden durch den Symmetriepunkt senkrechte Geschwindigkeitskomponente größer als 
die Schallgeschwindigkeitist (vgl. Ziffer 3, 4). Man kann daher dasMassausche Charakteristiken- 
verfahren i. a. nicht verwenden und ist, falls lineare Näherungsverfahren [14] [15] [16] nicht aus- % 
reichen, auf eine 3-dimensionale Gitterkonstruktion angewiesen. FR 


Bei Strömungen, die gleichzeitig dreh- und kegelsymmetrisch sind, bleibt nach Einführung 

von Polarkoordinaten nur die eine unabhängige Veränderliche 9 und man erhält an Stelle der 
„ partiellen Differentialgleichungen gewöhnliche Differentialgleichungen, also ein 1-dimensionales MN; 

Problem. Auf diese Weise haben Busemann [11], Bourquardt [12]und Taylor-Maceoll[13] 
die Überschallströmung um den Drehkegel mit Anstellwinke] Null und anliegender kegelförmiger ER 
Kopfwelle berechnet. ER 


2.2. Halblineare Verfahren für Nachbarlösungen drehsymmetrischer und kegelsymmeirischer Er ; 
Lösungen at > 


Bei Strömungen, welche von einer drehsymmetrischen Strömung mit der En als Dreh- 
achse nur wenig abweichen und zur Ebene & = O0 symmetrisch sind (Beispiel: Überschallströmung 
um einen Drehkörper bei kleinem Anstellwinke]l e), machen wir in einem Gallen 
system (Bild 1) für die Geschwindigkeitskomponenten den Ansatz [4] ERS 3 


ae EU, (2, r)sin ©, 
v—=v(2,r) + evılz, r)cos w, 
©=%y(e, r) + sw(zr)cosw, [' * 
s=H(&r)+e lan: co | 
Die Falltiollen v%(% 7), Wo(2, r), So (2 nl zusammen N 


eine drehsymmetrische ee und ae 4 sich Zi 
Nee Super 1 “ EN 


1} 
m 


ah > Ba ea Baal Ba BR N 
ee ar ua LEN Mg j 


NER 
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Der drehsymmetrischen. Strömung wird in dem Ansatz (3), in dem die Veränderliche & 
separiert ist, eine durch die Funktionen 


ur), vr), w(% Tr), 81(2, r) 
gekennzeichnete Zusatzströmung überlagert. Wir berücksichtigen e nur linear und erhalten dann 
für die gesuchten Funktionen %,,%,, %,s; aus den Gin. (1) ein System linearer Differential- 
gleichungen mit den 2 unabhängigen Veränderlichen z, r, welches die Charakteristiken der dreh- 
symmetrischen Lösung als Charakteristiken hat und wieder mit dem 2-dimensionalen Charakteri- 
stikenverfahren behandelt werden kann. 

Der spezielle Ansatz (3) mit der Abspaltung der Veränderlichen & wird dadurch gerecht- 
fertigt, daß sich mit ihm die Randbedingungen des Problems erfüllen lassen. Wir bezeichnen das 
Verfahren als ‚„‚halblinear“, da die drehsymmetrische Strömung streng, die überlagerte Strömung 
jedoch mit Beschränkung auf die 
in e linearen Glieder aus linearen 
Differentialgleichungen berechnet 
wird. 

Als Beispiel einer Reduktion 
auf 1 Dimension d.h. auf gewöhn- 
liche Differentialgleichungen bespre- 
chen wir die Überschallströmungen, 
die von einer zugleich drehsym- 
metrischen und konischen Strömung 
nur wenig abweichen und wieder zur 
Ebene ® = O symmetrisch sind (Über- 
schallströmung in der Umgebung 
der Spitze 0 eines spindelartigen 


$ & & Bild3. Überschallströmung um einen nahezu konischen Drehkörper 
Drehkörpers unter kleinem Anstell- unter kleinem Anstellwinkel e 


winkel, Bild 3). 
Wir machen jetzt in einem räumlichen Polarkoordinatensystem N 2) für die Geschwin- 
digkeitskomponenten den Ansatz 


U = U,(9) + & U,(9) cos ® + R U,(%) , 
V’=V,(d) + eV, d)cos®+ RV;(d), 
W= eW,(#)sin o, 
s = -+esıcs® + Rs(d), so=const, s7—=const. 
Durch diesen Ansatz sind die 3 unabhängigen Veränderlichen R, d, &separiert und man hat 
lediglich Funktionen der einen Veränderlichen % zu ermitteln. V 0(9), Vo(9) und sy = const 
liefern eine drehsymmetrische und kegelsymmetrische Strömung (Strömung um den nicht an- 
gestellten Drehkegel), U,(9), V,(9), W,(®) und s,; = const eine durch den Anstellwinkel s hervor- 
2 gerufene Zusatzströmung und U,(9), v0), s;(9) eine durch die Krümmung des Drehkörper- 
meridians an der Spitze 0 bedingte zweite Zusatzströmung. 


2 d 
Wir linearisieren sowohl hinsichtlich & als auch hinsichtlich R und erhalten ( = =) aus 
den Differentialgleichungen (1) nach elementaren Rechnungen 
a) für die drehsymmetrische Grundströmung: 


(— U2)(U+ U) +RUV+V,ctgd)—0, n=U, 


BY aa mit a5 (U24 U) ik 
E b) für die erste Zusatzströmung: 
UN HFRUHWU ALTO, = B 
% Er sin 3 (5b) 
is U, W,+ (U, 7 %% ctg 9) Ww, em UA ıt+ U, U, 1) a Ber ruagar>3: A 1) sr 0 
mit den Abkürzungen 
I Pas 0+ en alu’ + «—1)U8), 


Br U,+U;ctgd 


9) 20 — Vene dir: 


£ 


e 


TERN RT ZT ARTNET, 4 BIETER, 17 
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ce) für die zweite Zusatzströmung: 


U 


je RR, 3 
ES IBER E a]. 3 en 3 


A, 
1 ee ER M 
(@ — V,)V,+LW)U,;,+M) un U,a?s, i 
mit den Abkürzungen | 2 
5 MN 
9) 30 — (U34 VyLR It N), 
BB; 
M)=-—U, + (UF Voctg nn + 0 


v; 

Die nichtlineare Differentialgleichung zweiter Ordnung (5a) für U,(®) liefert die bekannte 
Busemannsche drehsymmetrische Kegelströmung [11] [12] [13]. Wenn diese ermittelt ist, y 
liefern die beiden linearen Differentialgleichungen zweiter und erster Ordnung (5b) für U,(9) 
und W;(®) die erste Zusatzströmung [4] [5] [6]. Gleichzeitig ergibt sich die zweite Zusatz- 
strömung aus den Gln. (5c); die erste dieser Gleichungen liefert s,(), die beiden übrigen Glei- 
chungen sind lineare Differentialgleichungen erster Ordnung für U,(d) und V,(9). 

Die bei den Integrationen auftretenden Integrationskonstanten sowie die im Ansatz (4) 
vorkommenden Konstanten lassen sich eindeutig festlegen durch die Randbedingungen an dem 
von der Körperspitze ausgehenden Verdichtungsstoß. Für die Durchführung im einzelnen sei 
auf die Originalarbeit [4] verwiesen, in der allerdings die Untersuchung der Kürze halber auf 
isentropische Strömungen (s} = 0) und auf kegelförmige Drehkörper beschränkt war. 


3. Gitterkonstruktion für 3-dimensionale Probleme j 

Wenn die für eine Reduktion erforderlichen Symmetrieeigenschaften nicht vorhanden sind, i 
also ein im eigentlichen Sinne 3-dimensionales Problem vorliegt, stehen für geringere, in vielen ü 
praktischen Fällen ausreichende Genauigkeitsansprüche lineare Näherungsverfahren [15] [16] i 
zur Verfügung. Bei höheren Genauigkeitsansprüchen ist man jedoch auf die 3-dimensionale 4 
Charakteristikentheorie angewiesen. Wir werden im folgenden zuerst die benötigten theoretischen pr 
Grundbegriffe dieser Theorie erläutern und dann die durch Verallgemeinerung der 2-dimensionalen Ö 
Verfahren [9] [10] sich ergebende räumliche Gitterkonstruktion behandeln. f 


3.1. Theoretische Grundbegriffe 


Wir erläutern die Theorie an dem System zweier quasilinearer Ditferentialgleichungen 


erster Ordnung 
Ad F AU + AU + didı + Biete + b13%3 = 91» (6) 
Ay + Ag, + Agua + da1Vı + Dead + da3V3 = 9 Ba 


für die 2 gesuchten Funktionen u,» der 3 unabhängigen Veränderlichen z, ys2. Bei « und v be- 
deuten die Indizes 1, 2, 3 Ableitungen nach x, y, 2. Die Koeffizienten a;,, bi, und g%=1,2; Ä 
; k=1,2,3) sind Funktionen von #, 
y, zund von u v. [ar 
Vorgegeben sei eine Fläche 
z2=f(«,y) und eine Werteverteilung 
u,v auf ihr. Hiermit sind auch die 
„inneren Ableitungen“ (Bild 4) 


0) 0 ö 
Fiir EEE 
5 oe 
TE 


_ derFunktionen u,v,d.h.dieAbleitun- 
gen der Funktionswerte u(&, y, f(x, Y 2 ME 
und v(z,y,f(z,y)) auf den Schni a | 

kurven /),m der Ebenen y=c t 
bzw. = Eng Re 


Bild4. Fläche mit Schnittkurven Z (y = const) und m (x = const) 


6) ° 
a y und qy= . bestimmen den Normalenvektor (2, 9; 2) 
Durch passende Wahl des Koordinatensystems jäßt sich stets e rreicl 


3 Er (Bild 8). 
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Fläche keine zur z-Achse parallele T angentenebene besitzt und daher durch z=f(«, y) dargestellt 
werden kann. Der Ansatz z— f(x, y) bedeutet also keine wesentliche Einschränkung. 


Wir versuchen nun zu den vorgegebenen inneren Ableitungen Re SEEN aus den 
dx" öy’ dw’ öy 


Differentialgleichungen (6) die aus der vorgegebenen Fläche herausführenden Ableitungen 
0) 

U= Se Dr — zu berechnen. Zu diesem Zweck schreiben wir die Differentialgleichungen (6) 

in folgender Form 


öu du 6% 
(AıP+ 9— 413) Wa + (di P + b127— bis) 9% = 2 Ei ng, ar n 2 bie Er 91» 
2 6% ÖV (7) 


(@1P + 4 — Ay) Ug+ (21 P + d22 9 — b>s) a ke 37 sr ünz, +57, + 6257 y — 9 


und erhalten dabei folgende Alternagive: j 


a) Wenn die Ausgangsfläche 2 = f(x, y) und die auf ihr vorgegebenen Werte u,» der 
Richtungsbedingung genügen 


” AHıP-+ 29 — Us b1P + b129— bis E, 
> Agı PT Qs2Q — Qsz bp + ba, a — bes 


Be 


haben die für, «;, v; linearen Gleichungen (7) nur dann Lösungen, wenn auch die Verträglich- 
keitsbedingung 


du 6V 6% 
MP 429 — 43 u ta 125, d15, Day, 


91 


= EEE, 
ou 077 ÖV ÖV . (9) 


UP Ag — Ans 2215, | 5, td 5, Hbaz, ga 


erfüllt wird; die gesuchten Ableitungen sind dann aber nicht eindeutig festgelegt, sondern es 
gibt unendlich viele Lösungen “3, v3. Die erste Spalte der Determinante V kann infolge von 
Gleichung (8) auch mit den b;;, statt mit den a;;, gebildet werden. 


b) Wenn die Ausgangsfläche z= f(z, y) und die auf ihr vorgegebenen Werte u, v der Rich- 
. tungsbedingung nicht genügen (R-+0), lassen sich die gesuchten Ableitungen 4,, v; aus den 
Gleichungen (7) eindeutig berechnen. Die Ausgangspunkte w, v bleiben in diesem Falle bis auf 
die Ungleichung R=+ 0 willkürlich und sind keiner Verträglichkeitsbedingung unterworfen. 

Wir bezeichnen eine Fläche bezüglich der auf ihr vorgegebenen Werte «, v als charakte- 
ristisch, wenn die Bedingungen R=V=-0 erfüllt sind, und für R=F0 als nicht charakte- 
ristisch. Die Richtungsbedingung (8) ordnet jedem Punkt P(z,y,2) mit vorgegebenen Werten 
u, v einen Kegel 2. Klasse (charakteristischer Kegel) mit den Normalenvektoren (9, 9, —1) 
zu. Wir bezeichnen die Differentialgleichungen (6) für das gegebene Wertesystem z, %, 2, u, v 
als hyperbolisch, wenn der charakteristische Kegel reell ist und nicht entartet. Alle nach- 
folgenden Betrachtungen beziehen sich auf hyperbolische Bereiche. Außerdem bezeichnen wir 
ein Linienelement als raumartig, wenn es außerhalb des charakteristischen Kegels liegt, wenn 
also 2 Tangentenebenen des charakteristischen Kegels (charakteristische Flächenele- 
mente) durch das Linienelement gelegt werden können (Bild 5). 

Die charakteristischen Flächen sind Hüllflächen charakteristischer Kegel, d.h. sie werden 
in jedem Punkt von einem charakteristischen Kegel berührt. Wenn die Koeffizienten a;, biz 
Funktionen von z, y, z und auch von «, v sind (quasilineare Differentialgleichungen (6)), 
ändern sich die charakteristischen Kegel und Flächen bei Änderung der vorgegebenen Funktions- 
werte u, v. Wenn dagegen die Koeffizienten a,;, dj; Funktionen nur von &, %, z, nicht aber von u, v 
sind (lineare Differentialgleichungen (6)), liegen die charakteristischen Kegel und Flächen von 
vornherein, d.h. unabhängig von den vorgegebenen Funktionswerten u, v fest. 

Die Charakteristikentheorie [17] [18] führt zu folgenden allgemeinen Sätzen für die in 
einem Raumbereich (R) existierenden Lösungen «, v der Differentialgleichungen (6): 


1. Bei vorgegebener Lösung u, v gehen durch jede Kurve mit raumartigen Linienelementen 
2 charakteristische Flächen. Wenn sich die Kurve auf einen Punkt P? zusammenzieht, entsteht 
eine charakteristische Fläche, welche in P den charakteristischen Kegel längs aller Erzeugenden 
berührt und wegen ihrer kegelartigen Gestalt als charakteristisches Konoid Bachnes 
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2. Eine Fläche (F) mit der Randkurve r sei bezüglich vorgegebener Funktionswerte u, ® 
auf (F) nicht charakteristisch; alle Linienelemente von (F) seien raumartig. Wenn dann eine 
Lösung der Differentialgleichungen (6) mit den auf (F) vorgegebenen Anfangswerten existiert, 
ist sie durch diese in einem Bestimmtheitsbereich (B) eindeutig festgelegt. Der Bereich (B) 
wird auf jeder der beiden Seiten von (7) von der ‚‚inneren‘ der 2 durch die Randkurve r gehenden 
charakteristischen Flächen begrenzt (— in Bild 7 ausgezogen —). 


N - ” 7. Ausgangsfläche F mit Bestimmtheits- 
m ” Bild 5. Raumartiges Linienelement 7 Bild 6, Charakteristisches Konoid/ und Einflußbereich 


3. Wir betrachten verschiedene Lösungen und ihre Anfangswerte auf einer Fläche (F*), 
welche die Fläche (F) als Teil enthält. Wenn die Anfangswerte auf (F*) nur außerhalb von (F) 
geändert werden, ändert sich die Lösung nur außerhalb des zu (F) gehörigen Bestimmtheits- 
bereichs. Wenn die Anfangswerte auf (F*) nur innerhalb von (F)) geändert werden, ändert sich 
die Lösung nur innerhalb des Einflußbereichs, der von den „äußeren“ durch r gehenden 
charakteristischen Flächen (— in Bild 7 strichpunktiert —) begrenzt wird. Wenn sich (F) auf 
einen Punkt P zusammenzieht, entartet der Bestimmtheitsbereich in den Punkt P und der 
Einflußbereich in das von P ausgehende charakteristische Konoid. 


3.2. Gitterkonstruktion bei Systemen von 2 quasilinearen Differentialgleiehungen erster 
Ordnung if 
Nach Satz 2 am Schluß von Ziff. 3,1 ist eine Lösung durch ihre Anfangswerte u, v auf einer 
Ausgangsfläche (F) mit lauter raumartigen Linienelementen in dem Bestimmtheitsbereich (B) 
eindeutig festgelegt. Wir wollen diese Lösung jetzt näherungsweise durch ein Differenzen- 
. verfahren mittels einer Gitterkonstruktion bestimmen: 
r “ Das Gitter besteht aus einer Folge von Ebenen 


YO, y=2e, y=4e, vn. y=2ne 


und aus einer Folge von Ebenen &= const sowie aus den charakteristischen Flächen der einen 
Schar, welche von den Schnittkurven l(y = const) der Ausgangsfläche ausgeben (Bild 8), 

Die Gitterpunkte P;, auf der Ausgangsfläche und ihre Funktionswerte (Pa), v(Pa) 
sind vorgegeben. 


In jedem Punkt P;, ergeben sich aus der Richtungsgleichung (8) 2 Richtungen dir, PA 


nach Einsetzen der durch die Schnittkurven vorgegebenen Richtungen p = --. An Stelle IE SS 
strengen Werte p (Tangenten) benützen wir dabei die Näherungswerte (Sehnen) 
Dr 2 (P,)—2 (Ps) 

u (Pa)—® (Pı2) i i 

Die Schnittpunkte P;, der durch die Richtungen gie und g4,a+1 Bestim timmten Geraden m 


mi,» +1 nehmen wir als Näherungen für Punkte der charakteris her hen Flächen des. ‚Gitter 


Näherungen der Funktionswerte v(P;.), v(B,) ergeben ich urch Anwendung der 
träglichkeitsbedingung WA a ine Linienelementeig und ı an Er wieder 
vhs 


usf. 


FE J ech 5 \ 4 
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h die N eriagswertei Pir und dir bzw. gi, und die Differentialquotienten ar usf, durch 
er enzenguotienten, z.B: ; Ei 


3) u (P..) —ıu (B, ER x=const 
ylmız Y (2) | 


R Be Er _ulP.)— UP) 
E;. Oylmis y(P Pe) upr) 


Br, fen) u(Pa%)— u(P}5) 
Beer) !., Fr (Pa) —% (Pı3) $ 


Dadurch erhält man für jeden 


Punkt P; je 2 lineare Bestim- 
mungsgleichungen zur Berechnung 
der gesuchten Näherungswerte 


Die Punkte P;; liegen ebenso 
wie die Ausgangspunkte Py,lür 
‘—=const in parallelen |Ebenen 
%= const, jedoch für k=const- 
nicht in parallelen Ebeneny=const. 
Wir gehen daher zu den Schnitt- Bild 8. Schema für die hir 


punkten P#, der Geraden P;; Pj, mit den Ebenen 


‘ 
‘ 


IE HB, Yesdle; y=edie, ey =lznl)e 


über und ordnen ihnen die durch lineare Interpolation sich ergebenden Funktionswerte u (P}.). 
v(Pi5) zu. Mit den Punkten Pj; und ihren Funktionswerten läßt sich dann die Gitterkonstruk- 
tion in derselben Weise wie bei den Ausgangspunkten P;; durchführen. 


3.3 Gitterkonstruktion bei einer quasilinearen Differentialgleiehung zweiter Ordnung 
Bei Systemen von 3 u: Differentialgleichungen 


en - r | Zaun + ben + Lou en 
ER e E Be ur + 2 0 u 2 69: ur —d4,> 110) 


BR 4 NR DER Say mt 3 bar 0 + I Cyr Wr — 9a 
für die 3 gesuchten Funktionen u (2, 9,2), v (2,2), w (x, y,2) treten an Stelle der Gleichungen (8) 
5 und (9) die ‚Riehtungsbedingung 


FE | ap + ag — Us bup +beg—bıs CuPt 21 — Cıs 
Ay Pt Ag, 9 — Gas b21,P + bee a — bs Pt 4 — Ca 
az P + Ag — Aya bp + bag — bs; 31 P + 632 1 — O3 
Er una die Verträglichkeitsbedingung A SR 


ri ou a 
ana ri 


UN ZSE = LATE ROEN, 


InaKen baten 


öw 


+ ehe 125, | 


un 


du 


ee Fessr En ed 


Öw 


U 92 


f 


Bu er "si, 


"iR LE 
. ng) 


SERDNT IN 
354 Sauer, Drei OEBERBHT? Pb lemd der Or akteristikontheone Ba 80 ee, ine ta, Sec, ku, 


Wir beschränken uns auf den praktisch wichtigsten Spezialfall, namlich die ee 
Differentialgleichung zweiter Ordnung 


Be Ark 1 G | | (Ax — Ari) RR ITE. > Orden (13) 


\ für die gesuchte Funktion @(z,%,2). Die Koeffizienten A4;,; und @ sind Funktionen von &, y,z 
und 91,3, 9; (— nicht aber von @ selbst —). Durch den Ansatz 


Up UPS |. ;.. %: Ve . . (14) 


wird die Differentialgleichung & äquivalent mit dem System der 3 Differentialgleichungen 
Ayıtı + Arte + Ass ts+ Aytı + Ag + Agsva + Au wı + Aa w + Ay w—=@, 


Ug — u] ==, 
© 0g — ==) 
R Hierzu folgt aus Gleichung (11) die Richtungsbedingung 
| R= AıuP+242Pg+ Aa —24A3P—2Ayg+ As — 0 | ..2.22022.(15) 
und aus Gleichung (12) die Verträglichkeitsbedirgung 
du du | 66V ÖW dw 
= Aus, 3 Aulse tg) 4 Ange + An An Are) 


Pr 
2 
re 
] 
h, 
ha 
1 


öw (16). 
(2 A ae al) —6=0 


Sie läßt sich mit Hilfe der aus 


Öu 6% öw 
day + U3g > 35 Aut WPp: zz utrWp 


su m ÖV öw 6% 
a 


N + (2A, —2 AP — An 0 


. Die EEE RRENAE,.. (15) ist wieder in 9, q quadratisch, die ee eh also 
"wie in Ziff.3,1 Kegel 2. Klasse. Die Zusatzbedingung (17) bringt zum Ausdruck, daß die Projescat h 
TAB des Vektors rot (u, v, w) auf die Flächennormale (p, ,, — 1) verschwindet. 

E a Die Gitterkonstruktion von Ziff. 3,2 ändert sich nur insofern, als jetzt in jedem Punkt. Pr 
Be. die 3 Unbekannten u(P,), v(B;), 29.) zu bestimmen sind. Da aber die Verträglichkeits- j 


bedingung (18) nur die beiden Ableitungen = r je nicht aber, enthält, liefert 8 jeweils 2 
lineare Bestimmungsgleichungen für die 2 Unbekannten v( Px), w(P;,). Nach deren erechnung 
erhält man aus der als Differenzengleichung umgeschriebenen Zusatzbedingung (17) eine lineare 
Bestimmungsgleichung für die dritte Unbekannte u(P). Se 
3.4. Erläuterung an der stationären Übersehallströmung RO 

Wirerläutern die 3-dimensionale Charakteristikentheorie am Beispiel der stationären räum- DE 
lichen Überschallströmung. Dabei beschränken wir uns der Kürze halber auf den Sonderfall de La 
wirbelfreien und daher isentropischen Strömung. Es existiert dann ein ( 
(2, y,z), aus dem sich die Geschwindigkeitskomponenten % u,v,W. 
” K Die erste der Gleichungen (1) liefert für dar Potenti lin cartesisch che 
4% Differentialgleichung Bar 


ER a (u — a?) P1 + (w — a?) Pa2 + (w? ea as) Pa | 


1  n 
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oder in Zylinderkoordinaten (Bild 1) die nichthomogene Ditferentialgleichung 

2 uU 2 uv Ve N 
Ra) gut (— a?) 994 (W— a?) p35+29Ww@pa5+ r' Pa zZ 7 Pa (w*+a?) (20); 


in Gleichung (20) sind die Indizes 1, 2, 3 den Koordinaten o, r, z zugeordnet und an Stelle von 
Gleichung (14) tritt 1 


ir nis 0 ea 


Die Gitterkonstruktion von Ziff.3,31äßt sich unmittelbar auf Gleichung (19) oder (20) anwen- 
den; im letzteren Fall wird das Gitter außer von charakteristischen Flächen von Ebenen ® —const 
durch die Achse desZylinderkoordinatensystems und von koaxialen Zylindern r= const aufgebaut. 
Die charakteristischen Kegel, Flächen und Konoide sind die Mach-Kegel, Mach-Flächen 
und Mach-Konoide der Überschallströmung. Sie sind dadurch gekennzeichnet, daß ihre Flächen- 
elemente den Mach-Winkel & (sin &= Schallgeschwindigkeit a: Strömungsgeschwindigkeit 
Yur + v2 + w2) mit der Strömungsrichtung einschließen oder, was dasselbe bedeutet, daß die 
zu den Flächenelementen senkrechte Komponente des Geschwindigkeitsvektors gleich der ört- 
lichen Schallgeschwindigkeit ist. Hieraus ergibt sich anschaulich-physikalisch in Übereinstim- 
mung mit. Satz’3 am Ende von Ziff. 3,1, daß eine in einem Punkt P wirksame kleine Störung sich 
nur im Innern des von P stromabwärts laufenden Mach-Konoids (Einflußbereich) ausbreitet. 

Die Gitterkonstruktion läßt sich sinngemäß auf Anfangswertprobleme übertragen, bei 
denen, wie dies bei Aufgaben der Gasdynamik oft der Fallist, Anfangswerte auf einer charakte- 
ristischen Fläche gegeben sind und außerdem eine Randbedingung (z. B. tangentiale Strömung) 
an einer anderen Fläche zu erfüllen ist. Ferner läßt sich die entsprechend verallgemeinerte Gitter- 
konstruktion auch auf die die 4 Differentialgleichungen (1) der nichtisentropischen Überschall- 
strömung mit den 4 gesuchten Funktionen u, v, w,s anwenden. Als charakteristische Flächen 
hat man hier einerseits die Mach-Flächen wie bei der isentropischen Strömung und andererseits 
die von Stromlinien erzeugten Flächen. 

Was wird aus der Gitterkonstruktion in den in Ziff. 2 behandelten Spezialfällen mit be- 
sonderen Symmetrieeigenschaften ? Bei drehsymmetrischen Strömungen braucht man die Gitter- 
konstruktion nur in einer Meridianebene (o = 0) durchzuführen. Sie geht dabei in die bekannte 
2-dimensionale Charakteristikenkonstruktion über; die drehsymmetrischen Mach-Flächen schnei- 
den die Meridianebene nach den Charakteristiken des 2-dimensionalen Problems mit den 2 unab- 
hängigen Veränderlichen r, 2. Bei Kegelsymmetrie wird der Arbeitsaufwand der Gitterkonstruk- 
tion sehr verringert, da die für einen Gitterpunkt Pf; gefundenen Werte zugleich für alle Punkte 
der Halbgeraden 0 Pf}, (0 = Symmetriezentrum) 8 Bölteh und daher nur eine charakteristische 
Fläche konstruiert werden muß. Man kommt jedoch in Polarkoordinaten R, d,© (Bild 2) nicht 
zu einem 2-dimensionalen Charakteristikenverfahren mit den 2 unabhängigen Veränderlichen 
®,@. Eine in einem Punkt P angreifende Störung geht nämlich’ wegen der Kegelsymmetrie 
sogleich auf dieganze Halbgerade O P über und breitet sich daher auf das Innere des Mach-Kegels 
mit der Spitze O aus. Alle Einflußgebiete fallen also mit diesem Kegel zusammen und man hat 
daher auf einer Schnittkugel R = const mit den 2 unabhängigen Veränderlichen 9, © im Innern 
des Mach-Kegels keine reellen Charakteristiken. 

In demin der Einleitung genannten Buch von A. Ferri [8] wird die Richtungsbedingung (15) 
und die Verträglichkeitsbedingung (18) für die Potentialgleichung (20) der wirbelfreien Über- 
schallströmung ohne Bezugnahme auf die allgemeine Charakteristikentheorie hergeleitet. Mit 
Hilfe einer Gitterkonstruktion wird dann die Strömung um einen aus Drehkegel und Drehzylinder 
zusammengesetzten Drehkörper unter kleinem Anstellwinkel näherungsweise berechnet. Das Er- 
gebnis bleibt jedoch aus 2 Gründen fragwürdig: Erstens wird die Strömung um den Kegel mit 
einem halblinearen Verfahren (vgl. Ziff. 2,2) ermittelt und die strengere dreidimensionale Gitter- 
'konstruktion erst bei der Expansionsströmung am Zylinder verwendet. Zweitens müßte die durch 
den Anstellwinkel hervorgerufene Wirbelbildung berücksichtigt und demgemäß das allgemeinere 
Differentialgleichungssystem (1) zugrunde gelegt werden. 


4. Zusammenfassung 

Im ersten Teil der Arbeit wurden 3-dimensionale Probleme der Charakteristikentheorie be- 
trachtet, welche sich wegen gewisser Symmetrien auf 2-dimensionale Probleme reduzieren lassen, 
und Nachbarlösungen solcher symmetrischen Lösungen (nahezu drehsymmetrische bzw. nahezu 
kegelsymmetrische Lösungen) mit einem halblinearen Näherungsverfahren ebenfalls auf 2-dimen- 
sionale Aufgaben zurückgeführt. . 

Im zweiten Teil wurde für die im eigentlichen Sinne 3-dimensionalen Aufgaben mit Hilfe 
‚einer Gitterkonstruktion ein Differenzenverfahren entwickelt, das eine naturgemäße Verall- 
Kr meipernug des bekannten 2-dimensionalen Massauschen Charakterjstikenverfahrens darstellt. 
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Zunächst wurden die für das Verständnis des Verfahrens notwendigen Grundbegriffe der Cha- 
rakteristikentheorie zusammengestellt und hierauf die einzelnen Schritte des Verfahrens im Sinne 
einer Gebrauchsanweisung ausführlich erörtert. 
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Eine verschärfte Fehlerabschäbuns 


zum Extrapolationsverfahren von Adams) 
Von Johannes Weißinger in Hamburg 


Für das A da m’sche Verfahren zur Integration gewöhnlicher Differentialgleichungen r-ter Ordnung sind 
zwei wesentlich verschiedene Fehlerabschätzungen bekannt ([1], [2)). Da die letztere zwar schärfer ist, aber- 
ee ie Lücke enthält, wird nach dem Verfahren von [1] eine ähnlich scharfe Schranke wie in [2] 

ergeleitet 

For the method of Adams performing the integration of ordinary differential equations of the r-th 

. order, two estimations of the error are known that differ essentially from one another (11), [2]). The latier 
e * being more effieient, but marked by a certain inaccuracy, limits similarly efiaen as Bla are deduced 
following the line of [1]. 

Pour la methode de Ada m s d’integration d’Equations differentielles ordinaires un ordre r on connadt 
deux proc&des ü l’estimation d’erreur ([1], [2]), essentiellement differents l’un de l’autre. [2], bien qu Eee 
soil plus efficace, manque de clart& dans un point principal. A l’aide du procede 1 on BoD: une gelpmantgn Meer: 
d’une efficacilE semblable & [2]. BE 

ra pemerna mpocTsx AuhepeunnansHex ypaBHeHnä nopsanka y no MeTony IE, 
UMEIOTCH ABA CYIINeCTBEHHO PA3AHYHLIe CIIOCOo0a eMeHKm morpemmocrtu ([1], [2)).. Tax KARO 
BTOpoü cIOoCo6 bonee TOYeH, HO CENEP:KUT ONHY UPHHHHIMANBHYIO HESICHOCTB, TO 3NeCh oe 

- IPoBoAuTeH 10 c6ocody [1] onenka NOrPe uHocTH TAKOH-IKEe TOIHOCTE, KaR “no emocody 2) 


fi 
A 


1, Bezeichnungen und Voraussetzun; en HR a NE a 
y(2) = exakte Lösung der Diiferentialgleichung y’= f(x, y); f(x, y) habe stetige partiel 

Ableitungen bis zur Ordnung r +1. Anfangsbedingung: y(z)=Y- 
K = Lipschitzkonstante von 2 a f(&, y) —f(®, yı)\< K 97 


F, +1 = Konstante, so daß F, Eyilz, say)l- 


y(2) = Näherungslösung nach dem n Bitrapeıeiion SY verfahren r-tei j 
% =-ıH+tih; h= Schrittweite. A 


!) Teileines auf der Jahrestagung der GaMM am 19. 
3. Abschnitt gemachten Bali ar en N er dgganf ge 


Ferner werde gesetzt r | 
2: ya = Yi» 
n Ya) = Dr Yir Üi, Y:) fh: fe: Yı) —f. 


Nach dem Extrapolationsverfahren ergeben sich die Y aus r--1 Anfangswerten %, Yıs - :: Yr 
mittels ) 


x i WERE NDR 


e=0 


B Dabei ist a $ > 
# % ” r ver | n 
Be PEN... () RR [ser - (u+o—1)du. Da 
Br‘. 2 Ri 
. Schließlich sei & ‚eine Schranke für den Rn aeNer‘ Ken“. " 
tn : e>|%—yl;, Et PRRSERRR > 5 Ni 
2 2 N b v 1 - i ae 
en. 8%. Die v. Misessche Fehlerabschätzung N 
E Um den Fehler “ra 
Br: 3 er r 4 “ +1 t = 
= | +1 = Nrı „Yirı = Yiyrı Ti ih: en ))dE 
# abzuschätzen, schreibt v. Mises 1] 
“rat oel ee ONE Ne OT); 
E2 „ wobei der ‚‚Quadraturfehler“ Er 

Er EHE « 
5 25% RER Ze ehren = f(x y(&)) 


‚durch 
lqil s gl +30, 11 By 


 abgeschätzt werden ‚kann. Aus der Lipschitzbedingung ergibt sich dann die Differenzen- 


eehung y 
le; wu zZ l&l E= hK 2 [&e| A + ES OR AS (2,2); 
‚ und jede Lösung E; der entsprechenden Gleichung 
N VAR = Eirı=EcthK I ledEi +4. ER Re 
Nefert eine Majorante Für die lei und damit eine Fohlrabschtzung, falls nur RE 
BET R | Er =0. na. a (2,4) van 
.. Dre zugehörige ‚charakteristische „eichung - 
5 nr ER De et #+Hi a —ıK 2, Bu... . (25) 
Re fe hen R sg en, 4 | za 
hat nach ser Cartesischen Zeichenregel genau eine BEE Wurzel 29, gie offenbar überdies 5. 
en als. ir ist. Dann ist | 2: de. 
I Eu ee | sr 10,6) 
e RESET. , EB, Oi ug w- 3 Rx P3 BET RE (2% ) | 


Pr 


einer a tenigen eine partikuläre Lösung von © 3). Wählt man 


Beier. ....(2,7) 


° Y 2 v x Pe d \ « a’: 
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An. dieser Abschatzune ist unbefriedigend, daß bei Festhaltung der von r habhängiäen‘ 
Größen die Fehlerschranke mit wachsendem r zunimmt, und zwar sehr kräftig. Aus der sofort 
ersichtlichen Ungleichung ; Rt. 


cs I o.= 1, 2,.\.. Wan (2,9) 


IN 


20 

ergibt sich nämlich 

9r 5 

= <1+2% — san ncıt: Le RE SE LU NE 
en, o= x 2 o=0 | 

Daraus folgt nicht nur, daß der vom Re... sherrührende Term in (2,8) mil r —oo gegen © 

. geht, sondern auch der vom Quadraturfehler herrührende, wenn man annimmt, daß F,; ı für 

r— © nicht oder jedenfalls nicht zu schnell gegen Null geht. Denn dieser Term wächst stärker 


als 
h\r u N 
. Se ie ) | 


was bei genügend großem « mit r gegen oo geht. 
Nun lassen sich zwar leicht Beispiele angeben, wo bei gleichem & das Verfahren höherer 
Ordnung wirklich eine schlechtere Näherungslösung liefert. Aber daß eine Vergrößerung von r 
sich so ungünstig auswirkt — insbesondere im Term des Quadraturfehlers — wie in der obigen 

3 Fehlerschranke, erscheint doch wenig plausibel. 
Bevor wir die Abschätzung zu verbessern suchen, noch zwei später Depkanche Bemerkungen ! 
Erstens sei darauf hingewiesen, daß man die Ungleichung (2,8) ganz trivial ohne Betrach- 
tung des charakteristischen Polynoms (2,5) erhält, wenn man die Differenzengleichung (2, 3) 
ersetzt durch die Gleichung 


= ”. 
By ı=(I+hRA)BI+g, n 
die offenbar erst recht eine Majorante für die |e;| liefert und gerade die rechte Seite von (2,8) - 
zur Lösung hat. Man kann (2,8) auf, diese Weise sogar,noch etwas verschärfen. Sicher gilt 


nämlich 
ler+ıl<(I+RKA)o+g 
f und durch Induktion folgt 


lei = = |Y Ye Yıl Se (1 +4 KAT a [(1 #RKA)NT-' u 1], (= r, TH ik ek (2,11). 
Zweitens gilt bekanntlich für beliebiges ce und Du | \ x ’ 


a 


und für genügend kleine c unterscheiden sich beide Seiten beliebig wenig. Wegen 32 == 


ar ne 


LT fo 2 
kann man daher statt (2,8) etwas weniger scharf auch schreiben 


| 4 ger u 2) [e ‚AR (e; »— %,) Et 


+kal N 
wobei sich diese Schranke von der alten für genügend kleines h beliebig wenig unterscheidet. 


SR u 


3. Die Abschätzung von Tollmien undFike Pi N 
Die zuletzt erhaltene Abschätzung erinnert stark an folgenden bekannten Satz : aus der 2 
Theorie der Differentialgleichungen: Genügen y(x) bzw. z(x) den Ditferentialgleichungen BIER 
y=/(n9); bzwizmo, 2) BR 2 SUR ER te 
mit den Anfangsbedingungen BE AR N 
Y(%) = Yo» 2(&)= 20,” 
und ist X die Lipschitzkonstante von f(x, y) sowie | 


M= Max [Fa =12)) Caleta), 
so gilt 


I a ee (il I % 


Bekhritiel 
sen Pie 
MR AN, 


DEs liegt nahe — und dieser Weg ist in [2] und [3] b 
mit En obige Biähorinaeit Pu EE 72) und ue 


« 
Fi 
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durch die Werte Fh-r=f&ü- —ı) (k=0,1,..., r) gelegten Interpolationspolynom 
P;(x) der Ordnung r zu identifizieren und dann auf Grund der bekannten Fehlerformel 


F N. a Ne ae) C 
’ Tan er Par % 2 VRR ee 
& I Y(2)) P;(x) (r +1)! rel %(8)) Pr (3,1) 
\ für den Interpolationsfehler zu setzen 
. n dr +1 1 23 ei m 
M = h"t1! Max gR sl (2, Y(a))| HER. L:. 
1 Die so erhaltene Formel 
%—yı < Karin .) AR I" al Pr +1 (ea u) ı) N RR (3,2) 
RK 
bzw. die durch eine Zusatzbetrachtung zu erhaltende verschärfte Formel 
| H=ylss "nt Kr ET ie BER se u (30 
wird mit (2,12) im wesentlichen identisch, wenn man dort F,.ı, durch F, +ı und A durch I 
ersetzt. Die Größe A, die für große r die Abschätzung so verschlechtert, tritt jetzt also gar 
nicht auf. 
Leider ist der eben skizzierte Weg nicht gangbar. Denn die Anwendung der Formel (3,1) 
setzt (r+1)-malige Differenzierbarkeit von Y(«) im Intervall (&—r, 2 ı) voraus. Andrerseits 
ist Y(x) = P;(e) in (&, 8% +1) %(@)= Pi ı(la)in (5 ı, %;) usw.;und da P;(e), P; —1(®) im 
Punkte «; im allgemeinen nicht übereinstimmen, so haben die erste und alle höheren Ableitun- 
> 2% Ir +1 
gen von Y(x) Sprünge in den &;. Die in (3,3) auftretende Größe F, + ı = Max ee %(ex)) 
ist also gar nicht definiert. In [2] wird darauf hingewiesen, daß man im Falle des Nicht- 
existierens der Ableitungen die Differenzenquotienten nehmen muß. Eine Abschätzung für 
diese liegt aber bisher nicht vor. 
Andrerseits wäre natürlich eine ähnliche Formel wie (3,3), die A nicht enthält, sehr er- 
wünscht, da sie — insbesondere für große Integrationsintervalle x; — x, — bedeutend schärfere 
Schranken als.die v. Misessche Formel liefern würde. 
Ar 4. Versehärfung der v. Misesschen Abschätzung 
4 Man kann nicht erwarten, daß die Formel (2,12) (und daher erst recht nicht die schärfere 

Formel (2,8)) richtig bleibt, wenn man einfach A = 1 setzt. Für das Problem 

y=-y 0-0, also ya)=0, q9=0, K=l 
würde man dann nämlich erhalten 
I Seth 
während man aus den Anfangswerten 
Male, Yr-i- Sl, Y=l-1N : 
wegen &, = (— 1)?|a,| errechnen würde 
Y 
IrHı=s+h2 |ale=(1+ Ahle, 
e=0 
so daß also . 
1+Ah<sertuh 

gelten müßte, was z.B. für A>r +1 und genügend kleine % sicher nicht richtig ist. Statt 

dessen beweisen wir den 
| Satz: Zu beliebig kleinem ö, > 0, 6, > 0 gibt es (bei beliebigem festem r) stetseein H >0, 
 . sodaßfür alle A mit O<h<H und für beliebig große « gilt 
Rt, h Y 
EEE EN 1l........ 


wobei 
Be. eF+—=sll+ö,), Ar=A(N-+ÖS,), 

Be: mit anderen Worten: In (2,12) darf man zwar A nicht genau durch 1, sondern nur (für entspre- 

chend kleine Schrittweite %) durch eine beliebig nahe bei 1 gelegene Zahl A* ersetzen, wenn man 

ar gleichzeitig BU; wenig vergrößert. a: A; 


U TA % Lu le ie 
d P Zur - ‘ 


“ e 


Weißinger, mon abschätzung z 


Zum Beweise summieren wir (2, 1) über die Eizes i bis Fe Ds er] alte 


% 5 I+r ; A 

fr 1 2 > be [Fi ze e—fi—o] +20. a ERRETT Zr ARE 

N N ! = 
Ri, $ mil | Bi © 
Be) | Kto 4:4 Go für 0<o<r—1 | 1 
Bar be WHO, tl . r<o<t. ae 


ee l<o zu \ 


Dabei soll > r eine beliebige, aber feste natürliche Zahl sein. 


- e ; ” ” 

Unter Benutzung der Lipecht Be und der früheren Bezeichnungen ergibt sich daraus en 
” als . 
lö+ıl <a +hE 2 ld.) lic] + (I A 1) 2 We Se A 2 


so daß man eine Majörante der |e;| in jeder Losung E,; der ERS. AN a 


Eiy=BıthE 3 BB ne. 0 

‚erhält, welche den Anfangsbedingungen FR 2 
nl. Drol...ır- Ve 

Bi genügt. Wir werden schärfer mit festem «* sogar 2 Br 
Di .BK2e*>|e|, De A 0Eı B S  e Nas .: (4,6) > 


fordern. Die Charakters sche Gleichung 
zit Hl — gr — ur "4 a Re 


hat wieder genau eine Wurzel 2, rechts von 1, so daß 


mit beliebigem Ü’' und 


vi i \ 

' a: LTR. er ee \ f 2 
i re et a 
auch die Anfangsbedingungen (4,6) befriedigt. Wir bekommen daher folgende Abschätzung 


, ) 4 — yı < € lee 1342 A LE ISTZ 0). 
u zu Aussagen über 27 zu kommen, schreibt man (4, 7) am besten) in der Form: vn ur BR 
Ye N : N '2r(zi+1 1) ö IR ; ver Eae 
Br» { ae : ON Re (4,10) ; 
PA aD SER Warsce Sr HR, 2: ® Er 


c=0 UP RER, I 
und untersucht die Funktion BER Setzt man noch j N 
| BE an head ln tr \ a ls 
50 daß also i Ri 
u: Hast, b55+b,=1} und damit Isel ns bel > 
0 gilt, so wird 


% 


B-i410 44 5 wann. | 


a 


Penner man unmittelbar, daß WIE bei festem 2>1 mit wachsendem 1 monoton zunimmt, so 
daß 27 bei a Kh mit wachsendem ! abnimmt, die Abschätzung (4,9) also verbessert wird 
(falls nicht &* zu groß wird). 

Als Ableitung von x(z) an der Stelle 1 findet man aus (4,10) sofort 


Eee 


| | B Ar’ 

so daß für kleine K’h näherungsweise gilt 

Bi 41-4 A*Kh. 

FE Für die zweite Ableitung ergibt sich - 

z B? el et Ele 

# WI {uH El dd} +2 S oflölt la 1) +2 rl. 
0 daß älso wegen (4,12) : 
B rc K1>0 
a gilt. Für Eulen Klemme Kh liegt die Kurve x (2) daher oberhalb ihrer Tangente in z=1, d.h. 

b es ist 

| BEER en ln (4,15) 
und folglich 
j i <{i N 
z i 


' woraus zusammen mit (4,9) die behauptete Formel (4,1) folgt. Es bleiben nur noch die Behaup- 
tungen über A* und s* zu beweisen. 
Zunächst kann man (bei festem, Er ee r) 


Ar + 7 Ziel + ID | 
offenbar durch Wahl eines genügend großen / beliebig nahe bei 1 N Bei nunmehr festem / \f | 
| beschränken wir uns auf so kleine h, daß (4,15) gilt, und bestimmten &* nach (2,11) zu Er 
E j x | R' 
$ =s(1+hKA) ht? Ben KBISRKAME IHN: . (4,16). R 
| 1 i Da L jetzt fest ist, kann s* durch Beschränkung auf genügend kleine h beliebig nahe bei & REsLEn ” 
> werden. Damit ist der Satz vollständig bewiesen. 
015 i" 
ri) 
00 
| a0s|— e 
E 0 105 110 115 z 120 ; 
1 1. Die Funkin (2) für verschiedene r und 2 RR Bild 2. Die nr er ua für r=4 Pe } ; 
$ Wr-% FR ; 
ücklich sei Bi: daß der x- bzw. z-Bereich, in dem (4,15) gilt, von l abhängt Bu 
ür 1— oo erhält man nämlich IR - 


ut ‚Null a Ri 
nr: 4 e BET Da 
| a: Mi 

a | eine 
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mit den Ableitungen 
2 r—1 > 
„=1, = zu En 
o=0 z 


x ee... / 

folgt. In der Praxis wird man aber ohnehin I nicht zu groß wählen, um &* nicht zu groß werden 
Br zu lassen, und da auch KA selten größer als 0,15 sein wird, wird die Ungleebunz (4,15) praktisch 
; meistens erfüllt sein, wie ein Blick auf O2 Bilder 1 und 2 lehrt. 


| 5. Zahlen werte }" 

E Beispiel: Wir stellen zunächst noch einmal die praktisch wichtigsten Formeln des vorigen 
Y Abschnitts zusammen und geben die Zahlenwerte der darin vorkommenden Konstanten an. 

Mit den in 1. eingeführten Bezeichnungen läßt sich für das Extrapolationsverfahren r-ter 


wi’ Ordnung or Fehler |; — %;| nach dem ;-ten Schritt abschätzen durch 
| HH —y|<s e* rein an 6,1) 
Bi; mit 
Ei VE FE RR NET RR LS Fa EA TE SRH RETET Ike 9 (5,2) 
A 4 Wet Sy 
& (+ RAY + aaa) I] ESS RS N He 


Hierin ist Zeine beliebige natürliche Zahl r< 1 < i und 2, die positive Wurzel der Gleichung (4,10), 
die für einige Fälle aus den Bildern 1 und 2 entnommen werden kann. In der Praxis wird meist 
2 <1-+hKA* sein, so daß man dann etwas weniger scharf 


IH yıl< ori -ARANH mr Hr RRAN 1] BE 


"hKA* N 
Er. f schreiben kann. Die Konstanten a, ı, 4, 4*, können der folgenden Tabelle entnommen 
Be werden: 
| dr | 4 | a+nar-n 
0 1.0000 1.0000 0 Pr: 
l. 0.5000 2.0000 1.0000 
2 0.4167 3.6667 1.8333 F 
3 0.3750 6.6667 3.6667 
4 0.3486 12.2444 8.5472 
5 0.3299 22.8222 17.8000 er N 
6 0.3169 N 
. Zum Schluß wollen wir ein Beispiel betrachten, das auch in [1] und [3] behandelt wordenist: 
Ye > 0)=1,- r— AIR 0.02, 0<x<200 
Be nr = ee 
e F,—= 8840, Kh=0.006, v.g = a F5 — 8654107 Fa 
ve.0.9:105% 


Nach v.Mises ergibt sich für den Fehler nach 100 Schritten die Schranke Si =19907 10-, = 
während sich nach (2,3), falls man unzulässigerweise F,;, statt F, +1 benutzt, die rund A 
80mal bessere Schranke | 90) < 25 : 10% ergeben würde (vel. [3))- Nach unseren obigen Formeln Er 
Ice erhält man für > =7 mit . = eh 


010 10 Fe ae 


Beer und a 
A* — 2,0684, A*Kh= 0,01241, 2 1,01241 12 ö 
die Schranke | Re 


h * 
# . r f 
‚ .% l h 2 
pr „ 
ar: > E 
AR \ ’ 
rn. 
u 4 > 1 
2 4 
*, hi ' } 
. 
= ” > 


1,865 - 107 
< 2,454 - 10-8 - (1,01241)700 | SET 

| iool < 2,45 0 ( 0 220 0,0124 
— 8,44 :10-% + 36,62 10-°— 45,1. 10, 
die nicht wesentlich schlechter als 25 - 10-8 ist. 
Gelegentlich kann man die Abschätzung noch v erbesser 
Berechnung von e* nicht die für das ganze Anteg 
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für 2, S@z <a,+, gültigen Konstanten X, F,.., zu nehmen braucht. In unserm Beispiel ist 
an 
dy (y+»)? 1,40? 
so daß sich mit AK = 0,0045, AhK = 0,0551 
sr 1,886 --10-6, lan Asshrid*t, 
also nur eine unwesentliche Verbesserung ergibt. 
Abschließend sei noch darauf hingewiesen, daß der Grundgedanke der geschilderten Feh- 


lerabschätzung offenbar allgemeiner anwendbar ist, wofür in einer weiteren Mitteilung einige 
Beispiele gegeben werden sollen. 


0,225 für 0% <S022% 
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Zur Berechnung der Torsionseigenschwingungen von 
Maschinenwellen 
Von Alfred Vogel in Stuttgart 


Es wird die dynamische Matrix aufgestellt, deren charakteristische Zahlen (Bigenwerte) mit den Qua- 
draten der Eigenfrequenzen des vorgelegten Schwingungssystems übereinstimmen. Für homogene Ma- 
schinen werden alle Eigenfrequenzen explizit, für teilhomogene Maschinen implizit mit Hilfe der Frequenz- 
funktionen angegeben, für inhomogene Maschinen werden die Eigenfrequenzen nach konvergenten Matrizen- 
verfahren numerisch ermittelt. 


The dynamical matriz is set up, the characteristic values of which are the squares of the natural frequen- 
cies of the given oscillating system. For homogenous machines, the natural frequencies are altogether given 
ewplicitly, for partly homogenous machines, ihey are given implicitly by aid of frequency functions, and for 
non-homogenous machines, they are found by converging matrix methods. 

On definit la matrice dynamique, dont les valeurs propres sont les carr&s des frequences naturelles d’un 
sysieme oscillant donne. Pour des machines homogenes toutes les frequences naturelles sont calculees expli- 
citement; pour des machines homogenes par parties elles sont donnees implicitement a l’aide des fonctions 
de frequence. Enfin, pour des machines inhomogenes, elles sont irouvees numeriquement par des methodes 
convergeantes & malrices. 

Onpepenaetca ‚‚AUHAMMUyecuaa‘‘ MATPHNA, xapakTepucTuyeckue (COÖCTBEHHEIE) SHAYEHUA 
KOTOPOH COOTBETCTBYET KBANPATAM COÖCTPEHHEIX YACTOT MAHHOK KOlebaTeIbHOH CHCTEME. 
B cuyyae ONHOPOAHEIX MAIIHH BCe COÖCTBEHHBIE YACTOTEI BEIPAIKAWTCH B ABHOM BUN, B 
cHyyae YACTHYHO ONHOPOAHEIX MAIIHH —B HeABHOH dopMe IpH NOMOIMH PyHunNuä yacToTsı, 
& B CAHY4ae HEONHOPOAHEIX MAIHHH COÖCTBEHHLIE YACTOTEI HAXONATCH YHCJHIEHHEIM IIYTeM IPu 
HOMOIIH CXOAAIMeTOCH MATPHYHOTO METONa. 


1. Einleitung “ 


Es gibt eine große Anzahl von Verfahren zur Berechnung der Torsionsschwingungen von 
Maschinenwellen. K. Klotter!) hat alle bisher bekannt gewordenen Verfahren einer gründ- 
lichen Analyse unterzogen. In seinem Bericht findet sich keine Methode, die von der Matrizen- 
rechnung Gebrauch macht. Da sich die Eigenfrequenzquadrate von ungedämpften Torsions- 
schwingungen als Eigenwerte einer Matrix darstellen lassen, liegt es nahe zu fragen, ob die Itera- 
tionsverfahren zur Berechnung der Maximalwurzel einer Matrix auch Verwendung finden können 
für die Berechnung der Eigenfrequenzen. Dem steht zunächst die Schwierigkeit entgegen, daß 
in der Praxis gar nicht die größten, sondern nur einige der kleinsten Eigenfrequenzen interessieren. 
Also müßte vor Anwendung eines Matrizenverfahrens zur Berechnung des größten Eigenwertes 


_ zuerst der Übergang zur reziproken Matrix geschaffen werden, und das ist im allgemeinen bei 


Matrizen höherer Ordnung lästig. Man kann jedoch den kleinsten Eigenwert auch ohne Umkeh- 
rung der Matrix dadurch gewinnen, daß man zweimal das Iterationsverfahren zur Berechnung der 
Maximalwurzel durchführt.?) A 


1) K, Klotter: Analyse der verschiedenen Verfahren zur Berechnung der Torsionseigenschwin- 


gungen von Maschinenwellen. Ingenieur-Arch. Bd. 17 (1949), 8.1. 
2) A. Vogel: Zur Bestimmung der Eigenwerte einer Matrix durch Iteration. Z. angew. Math. 


Mech. Bd. 30 (1950), 8. 174. 
; ; 2 Fi 
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Die bekannten Matrizenverfahren können I maßen auf homogene, Yirıhafrhgene und 28 
inhomogene Maschinen angewandt werden, haben jedoch nur für inhomogene Maschinen Bedeu- 
tung, w eilsich für homogene Maschinen die Eigenfrequenzen explizit anschreiben lassen, während EN 
man für Maschinen mit homogenem Hauptteil und nur wenig Zusatzdrehmassen nach der von 
R. Grammel?) angegebenen Methode der tabulierten Frequenzfunktionen bequem alle Eigen- 
frequenzen bestimmen kann. 

In 2. wird die dynamische Matrix (Frequenzmatrix) für ungedämpfte Torsionseigen- 
schwingungen von Maschinenwellen mit beliebig besetzten Drehmassen und beliebigen. Dreh- 
N, steifigkeiten aufgestellt. Im Falle der homogenen Maschine werden alle Eigenfrequenzen explizit 
angegeben (3.). Für Maschinen mit homogenem Hauptteil und einzelnen Zusatzdrehmassen wird 
die Frequenzfunktion nach einem Satz von J. J. Sylvester mit Hilfe der Frequenzfunktionen 


2 ae 
OSEIRUURGHTIEREN ni 


x für den homogenen Teil der Maschine dargestellt (4.). Schließlich werden für eine mit fünf be- b 
b: liebig verteilten Drehmassen arbeitende inhomogene Maschine sämtliche Eigenfrequenzen nach 6 hi 
2 einem Iterationsverfahren berechnet (5.). 7 
1 2. Die Frequenzmatrix e 
Re Die auf ihre Torsionsschwingungen zu untersuchende Kurbelwelle mit ihren Arbeitszylin- h 
Br dern und sonstigen Drehmassen seiin bekannter Weise®) auf eine masselose glatte Welle mit n + 1 
: j einfachen Scheiben von kreisförmigem Querschnitt abgebildet. 
et} * l 
y  Mı yet yr-ı Mk Yk j Yn-ı Mn Yn : 
C ". Ok r . On ; 
I, Ik RO m | - 
9% 9 Or-ı [77 On-ı On 
0, (k=0,1,...,n)sei das polare Trägheitsmoment der %-ten Scheibe, kurz k-te Drehmasse 
genannt. Die Länge des k-ten Wellenstücks zwischen der (k — 1) -ten und k-ten Drehmasse sei /;, 
die Torsionssteifigkeit eines solchen Wellenstücks sei C,(k=1,2,...,n). Bezeichnen wir noch 
die über den stationären Lauf der Maschine überlagerten Schwingungsausschläge der Drehmassen 
in ’bezug auf eine beliebige, raumfeste Richtung mit y, (k=0,1,...,n) und die sie hervorrufen- 
den Torsionsmomente mit M,(k=1,2,...,n), so lauten ph 
a) die statischen Gleichungen 3 
28 Or f 
a) =, (m — wm) Ik 4; 2 BR n) Eu tale ehe m (A), £ 
b) die dynamischen Gleichungen 5 RE e 
9,9 = Mi — Mı+1 ni | so 
Ei Ar Sun Take ee | 
Mit Hilfe der „‚Verdrehwinkel“®) Be BE N, as Naeh 
- jäßt sich die Anzahl der Bewegungsgleichungen um eine Ba Aus Ur wird En ONE | 
I; | a 
nun M;. (k=1, RL GR n) ae Elan ” 
aus (2) wird | A RE Ei 
} 5 . M.-ı ( 1 1 ) Mi+ı a RT ee 
s ‚= MM — Mm — — + —|M,;, — k=1, 2 Jan): Er 
a A if ! 9-1 = 9;-ı iR 9: \ 9; ( ei ) s\ ) i ET si 
’ Wenn wir jetzt die Matrizen & a N 
ß I, 
Yı | C, M, 
® iz ’ st 3 ’ IM Tr ’ 
ER ) RE AR 


®») R. Grammel: Ein neues Verfahren zur Bere 
wellen. Ingenieur-Arch. Bd. 2 (1931), S. 228. | j ! 

#0. B. Biezeno— R. Grammel: Technise 

5) T. PöschlundL. Br: Über die Bere 
 gener Maschineuinäit un m . Z. angew. Math. 


v 


ae 


4 ’ 2. s 
A 
MN Pu R wo 
ve KR u a 


Ks Dr  .. Bent 0 \ 
a je ı = ER ER 9, n ” 
“ einführen, so lassen sich die Gleichungssysteme (3) und (4) auf die einfache Form I 
# chr) 3% 5 DL UNE. a. (0) 
bringen. ‚Nun kann man M eliminieren und entweder nach ® oder ® auflösen. Im ersten Fall Ai n 

| erhält man aus (5) und (6) Ri 

3. =—LT1®, | I 
im zweiten Fall ö 

RER ER Sk LET KR NEE Sr ER (2) / 


Die er der Matrix 3 ist umständlich, dagegen läßt sich die DR ix $ einfach 
dadurch umkehren, daß man die Elemente in der Hauptdiagonale umkehrt. Wir setzen 


' TmI=9 Ai 
{ und nennen D die „dynamische Matrix‘. Sie hat die Form 
[° r 6 203 \ 
F (ve % Ge, —6yz 
| Ba | NS 
Er Ä ; 
; ; \ a‘ Ca 
wenn man die Abkürzungen 

> 6 ErHbeN, und“ ao A (1,2 n) ! (9) 

zer 2 ; ENTF ER SE RE) Er ee TR eh, 


BR verwendet. Die Koeffizienten c, und cz heißen die Wellenkoeffizienten der Maschine. 
FE ER Mit ® wird aus Mm A 


N nV \ | ö ea, DE RR Rd (10). 
act man für Y% Ken bei stehenden Schwingungen üblichen Ansatz 


En ar ni a at 
At, mit der Belehendene d,so wird E | 
2 DENDE We} N “ A = —ag 
BR EEE NE er EN N 
t man noch a =h so wird aus (10) und (u) 
a D0d=10 
. i Kae: a 
Br ne— D)B-9. A 


ende unserer Schwingungen sind also identisch mit den Eigenwerten 
en Matrix ». a der unsyanglllehen Matrix © kann man auch die symme- 


366 


nehmen, die aus ® durch die Ähnlichkeitstranstormation Dee R-19 SB mit ser nicht singulären. 


Transformationsmatrix 


4 


ER ee 
n ; G;; 


hervorgeht. In def dynamischen Matrix D oder ersalzweise in der symmetrischen Matrix & 
treten positive und negative Elemente auf. Für die praktische Rechnung ist es bequemer, wenn 
man nicht auf die Vorzeichen zu achten braucht. Wir ersetzen deshalb ® bzw. © durch die 
ähnliche Matrix 


\ 


. . J ’ VETER R { 
ae | Re +Jeie, 2 NA 
RR HY 6, Te 0, +-Ye, DR Be 
G- ar bzw. 3 = 6 \. 142), 
— in ı ar cn | Be Vaıch | nt Ca 
die aus D bzw. © vermöge der Ähnlichkeitstransformation mit der Matrix 
PAST RT, 
+1 | | N} 
Ku BR | 

| : (— 1% 


entsteht und keine negativen Elemente mehr enthält. Die Matrix & bzw. 3, deren Eigenwerte 4; 


die aa ö;liefern, nennen wir die „Frequenzmatrix‘, ihre charakteristische Funk- 


tion . \ a)=|raE—E| . SER. (13). 


- die „Frequenzfunktion‘ des Systems, dessen freie TOR Durch erden 
sollen. ! hi 5 


3. Die homogene Maschine 
Unter einer homogenen Maschine versteht man eine Maschine mit gleichen Drehmassen ©, 


gleichen Wellenstücken von der Länge I; und gleichen Drehsteifigkeiten C;. Es ist also 


%=G%—=c(k=1,2,...,n). Die Frequenzmatrix der homogenen Maschze ‚mit N En 1 Dreh- 
massen ist nach (12) 

LANG Vrasne 

er ee Dt BERN 


x KERN, 


“ 


Die Frequenzfunktion ist nach (13) ir \ | 
| | Male s. EERRES, 


Diese spezielle Determinante ist eine Kontinuante, Die re Ko 
_ durch die 2n + Ü)- u Kontinuante 


ERENENHERNENEENEERER 


a 


Führt man statt A die a Größe © = . oz so entsteht die reduzierte Frequenzlunktion 


D 


ae 


u ee 


a a ie u Pu . ı A da A}, 
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ersetzen®), die sich ihrerseits als Tschebyscheffsches Polynom zweiter Art in der Form 
sin(2n--1)r 


} a. & === 
| darstellen läßt‘), wenn man Ve= cos rsetzt. Die Nullstellen a =1,2,...,n) von 9,(£) er- 
geben sich aus den Nullstellen des Tschebyscheffschen Polynoms zu | 
«(n) \ 5 y IT £ 
En) 14 le E et 
ar cos ey (3 { n) 
und damit die Figenfrequenzen @, der homogenen Maschine zu 
Bun 2 V ccos \ — = se ne urn). 


Dieses Ergebnis ist bekannt. 


4. Die Maschine mit homogenem Hauptteil und einzelnen Zusatzdrehmassen 


Wenn zu den n gleichen Wellenstücken mit »--1 Zylindern einer homogenen Maschine 
weitere Drehmassen, z. B. Schwungrad, Dynamorotor, Kupplung usw. hinzutreten, spricht man 
von einem Maschinenaggregat mit homogenem Hauptteil und einzelnen Zusatzdrehmassen. 

' Wir betrachten beispielsweise ein solches Maschinenaggregat mit n gleichen Wellenstücken 
(mit a + I Zylindern), einer zusätzlichen Drehmasse links davon und zwei zusätzlichen Dreh- 
massen rechts davon. 


Cnr+1 n+3 
n Wellenstücke 
Die Frequenzmatrix dieses Systems ist 
, 
[fr r En+1 
LG 
| En+1 
C= 
En+2 
[4 
Cn+2 t On+2 Cn+3 
‘ ! 

On+2 " Cn+3 Sr Cn+3 


die Frequenzfunktion ist /n+3 (A) = | & —E|, die reduzierte Frequenäfunktion 


mtr) —1 


u u a Fe en 


—y [2 —1 h 
2 
er EEE 
| Ze Be 1 (15) 
4 5 ü all —2 SEE | 
E ® —1{—(»n+r)  —% 
| u a a) 
mit den Abkürzungen 
er mi ana, 3). 
G G G 


Da man die Eigenfrequenzen einer homogenen Maschine kennt, wird man bestrebt sein, 

_ möglichst viel von dieser Kenntnis zu verwerten, wenn es darum geht, die Eigenfrequenzen einer 
teilhomogenen Maschine zu bestimmen. R. Grammel’) hat die Frequenzfunktionen 9,(C) für 

homogene Maschinen tabuliert und.die reduzierten Frequenzfunktionen 9,47 (£) für teilhomogene 

Maschinen mit p Zusatzdrehmassen mit Hilfe von ,(£) und @,-1(&) für homogene Maschinen 


Tu 6) A. Vogel: Der Sättigungswert. des Kurzschlußstroms in einem Netz mit beliebig vielen Kraft- 
 _ werken. Z. angew. Math. Mech. Bd. 28 (1948), S. 218. 

ar ”)R. Grammel:a.a.O,., Eee 8) undR. Grammel: Die Berechnung der Drehschwingungen 
_ von Kurbelwellen mittels der Frequenzfunktionentafel. Ingenieur-Arch. Bd. 3 (1932), 8.277. 


- / 


en n heiten L Ingenieur-Arch, Bu. 14 (1943) 


"Funktionieren der Iterationsverfahren garantieren. Es ist nicht zu befürchten, daß Mehrfach- 
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dargestellt. Der Zusammenhang zwischen Pn+2 (2): Pn(E) und BR = eur sich sofort aus 
einer bekannten Eigenschaft von Kontinuanten®): 


ml 
N," 
ir! 
Mm; | l; 
N; My+ı li+ı 
Hera? 
i ol 
i FSK I 
m Ih Mi+ı lit m; Iı Mi+z li+2 
a Dat 10, vr er 
ee ; I; N; i 
A ie N en, Be 
Ni MM; N My Na Mi) Ny—ı My 2 


Nach diesem Satz wird die Frequenzfunktion (15) unserer teilhomogenen | 


In+3(d)= EAnbR): ee v2) ER Ka le + )) Ze) | 


RR Ze 9 

mr ‚ +92 [E- (nt Y)IIE-Ptr)) vr - DE (rsFr)PR- ld) 
\ NEE) 

Da ssich die Funktionen 9, (£) und 9,-,(£) nach (14) als Tschebyscheffsche Polynome explizit 

darstellen lassen, hatman jetzt auch für 9,+3(£) eine explizite Darstellungin ADBOIAESTE von 

allein. 


Die gesuchten reduzierten Eigenwerte Z, als Nullstellen der reduzierten Frequenzfunktion 
Pn+3(£) genügen einer transzendenten Gleichung von der Form f (£)9, (2) = 9(£) 9n-ı (2). Daraus 
ergeben sich unter Verwendung der Tafeln für 9, (2) die Nullstellen selbst dadurch, daß man die 
Abszissen der Schnittpunkte der beiden Kurven 9, = f(£) 9„(&) und y,= 9(£) „1 (£) bestimmt. 

Wenn nun zu den n--1 gleichartigen Drehmassen des homogenen Häuptteils einer Ma- 
schine noch zahlreiche Zusatzdrehmassen hinzutreten, leidet natürlich die Übersichtlichkeit des 
sonst sehr bequemen Verfahrens der Frequenzfunktionen. Man kann dann ein automatisch 
arbeitendes Matrizenverfahren anwenden, das die Eigenfrequenzen numerisch zu bestimmen er- 
möglicht. In der nächsten Nummer wird die Berechnung der Eigenfrequenzen nach einem ‚solchen 
Verfahren für eine inhomogene Maschine durchgeführt. 


ne re in ag in 


‚5. Die inhomogene Maschine 


Wenn man die Eigenwerte A; der dynamischen Matrix ® oder der dazu ähnliähen Matrizen © 
bzw. & oder 3 kennt, berechnen sich die gesuchten Eigenfrequenzen der Maschine aus 


Be /Fos BO i 


Die Berechnung der Eigenwerte A; kann nach den früher beschriebenen "Iterationsverfahren®) 
erfolgen. Man weiß, daß alle Eigenfrequenzen reell, positiv und verschieden sind!P). Die dyna- 
mische Matrix ® oder die Frequenzmatrix € weist also genau die Eigenschaften auf, dieeinglaties 


wurzeln, komplexe Wurzeln oder zwei reelle Wurzeln mit gleichem absolutem Betrag auftreten, 
die eine Modifikation der Iterationsverfahren erforderlich machen würden. 

Wenn es nur auf die Bestimmung der Eigenfrequenzen ankommt, nicht auf [die chwingungs- & 
form, verwendet man zweckmäßig die Matrix & mit lauter positiven Elemente bzw. die dazı 
ähnliche symmetrische Matrix 3. Man braucht dann bei der Rechnung ga; nicht auf die Vor ® 


an 


zeichen zu achten. Möchte man aber außer der Eigenfrequenz auch noch ie Schwingı gsiorm 
kennen, verwendet man zweckmäßig die 2 DO bzv symme tische Form ©. 


Eye 


Ser. (4), Bd. 6 (1853), 8. 297. 
®, A. Vogel: a.a.0. Fußnote 2), 
10) Über diese und weitere Eigenschaften der Eige 


°) J. J. Sylvester: Ona fundamental rule i AlenrBE 
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N \ 3 
In allen Fällen erhält man durch die beschriebenen Iterationsverfahren zunächst den 
‚größten Eigenwert /,, der praktisch allerdings am wenigsten interessiert. Zum kleinsten 
Eigenwert A, könnte man dadurch gelangen, daß man die reziproke Matrix zugrunde legt. Dieses 
Verfahren ist jedoch nicht zu empfehlen, weil erstens die Umkehrung einer Matrix von hoher 
Ordnung recht umständlich und zweitens die reziproke Matrix eine Vollmatrix ist, also nicht mehr 
i die für das Iterationsverfahren recht günstige Struktur der Frequenzmatrix aufweist, die eine 
Dreiermatrix ist, wodurch die Rechenarbeit für jeden Iterationsschritt erheblich erleichtert wird. 
Wir berechnen deshalb lieber erst den größten Eigenwert A, von € und dann nacheinander die 
größten Eigenwerte 05; On-1,;,°*' der Matrix R=(A,©— €). Die kleinsten Eigenwerte von & 

sind dann 

An) =4 


An-ı 2: 4, — On-1 


Beispiel: Es sollen die Drehschwingungszahlen der durch die Figur beschriebenen Ma- 
schine ermittelt werden. 


©; 1,5 . 1010 1,2.10°° | 0,6.102%. |o,3. 10%0 
1, 50 | 20 | 
9; 120 40 20 30 10 


In die Figur sind die Zahlenwerte für die fünf Drehmassen ©; [cm kgsec?], die Längen /; [cm] 
‚und die Drehsteifigkeiten C, [em?kg] eingetragen. Damit können die Wellenkoeffizienten der 
Maschine berechnet werden. Nach (9) ist 

ee 0,25» 10° = 0,75. 107: 

u a BR LE RN VA 1185 

EN UN Re EA EHE RE LEER 

0 al,9,,5107: ie Zeil 10T 


Die Frequenzmatrix € ist nach (12) aufzustellen: 


1 1 
«| 975 : a -107. 
0002) 


Die Frequenzfunktion ist fı(A)= |AE—€|. Um dimensionslose Elemente zu erhalten, divi- 
dieren wir & durch er +ej= 1-10" und bekommen so die reduzierte Frequenzmatrix 


| (! { '% 1 0,866.025 
BU 0,75 : 0,6 oder die sym- | 0,866025 3 1,095 445 
I=-—_ = R Z= 
10? 2 1 0,5 |metrische Matrix 1,095445 1 0,447 214 
a | 0,447214 2 


und dazu die reduzierte Frequenzfunktion F,(£) = |E& — Z) mit &= 10-74. Die Nullstellen £; 
von F,(£)sind die Eigenwerte von Z. Mit Hilfe der Z, ergeben sich die gesuchten Eigenfrequenzen 
der Maschine zu | 
6. = = ZU 100 Cr (k= 1; 2: 3, 4) Er LI (16). 
Man bestimmt zunächst die Maximalwurzel &, von Z und dann die beiden größten Eigenwerte o, 
und 0, von P=(£, € —Z). Die vollständige Rechnung ist‘in der mehrfach zitierten Arbeit!!) 
für die oben vorliegende Matrix Z durchgeführt worden. Dort ergab sich £, = 3,732 891; 
| 04 3,508 773; 05 2,7% 801. Damit wird 
. .  b=ıb-—- 00,224 118 
- } GL — 03 = 0,96 09% . 


4 
Wegen z &=7 muß schließlich der noch fehlende Eigenwert &,—=2,106901 sein. Die direkte 


Bene dieses zweitgrößten Eigenwertes ergab den Wert 2,106 900. 
Buch (16) sind die gesuchten a der Maschine 


i S = 1497 sec7!; = 4590 sec-! 

he ü; = 3060 sec-!; = 6110 secr!, 
TR Vogel. 2.2.0. Fußnote 2), R 
en, am KR ‚10. 1999. | Pe u 
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Berechnung an den Rändern belasteter, allgemeiner Schalen 
f Von W. Zerna in Hannover 


Nach einer kurzen Zusammenfassung der Ergebnisse der allgemeinen Theorie dünner elastischer Schalen 
in der von Green und Zerna entwickelten Form wird der Fall der an den Rändern belasteten, beliebigen 
Schale behandelt. Das Problem wird auf eine einzige komplexe Differentialgleichung vierter Ordnung 
zurückgeführt. Der Realteil der Lösungsfunktion liefert eine mit den Normalverschiebungen zusammen- 
hängende Größe, der Imaginärteil eine Spannungsfunktion. Mittels dieser beiden Größen lassen sich die übri- 
gen Verschiebungskomponenten, sowie die Längskräfte, Querkräfte und Momente der Schale leicht berechnen. 

After a brief summary concerning the results of the general theory of ihin elastic shells, as it has been 
developed by Green and Zerna, Ihe author discusses the case of the shell of arbitrary shape that is loaded 
at its boundary. The problem is reduced to a single complex differential equalion of the fourth order. The 
real part of the solution is connected with the normal displacement, the imaginary part wülh a siress 
function. By aid of these two quantities, the other components of the displacements as well as the longi- 
tudinal and transversal forces and ihe moments of the shell can be easily comphted. 

Apres avoir retrace les r£sultats de la theorie gen£rale des nappes minces &lastiques, yui a ElE dEvelopp£e 
par Green et Zerna, l’auteur traite le cas de nappes d’une forme arbitraire et qui sont chargees au bord 
seulement. Le probleme est reduit a une seule Equation differentielle complexe de quatrieme ordre. La part 
reelle de la solution fournit d’une grandeur, qui est en connexion avec les deplacements normals, la part 
imaginaire d’une fonetion de tension. A l’aide de cettes deux grandeurs on calcule les autres composantes 
du deplacement, les forces longitudinales et transversales et les moments de la nappe. 

Ilocııe kpaTKoro ONHCAHUA Pe3yAIBTATOB TEOPUH TOHKUX YHPYTUX 000JI0Yek, paspa60TaHHoH 
I'punom u Ilepna, uecsenyerca sarpyskeHHaag TO KOHTypy O00N10yka „m000ü hopMBL. ITa 
3amaya CBOAUTCHA K OMHOMY KOMILNIeEECHOMy AudepeHNmalsHoMy YPaBHeHuWw yeTBepToro 
mopanka. JleiiersnteipHan yacTb PeIeHHA 9TOTO YPaBHeHHA NAeT BEeIHYHHY, CBABAHHYI C 
HOPMA.IBHOÜ KOMIIOHEHTOH CMEINeHHA, 3 MHHMAA YACTb — PyHRINO Haupsskenua. Ilpn momomm 
9TUX ABYX BEJIHYUH MOTYT ÖBITb JIETKO BEIYHCHEeHBI APYTHe KOMIIOHEHTEI C MeIICHHA, H TAKIKE 
IPONONBHEIe HU MOITePeYHLIe CHIIEI H HSTHÖAIIIUEe MOMEHTEI ODONOYEN. 


1. Einleitung 

Die Berechnung von elastischen Schalen läßt sich in den meisten praktischen Anwendungs- 
fällen in zwei wesentlichen Schritten durchführen. Zunächst wird der Membranspannungs- 
zustand ermittelt, der im allgemeinen jedoch nicht die vorgegebenen Bedingungen an den 
Rändern befriedigt. Die von den Rändern ausgehenden Störungen des Membranspannungs- 
zustandes werden dann durch die Biegetheorie der nur durch Randkräfte bzw. Randmomente 
belasteten Schale erfaßt. Während sich der Membranspannungszustand in vielen Fällen recht 
einfach bestimmen läßt, bereitet die Ermittlung der Biegespannungen nach der üblichen Theorie!) 
fast stets erhebliche Schwierigkeiten. Es ist der Zweck der vorliegenden Arbeit, einen ein- 
facheren, allgemeinen Weg zur Berechnung der an den Rändern belasteten, beliebigen Schalen 
aufzuzeigen. Den Ausgangspunkt dazu bildet die von Green und Zerna?) aufgestellte 
Theorie dünner elastischer Schalen, die eine systematische erste Approximation des Schalen- 
problems darstellt, und deren Ergebnisse zunächst dargelegt werden sollen. 

Es werden völlig allgemeine Koordinaten zugrunde gelegt und für die Herleitungen die 
Vektor- und Tensorrechnung in der Schreibweise des Ricei-Kalküls benutzt. 

Demgemäß werde verabredet, daß griechische Buchstaben als Indizes die Werte 1, 2, 3 
und lateinische Buchstaben die Werte 1, 2 annehmen. Sofern nicht ausdrücklich anders ver- 
merkt, gilt die übliche Summationsregel. daß über doppelt auftretende gleiche Indizes zu sum- 
mieren ist. Partielle Differentiation soll durch ein Komma, kovariante Differentiation einer 
Tensorgröße durch einen senkrechten Strich angezeigt werden. Es ist beispielsweise 

Odyx 


y 2 2 
= 5, AB = A, — IQ, , Aypiy = Aapıy — I sy ap — I pyQaı ; 


Ax,B 


!) Die übliche Schalenbiegetheorie baut auf den grundlegenden Annahmen auf, daß die Normalen 
zur unverformten Mittelfläche auch Normalen zur verformten Mittelfläche, und daß die Spannungskompo - 
nenten normal zur Mittelfläche vernachlässigbar sind. Mit diesen Annahmen wurden allgemeine Schalen - 
' biegetheorien, die völlig allgemeine Koordinaten verwenden, entwickelt von A. L. Goldenweiser: 
Akad. Nauk. S. S. S. R. Prikl. Mat. Mek, 8 (1945) S.3; J. N. Rabotnov: C. R. (Doklady) Acad. Sci. 
URSS (N. S.) 47 (1945) S. 87; W. Zerna: Ingenieur-Arch, 17 (1949) S.17;H. Neuber: Z. angew. Math. 
Mech. 29 (1949) 8. 98 und 142, 

In der letztgenannten Arbeit von Neuber wird eine abweichende Darstellung benutzt. Es werden 
zwar auch noch die oben erwähnten grundlegenden Annahmen verwendet, aber einige Besonderheiten ein- 
geführt. So werden vor allem die Grundgleichungen nicht durch die üblichen Schalenschnittkräfte und 
-momente, sondern in anderer Weise ausgedrückt. Der dort eingeschlagene Weg deckt sich teilweise mit einer 
von W.Z.Chien: Quart. appl. Math. 1 und 2 (1944) S. 297, 43, 120 gegebenen Darstellung, die allerdings 
wesentlich allgemeiner ist, da dort keine der grundlegenden Annahmen vorausgesetzt werden. Doch wie 
A.E.GreenundW.Zerna: Quart. J. Mech. appl. Math. (im Druck) ausgeführt haben, scheint diese Art 
der Behandlung des Problems keinerlei Vorteile für praktische Anwendungen zu bringen, sondern gibt darüber 
hinaus zu einigen Zweifeln Anlaß. 

®)A. E. Green und W. Zerna: Quart. J. Mech. appl. Math. (im Druck). 


* 
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worin [’5,,dieChristoffelschen Dreizeigersymbole zweiter Art und 9, unabhängige Para- 
meter sind. 
Auch die Kronecker Deltas haben die übliche Bedeutung: 


K=lli=h, K=0h+h). 


2. Allgemeine Schalentheorie 


Die Punkte der unverformten Schale mögen durch allgemeine krummlinige Koordinaten &, 
beschrieben werden. Die Fläche, die durch &, = 0 bestimmt ist, heißt Mittelfläche; x, ist der 
Abstand eines Punktes von der Mittelfläche und wenn ? die Dicke der Schale angibt, so ist 


t 37 
Die Schalenlaibungen werden durch die Gleichungen x, = +5 bestimmt. Die Parameter &; 


stellen irgendein allgemeines Koordinatensystem auf der Mittelfläche dar. 
Es werden nun dimensionslose Koordinaten eingeführt, indem gesetzt wird 


= — u, > = (nicht über © summieren!), 


worin L, für die Längenabmessungen der Schale charakteristische Konstanten sind, die derart 
gewählt werden mögen, daß 


REdUE N: 
worin das Intervall ($&;, 7;) gleich der Einheit ist. Auch d, läuft gemäß (2.1) im Intervall Eins. 
Der Ortsvektor R* eines Punktes der Schale werde durch 


Rr=ır +08 h 
angegeben, worin r* der Ortsvektor der Mittelfläche und e, der Einheitsvektor der Normalen- 
richtung zur Mittelfläche ist. Dimensionslose Vektoren R, t werden nun durch 


RD en | 
definiert, worin L eine charakteristische Länge angibt, die unter den Z, und den Krümmungs- 
radien der Mittelfläche die kleinste Größe darstellt. Wird 


t 
Bez 
gesetzt, so bestimmen sich die Punkte der Schale durch den dimensionslosen Ausdruck 
R =Y +19, eg TU are ee ea heil oa (2.2). 


Von nun ab werden nur noch R und r benutzt, die als reduzierte Ortsvektoren bezeichnet werden 
mögen. Weiterhin werde angenommen, daß A =const., obwohl es keine grundsätzlichen 
Schwierigkeiten bereiten würde, veränderliches A zu berücksichtigen. 

Der kovariante bzw. kontravariante Maßtensor der reduzierten Mittelfläche werde mit 
gi» bzw. g‘* bezeichnet. Die Koeffizienten der zweiten Grundform seien b;z; sie bilden die kovari- 
anten Komponenten des Krümmungstensors, dessen gemischte Komponenten bj; sich zu 

br = g'! bir 
errechnen. Die mittlere Krümmung H und das Gaußsche Krümmungsmaß K lassen sich 
in der Form 


1a, | 
=55, K=b! 2 —bı b2 


schreiben. Alle hier angegebenen geometrischen Größen beziehen sich auf die reduzierte Mittel- 
fläche und sind demgemäß dimensionslos. 

Für die in (2.2) gegebene Darstellung der Schale ist von A.E. Green undW.Zerna?) 
‘eine Schalentheorie entwickelt worden, in der nicht von den üblichen Annahmen!) aus- 
gegangen wird, sondern die sich auf die Voraussetzung Arünget, daß die Dicke der Schale hin- 


reichend klein ist, so daß 
IM ...32) 1, SSEB7 ER FREE Ber (2.3). 


Eine damit Anighehte systematische Approximation führt auf die nachstehend ee 
Pigebüisie: 


27* 


ae 
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Wird mit s*? der a varianie Spannungstensor bezeichnet und die Spannungsgröße 
ok (5 — 20 6) ser, 
sowie der geometrische Ausdruck 
m =A(l—2A0H +2°02K) 


eingeführt, so lassen sich zweckmäßig folgende Schnittgrößen definieren, die in enger Beziehung R 
zu den physikalischen Längskräften, Momenten und Querkräften stehen. % 


Rn Er REN 3% 
2 | Nik— [ moikdd®, Ä 5 
' — 1/2 , 3 
j R 
Ft +1/2 | 
N | Mir — | md. 0% do, een... .: 
2 ne. e ® | n Z 
nV +1/2 
AN Qi = [ msidds. 
I —ı/2 f 
R | Ferner werden die Ausdrücke 
Br malt | 

; RR a (2.5) 

R* -/md Ra 


verwendet, die durch die auf die Schalenlaibungen wirkenden Belastungen bestimmt sind. 
Be: Unter Benutzung der Größen (2.4) und (2.5) nehmen die Gleichgewichtsbedingungen am 
, Schalenelement folgende Form an 


SE u BR >: 0 % 
Nr, —abEQ—abERr+ Pr =0,] . (2.6), 


Nr!b, +AQr|, +APS=0 
wg ibm mo 0. 
Die kovarianten Differentiationen darin bedeuten 
Nrkl,— Nrk,, + T*,Nriı Ti, Net, | / 

Mr, — Mk, + Er Mrı LP, MEN ee R 
% 
Me En ze 5 QJ | I 
worin die Christoffel-Symbole T‘}.ı mittels der Metrik der Mittelfläche zu berechnen sind. 2% 
Die Gleichungen (2.6) und (2.7) stellen das Gleichgewicht der Längskräfte und Momente E 
dar und sind exakt im Sinne der klassischen Elastizitätstheorie. Die Näherung kommt nun erst. - 78 
bei Aufstellung des Elastizitätsgesetzes herein. H u 


0 Die Verschiebung eines beliebigen Punktes der durch (2.2) gegebenen Schale sei durch die 0034 
=  dimensionslosen Vektorkomponenten % = %y (9 ®, 9,) bestimmt, aus Be die flgenden 


“ Verschiebungsgrößen gebildet werden. x BE 3 
+1/2 +1/2 3 +1/2 % ir = r 
n= [ vwd», W= [| ud, W== | (1493) 0,40 es). 
8 —1/2 2 1/2 


Be 
Es kann vorausgesetzt werden, daß V, und w, etwa von der gleichen Größen erden = 
sind. Dann läßt sich mit (2.3) das Elastizitätsgesetz näherungsweise in die A Br FAN: N 


2 ve 
ve — ER ge (Rr4 25 Pi) = DER 
| 22 u 
ER AT RE y Y 
RE im) er 
1 y_5D | 


_ Darin bedeutet ie die, Querkontraktionsziffer und E ge 


Math. ‚dien, 1950 Zerna, Berechnung an den Rändern belasteter, allgemeiner Schalen 373 


= sung 
Be a") | 


— 5 [g {rat + ge! gi’ + u(et! gkr + gkl e'r)] | 


Der Tensor vierter Ordnung Eik! ist gegeben durch 


\ Bikir — Br ge 


mit 


1 
sell = 23V, sl? — — el —, 9= |gir| = 911 922 — Yı2 9ı2- 


Ar 


B Es gelten die Symmetrieeigenschaften 

= Eiklr — Hkilr — Hikrl _ Hlrik, 

Die Ausdrücke (2.10) bis (2.12) stellen den Zusammenhang dar zwischen den in (2.4) definierten 
' Schnittgrößen und den in (2.9) eingeführten Verschiebungsgrößen. Zusammen mit den Gleich- 
gewichtsbedingungen (2.6) und (2.7) ergeben sie das vollständige System der Gleichungen zur 
Lösung des Schalenproblems. 


3. Schale mit Randbelastungen 


Es werde nun der Fall betrachtet, daß die Schale nur an den Rändern durch Kräfte oder 
Momente beansprucht wird, im übrigen aber unbelastet ist. Wird W; aus (2.11) mittels (2.12) 
eliminiert, so ergibt sich mit P* = R*=0 
b y 2 
® Mik _ ars [gli Mrkl,, + u(gi® Mril,, — gi Mti],,)] = — BEiku P,,, . . (81), 
worin. eins | Er 
} Di EA: . 


PATE 
gesetzt ist. Der Ausdruck (2.10) lautet jetzt 
Se \ NEID En N NE et (9:2); 
wenn 
Ylr r—— Yrl:— Fish Vu — by1 V, ANETTE a LTR 0 (3.3) 


j Me bedeutet. 
RE Wird (3.2) nach y;, aufgelöst, so ergibt sich 


ya 4 Frur TE Ren N N LA Re RR (3.4), 
worin 
£ Fir = 5 [9:1 9er + Irı Gier — MlEir net em ir)) » > + (83) 
mit 
1 = 89 =, 22 = — &1 =V9 
Es werde der Ausdruck 
Met gih Yaıt — an NEN N NEE NR (3.6) 


2; eiet. Bei Einführung von (3.3) und Beachtung von (2.12) zeigt sich, daß mit der Genauig- 
_ keit der hier benutzten Näherung (3.6) identisch verschwindet und geschrieben werden kann 


a a ON a aRlend: 


Mit (3.4) wird daraus, wenn berücksichtigt wa ‚daß die kovariante Ableitung des Maßtensors 
Ul ist: 


di en keöchben Ni Krlen ne durch eine skalare Spannungsfunktion ® 
ktorielle Be K; a. wi Es wird der Ansatz 


EN 


ke, 
5 


Taler 
2 y \ 
fe 2 % 
IV RE L r ) 
Ev Ran De 


a a a a EB A Zu Dt Ne 
A ne 
ar F u 
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gemacht. Einführung in die erste Gleichung (3.9) bringt nach einigen Zwischenrechnungen 
Ta - LP uiı+ 1 Tamm - lin) dr + et ek nr 0 2. let), 


woraus sich schließen läßt, daß %,]„ die Größenordnung von ® hat. 

Es ist bekannt, daß infolge der Randbelastung die Spannungen in der Schale in der Form 
von gedämpften Schwingungen verlaufen. Aus der zweiten Gleichung (3.9) kann nun gefolgert 
werden, daß die betreffenden Funktionen beim Differentieren um A-V/? mal zunehmen. Dann 
wird aus (3.10) und (3.1) in konsequenter Durchführung der eingeführten Näherungen 


Nik — gir ekl Bl, Reh Ve (3.12a, b). 
Werden die Beziehungen 
ek, gel 
berücksichtigt und (3.12a) in (3.8) eingeführt, so ergibt sich 
gi HE Diet PA Ha ee Van 0... 2 (3.13), 


womit eine partielle Differentialgleichung vierter Ordnung für die Spannungsfunktion ® und 
die ‚Verschiebung V, gewonnen ist. 1 


Um nun eine zweite Gleichung zu erhalten, wird (3.12b) in die zweite Gleichung (3.9) 4 
eingesetzt. Bei Berücksichtigung der eingeführten Näherungen liefert dies z 
ABgE gan —buet ed n=0 ....2.2.0.0@14). 3 
Die beiden Differentialgleichungen (3.13) und (3.14) ermöglichen die Bestimmung von ® und V.. h 
Die N‘* folgen dann aus (3.12), die M‘* aus (3.12b), die Q® aus (2.7) und die Verschiebungs- 
größen V; lassen sich aus (3.4) ermitteln. 
Diese beiden Gleichungen lassen sich natürlich auch durch eine einzige ersetzen. Um sie 
zu erhalten, wird zweckmäßig ein Weg beschritten, den E. Reißner?) in einem ähnlichen 
Fall benutzt hat. Die Gleichung (3.13) wird mit einem konstanten Faktor k multipliziert, (3.14) 
wird durch A B dividiert, dann werden beide Gleichungen addiert. Es ergibt sich 
gi gi (Vaikir + k Dru)+ kAE b;ı eir EK(Vzjr% —4A?BEk ®|,x) AR, une (8.15). 
1 ; 
r } a = Re x 
Der Faktor % soll nun so bestimmt werden, daß k = mBER' | 
daraus folgt 3 
k.=iK - 
wenn 
Y12(1—u®) 
| a Be Ve In (3.16) 
Mit Einführung der komplexen Funktion en 
y=VY;tkd=V; 4iKEd...... ARE TEN (3.17) 
und 
E ER, 
BE 7 
geht (3.15) über in \ 
gei Ar ylart bee. (3.18). 


Dies ist eine komplexe Differentialgleichung vierter Ordnung, die das hier betrachtete Problem 
beschreibt. Die kovariante Differentiation braucht dabei mit Rücksicht auf die in der ganzen 
Theorie enthaltenen Näherung natürlich nur bis einschließlich der zweiten Differentialquotienten 


entwickelt zu werden. Vernachlässigung von Gliedern der Größenordnung VA ergibt eine weitere 
Näherung und aus (3.18) wird dann einfach 


gl GEM ikirt 7: EA A) Se 819). R- 
Ist die Funktion 9 gefunden, so liefert gemäß (3.17) ihr Realteil die Verschiebung Z er 
und ihr Imaginärteil die Spannungsfunktion ®. vs . 


®) E.Reißner: J. Math. Physies 25 (1946) S. 80 und 179 behandelt den Sondorfall der flachen 4 ’ R 
Kugelschale und gelangt zu entsprechenden Gleichungen. | 


‘ Eingegangen am 8.12, 1949, L 


} 
- 
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‚Die Plattengleichungen für den elastisch-plastischen Zustand" 


Von W. Swida in Karlsruhe 


Die Grunddifferentialgleichung zur Untersuchung der Plattenbiegung im elastisch-plastischen Zustand 
wird aufgestellt. Für einen Sonderfall, in dem die Biegefunktion durch einen Ansatz angegeben werden 
kann, werden die Spannungsverhältnisse gezeigt. 

For the investigation of the bending of plates in the elastic-plastic state, the governing differential equation 
is stated. In a special case, Ihe bending function is obtainable by way of trial, and the distribution of the 
iensions is explained. 

Pour le probleme de la flewion de plaques dans l’ötat Elastique-plastique L’Equation differentielle fonda- 
mentale est Etablie. Dans un cas special on peut faire une position pour la fonction de flewion ei determiner 
la distribution des tensions. 


BsrBoauTca OCHOBHOE AuhepeHmanpHoe YpaBHeHne AI UCCHEeNOBAHHA U3TUÖA INTACTUHKU B 
yupyro-MIacTuyeckoM COCTOAHUM. 1 HEKOTOPOTO YACTHOTO CAYYAA, IA KOTOPOoro PyHkuus 
u3ru6a MO3KET ÖBITB BhIpaskeHa dopMmyoü, ompenenserca pacıpeperenue Hampaskeunü. 


1. Einleitung 
Es sollen für eine Platte, die elastisch-plastisch deformiert ist (Bild 1), einige Grund- 
beziehungen aufgestellt werden. Dabei werden folgende Annahmen zugrundegelegt: 

a) Die Durchbiegung der Platte ist gering 
gegen die Plattenstärke. 

b) Die Punkte einer zur Mittelebene der Platte 
senkrechten Geraden liegen nach der Form- 
änderung wieder auf einer Geraden und 
bilden eine Normale zur Biegefläche. Diese 
Annahme gilt sowohl für die elastische als 
auch für die plastische Zone. 

c) Die Platte besteht aus einem elastisch- 
ideal-plastischen Werkstoff. Dabei sind die 
Zug- und Druckstreckgrenzen gleich. 


Bild 1 n Bild 2 


2. Plattengleichungen in rechtwinklichen Koordinaten 
Wir betrachten einen beliebigen Punkt K der elastischen oder der plastischen Zone (Bild 2) 
mit den Koordinaten x, y und 2 und bezeichnen die Verschiebung dieses Punktes in der Richtung 
der &- bzw. der y-Achse mit & bzw. mit n. 
Auf Grund der zweiten Annahme ergeben sich für die Verschiebungskomponenten & und 7 


- sowie für die Verzerrungskomponenten &,, &, und y,, die gleichen Ausdrücke wie bei der rein 


elastischen Biegung?) und zwar: 


E— =. # ow 
2 0%’ 1 4’ 

Ba a2 080 a N ar 
ir R as ROY oy? 
_3, m ,, u 

u Bir PT" 9x 9y 


wobei w(z, y) die Durchbiegung der Platte ist. 


1) Nach einem Vortrag gehalten auf der Tagung für angewandte Mathematik und Mechanik in Darmstadt 
(April 1,950). Zu gleicher Zeit stellt die vorliegende Arbeit einen Abschnitt aus der Habilitationsschrift des 
Verfassers (Techn. Hochschule Karlsruhe, Juli 1950) dar. Die anderen Teile der Habilitationsschrift „Die 
elastisch-plastische Biegung des krummen Stabes und der Platte‘ sind wie folgt veröffentlicht: 1. Die elastisch- 
lastische Biegung des krummen Stabes, Ing.-Arch., H. 5/6 (1948). — 2. Die elastisch-plastische Biegung des 
ale Stabes unter Berücksichtigung der Materialverfestigung, Ing.-Arch., H.4 (1949). — 3. Über die 
Restspannungen bei der elastisch-plastischen Biegung des krummen Stabes, Ing.-Arch., H. 2 (1950). — 4. Die 
Berechn von stählernen Bögen unter Berücksichtigung der Tragfähigkeitsreserve im elastisch-plastischen 
N Eiteohn,, H. 10 (1950). 
2) Siehe z.B. E. Nadai: Elastische Platten. Berlin 1925, S. 19. 
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Für die Spannungen in der elastischen Zone gelten daher die Dekor Formeln: 


R Ez 2 u 
ee vB og: | . 
Ez ee 0? w — Er 
a y? TuS RE u \ (1) 
x 0: w Ez 0:w . / Y 
wobei @ = der Schubmodul und u die Poissonsche Konstante ist. Man ersieht, E 


2(1+u 
daß die en) proportional dem Abstand z von der Mittelfläche anwachsen. 
Nach der Annahme von Mises-Hencky (Annahme der unveränderlichen Gestalt- 
änderungsarbeit) muß das Fließen in denjenigen Fasern der Platte beginnen, in denen folgende 
Plastizitätsbedingung erfüllt wird: 


» e 


ty, +3 Try? =} N se EEE NE 


Hierbei ist o, die Fließgrenze beim einachsigen Spannungszustand. 
Die Spannungen in den Punkten der Grenzfläche C zwischen derelastischen und Hlastirchen 


hi. Zone der Platte (Bild 1) müssen also diese Bedingung befriedigen. 
Substituiert man in die Gl. (2) an Stelle von o,, oy und 7,, deren Werte (1) und ersetzt 
MR man z durch Z, wobei & die Ordinate der Grenzfläche zwischen der elastischen und plastischen 


Zone ist, so erhält man: 


ee nee 
el a Nager Ban ERDE Fr NETT +) 


8). 
9w \? s 
zus (5 a) 
'Daraus folgt: 
Er a RICH ec Een, 
= 2 & “)] 2 FEIN Br a re 
2) “ le ar Br 94? ap? 92 0yl j } 


worin 


u 


EEE ELLE ERENTO VE 


-N Zu4t Betn—i-n,  y=3l—u). 
Aus der letzteren I: kann (x, y) gefunden werden, wenn die Biegefläche w (z, Y) be- 
kannt ist. 
Nun wollen wir die Grunddifferentialgleichung zur Ermittlung der Biegefunktion w auf- 
- stellen. Bekanntlich lauten die Gleichgewichtsgleichungen für ein Plattenelement (Bild 3) 


. 00, 


Ep Ta 0 RR u ee. . Er a 

oM, , 2M,, 

0% Ei; °y 0% a ke OR. ” 
) u Te Vak Fr E 


Dabei bedeuten Q, und Q, die Querkräfte, M, und M, die Biegemomente, EUER das Torsions- er 
moment und p die äußere Belastung pro Flächeneinheit. Bene 


N ee AN 
Eliminiert man die Querkräfte Q, und Q, aus ar mit Hilfe von 6) und (7,0 so erhält 
0? M, +2 Hm ®M,, | 9M, er Me hi 


0%* 0% Br ‚ay? 


Wie wir schon gesehen hobnı wachsen die Span nungen O2 6 
entlang einer beliebigen Normalen zur Mittelfli 
60: Ber an. In der apa Abe een nen 


ath, 


‘a Er maßen angeschrieben werden (Bild 4): 
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+ 


il. ) Bestledigen. Deswegen kann die Definitionsgleichung für das Moment M, folgender- 


Er E 2? (02 w 02 w BEE 02 w\. (Ra TaH 
Br. ee Me Fan) de- et mer: ee) | s 
4 Eh BL a3)" (+ 02 “) | 
R u 8) N ge 
5 a 2 F 
Bi ; 
R Sa dx 
k 
E al din 
Bi 
2 Er: zZ 
2 N 
j dy y je 
Bild 3 GE 
Durch Vertauschen von x und y erhält man: N 
at Es ) Br = Er 
h mi Kar 1a 4 Earl day: + u 928 a ee Be: he (10) Par 5 A 
Für das Torsionsmoment ergibt sich dementsprechend: | A “x 
KL | | R 
E 22 ow EE 0° w /[h? a 
M;y = [ny20: an ah: Bi ug 2) Be 
1+u2 er: oy 1-+u 2%x- dy\4 ir 
ee / | . (11) IR 
# ; 4, u | ’ RE Be ke &: 02 w | ar 
Be |  TIrta\a 3) dw 0y Ben 
Man setzt die erhältenen Werte von M,, M, und M,,in die Gl. (8) ein, wobei zu beachten ist, DAR DIS 
} daß nicht nur w, sondern auch & mit « und y veränderlich ist. N 
Infolgedessen muß £ ebenfalls mit differenziert werden. Dann ergibt sich nach einer be- 0 
Fir, kannten Regel: | Li 
it” 4 3 che 
aM‘ E a -) Be 03 w pt Ur ) Bi 
a We eg en Vrrislsrrrr) Erw E 
: form SE - 2: (&) 1. # 
er Es oy? le - S 002 0x) I)’ f 
Ba 1 N: A 5% er er En) ye 2 ®) @ w 03 w e.. h? ) 3 
y2 TI Ei Oy? 0x2 3 "B% Ay oy Om: a rk N 


ar : ES aut 


cl wa 
ae ay oy 
En yy 


an a: za DB 2) 


“7 
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Nach dem Einsetzen dieser Werte in die Gl. (8) erhält man folgende Differentialgleichung: & 7 


Ma a 
. me Be ae 
= en ee 
»a-0[e-2)- 0° L oda 


en 
90y 92 aylonay m Be) F 


abi, 


Die Gleichungen (12) und (4) dienen zur Ermittlung der Biegefläche für den elastisch-plastischen 
Bereich ® (Bild 1), sowie der Grenzfläche C zwischen der elastischen und plastischen Zone. 


“ Unter Benützung von % (10) und (11) ergeben sich aus (6) und (7) folgende Ausdrücke 
4 für die Querkräfte Q, und 


° 
[4 


TE 
E ee & 
EN. ; 
u Eee 
E. ag za 
2 Für den elastischen Bereich X gilt die übliche Gleichung | “ / 
tlg tn 5 NN, as) 


a | En ES 
ı 3 : De 121) ae Er 
ist. An der Zonengrenze müssen die Bedingungen des stetigen Überganges für die Spannungen | u 


und Formänderungen erfüllt werden. 
a Bedingung (4) kann man durch die Momente ausdrücken. Man sieht aus © 00) und 
rg 


e) 0, Per ken 1.0 sah ra eek or Lie 
. . . DE ee a ae . . (IQ), 


Ei 


eh 3 a " 


a 
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das statische Moment der oberen bzw. der unteren plastischen Zone. 
Setzt man die Werte (14), (15) und (16) in Gl.(3) ein, so erhält man: 


| :-V EM MEINEM 3 M3, 222,022 208), 

x & 

} Nun wollen wir die Bedingung aufstellen, die die Momente befriedigen müssen, damit iu 
der Platte ein elastisch-plastischer Zustand auftritt. Wenn man annimmt, daß die Plastizierung 

{ an der Stelle, an der die Momente M,, M, und Mu wirken, über die ganze Plattenstärke verbreitet 


wird, dann wird in (18)& =0 sein und die Momente müssen folgende Bedingung befriedigen: 
2 
VIEH; MM, 30, =4-=W,:0. a AN) 


Hierbei. ist W, SET das plastische Widerstandsmoment je Längeneinheit. Natürlich hat die 


Formel (19) nur Danngien Wert, da die "u behandelte Theorie nur für kleine Formänderungen 
anwendbar ist. 


Setzt man in (18)& — nn 
der oberen bzw. der unteren Oberfläche der Platte: 


M}; +M}) — M. M, +3 MM}, Sa ee An eil (20), 
Y s s 


‚so erhält man die Bedingung des Fließanfanges in den Punkten 


2 
wobei W, -* das elastische Widerstandsmoment ist. 


Wir bezeichnen im folgenden den Ausdruck: 
YM: +M:2— M,M, +3M:, 


als reduziertes Moment (M,.a). Der elastisch-plastische Zustand ist also in irgendeinem Platten- 
bereich dann vorhanden, wenn in diesem Bereich folgende Bedingung erfüllt wird: 


W, O5 >= Mzea >> W, O5 . . . . . . . . (21). 


Die Bedingung M,.a < W, o, charakterisiert den elastischen Zustand, und die Re dinrun | 


M;,ca > W» 0, den nach der Plastizierung auftretenden Verfestigungszustand. 
Nun wollen wir noch die Formeln für die Spannungen aufstellen. Die Formeln (1) lassen 
sich mit Hilfe von (14), (15) und (16) folgendermaßen darstellen: 


M,% 
es) 


2 ; i 
Man führt die Bezeichnung J, = 3 &? ein, wobei J, das Trägheitsmoment der elastischen 


0 Pu 


‘ Zone ist. Dann wird 


6 er 22 
m MD Byt El 
M,2z 
ai Yy 99 
a De. 
My 2 
ey = n . (24 
| Try 4 2 SC ( ) 
Setzt man.in diesen Gleichungen 2 = +{ und benützt man (18), so ergeben sich folgende 
4 Ausdrücke für die Spannungen in der plastischen Zone: 
h Be, een, Moos 
sr MASTER "NEM M,M,+3M:,  Mya - : - (28), 
AN 
. M. A 
ee 1 BSR (26), 
P j red 
E M [07 
ne ER; > 27). 
zyl M,yea ( ) 


ni 


EUNERR n ‘ TE 


= 


2 ven TE di . LEN - a $ ( e 
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Ä | 3. Einige Sonderfälle NE, en 7-0 

“ Die Auflösung der Gleichung (12) bietet erhebliche Schwierigkeiten. Für den weiter unten ” m 
er behandelten Fall wird deshalb die Biegefunktion durch einen Ansatz angegeben, damit de 
Spannungsverhältnisse an Hand eines Zahlenbeispiels gezeigt werden könnten. 1‘ 


Ir a)Reine Biegung | F 
Es wird angenommen, daß die Biegefunktion für den elastischen Zustand eine quadratische 
Funktion 


A vw=02% +206%%Y + 0,Y° Ä Dt 


ist. Dann werden die Momente M,, M,und M,;, von den Koordinaten »und y unabhängig und 
/  Querkräfte sind nicht vorhanden. et 
ER Bei der Zunahme der Momente müssen in diesem Falle plastische Formänderungen u 
gleicher Zeit in allen Punkten der oberen sowie der unteren Oberfläche auftreten. Die bei dem 
weiteren Anwachsen der Momente entstehenden plastischen Zonen haben überall dieselbe Stärke, 
d.h. Zist von x und % unabhängig. % 
Für den elastisch-plastischen Zustand kann die Biegefläche folgendermaßen ‚dargestellt 
werden: j 


w=c()a? +2,(ö)ay +(d)y. 
Auf Grund von (9), (10) und (11) wird R 
NM, =— ya) Fuß)» 
M, =— vl) led) ra )], 
NMy,=—yO)A—u)al), 


u Tr I 


Bi worin Br 
Bi. er, | {8 
„ist. Aus diesen Gleichungen ergibt sich: | 
M ” 
ae, 
I 1—u)yl) 
1 M.—u M, 


m Se gt x? +2(1 +4) My’ ©Y +(M, — uM,) y} 5 KH 28). i 


Die halbe Höhe der elastischen Zone kann dabei aus (18) ermittelt werden. nur, 


b) Biegung nach einer Zylinderfläche Üs 


In diesem Falle wird 


ww 0° w Bi; Fe. 1:02 wi 
Dy®, ..:'.; nor 0 ya oYy: 


Dann ergibt sich unter entsprechender Verwendung der früheren. 


. (8), 
he See | (22) Er, ae iu 3 LESEN LS LAEE E TU 
2. 2 al N 1 A > N en | 2: ; 
de Ve an) le Er), 
B \ d? w iR M, M; } 
a dan EL KB a7 EM, 
IR i 
M,= 6 Os, 7 5 
3 3M, u 
| Sale are v3, 
Maa=»M, M 
f / | nu 
| Wan SAME .W,o, 21) N 
| 
| x u M, zZ 99 Y 
| Sa RER re (22) Pr 
[07 
u a 0 
Die lee (12) und a4) sind die Differentialgleichungen der Biegefläche. In der ersten 8 e; 
Gleichung wird w durch p und in der zweiten durch M, ausgedrückt. a 
4. Beispiel Ba 
Eine stählerne Platte mit A =1cm wird durch die konstanten Biegemomente u; — : RR 
580 kg/cm und M, =300 kg/em beansprucht (Bild 5). Es wird angenommen, daß M, Se 


= 400 kg/em:, E—2,1- 10° kg/em2 und u = 0,30 sd. vi av, 7 a a 
H Das reduzierte Moment ergibt sich zu A u) ? 65 
I; ; TA: 3: SE 

| Mya = YM} + M} —M, M, = 502 kgjem. | |, A 


ER ° Man ersieht aus der Bedingung (21), daß in der Platte ae 

sein N astlseheplestincher arg auftritt. Die Formel (18) My 

En lielertz = Ju «.%. 1 Er 
EEE 2 =0,350 cm. A BERE 

a Au Die Gleichung der Biegefläche (28) wird: $ 
es w = — (0,00159/em) 2° — (0,00041/cm) y3. y B: 


Sie stellt ein elliptisches Paraboloid dar. Pe E “ 


kn nr A Die ‚Spannungen i in der elastischen Zone ergeben sich nach (22) und (23) zu: 
Sc E: 
em?) ABER NH 


BE = 2 2770 kg/cm?; | Oyp = MR 1430 kgfem. 


, ER die Spannung Op die Fließgrenze o, wesentlich überschreitet. 


= 

u 

a 

203 

“ 

“ £ 
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Über einige unendliche Reihen von Bessel- 
schen Funktionen. 

Essollen Summenformeln für einige unendliche Reihen 
von Besselschen Funktionen mitgeteilt werden, auf 
die bisher in der Literatur nicht aufmerksam gemacht 
worden zu sein scheint. Neben den Besselschen Funk- 
tionen Jn (2) mit ganzzahligem Index n treten in ihnen 
die sog. modifizierten Besselschen Funktionen y(2) 
auf, die folgendermaßen definiert sind 


In) = iR In(iz) ). 


Die Funktionen J„ lassen sich bekanntlich durch eine 
erzeugende Funktion darstellen 


s Ve N oe 


Sun 
[A 


— irn (e) 27. 


(Alle Summen, bei denen nicht ausdrücklich Grenzen 
‘angegeben sind, sollen hier stets von — © bis + © 
erstreckt werden.) Wegen (1) gewinnt man unter Er- 
setzung von 2 durch iz , und von it durch 7 für die /„ 
die Darstellung 
Ztr+r) 


rn. tal. 
N 


Multipliziert: man (2) mit (3), setzt 7 = i”1 und ordnet 
das Produkt der Reihen in geeigneter Weise, so erhält 


man 
ei _ 2 '2 Inka) Inzrei.t  {aB; 
(Auf die ae der Umformung soll hier nicht 


eingegangen werden.) Entwickelt man die Exponen- 
tialfunktion in eine Reihe nach Potenzen von t 


2 


und vergleicht: (5) mit: (4), so gewinnt man folgende 
Beziehung 


. 6) 


DIE Te ae (6), 
N 
; 
= In) In-ı@) füri=1,2,3,.....(M), 
& n 
0 = YIn(z) Inı(e) = I(—1)" In(2) In-ı 
n füri=1,2,3,. Ko: 


Das zweite Gleichheitszeichen in (7’’) folgt daraus, daß 
In) = (— DV" In); In) = Ink) . . (8). 


Durch Addieren und Subtrahieren von (7°) und (7) 
erhält man schließlich 


12° 
31 =, ni m(?) Ey 2) = 377 rm4ı®) T5 m-+1-1(%) 
für DL E00). 


Es ist wegen (8) zweckmäßig, nur Indizes >0 en 
‚zu lassen. Man gewinnt so die Summenformeln (unter 


Fortlassung der Argumente 2): -_ 

Ih 2 Fonda -1 (oa), 
Jul; + ZIom (len 1 + Iamtı) > (10b'), 
Jlo+ B> Ja m+ı 7 Jam-1) Izm = = (10b’’), 
It Z Famllameet I, m+2) -£ (10e°), 


FR 
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Jıl, Et: > (Jam | Io m-ı 
L 


5 Jon Tom+ı) = =: (10e’’), 
Jols Ei J5 (I, Tr 15) 
oo 35 N ’ 
z 22 J5 m+2 (1; m-1 % I; m +5) — ip) (104), 
1 


J(— 1,) Er > Jam ‚ı Jam 
} 


15 m+4) 12 (104”) 


USW. 


Für den allgemeinen Fall lassen sich nicht einfache, 
übersichtliche Formeln gewinnen, die nur Indizes =0 
enthalten. Man muß derartige Formeln für jedes 
einzelne / aus (9) ableiten. 

Die angegebenen Formeln lassen sich auch durch 
einen etwas mühsamen Grenzübergang aus dem Addi- 
tionstheorem für Zylinderfunktionen gewinnen. Jedoch 
scheint uns die hier gegebene Ableitung angemessener 
zu sein. 

Die Gl. (10) treten bei der Untersuchung von großen 
Amplituden in Magnetfeld- und Traveling- Wave- 
Röhren!) auf. 

Paris, Centre des Recherches de la Compagnie gen6- 
rale de T. S. F. 

Hamburg, Institut für Theoretische Physik der Uni- 


versität. 
[07 Döhler und G. Lüders., 


!) J. Brossart und O. Döhler: Annales de Radiodlectr,, im 
Erscheinen, 


Zur Frage, ob alle Wurzeln einer algebraischen 
Gleichung einen negativen Realteil haben 
(Stabilitätsfrage). 

Die verschiedenen algebraischen und geometrischen 
Stabilitätskriterien setzen insbesondere algebraische 
und funktionentheoretische Kenntnisse voraus, die 
nicht jedem Ingenieur geläufig sind. Überdies ist der 
Zusammenhang des Schemas der Anwendung dieser 
Kriterien mit ihrer Begründung nicht so einfach, daß 
es möglich wäre, diesen Zusammenhang während ihrer 
Anwendung gegenwärtig zu haben. Bei dem folgenden 
einfach zu begründenden Stabilitätskriterium ist der 
Beitrag, den jeder Schritt der Rechnung zur Be. 
antwortung der Stabilitätsfrage leistet, besonders 
leicht zu vergegenwärtigen. 

Der enge Zusammenhang des Kriteriums mit den 
Kriterien von J. Schur [l]ist am Ende angegeben, 

Wir beschränken uns auf den technisch wichtigsten 
Fall positiver reeller Koeffizienten und führen in die 


mit 
ım)= ntrapteptt: 
hp)=aPptep’tap+t' 
gebildete Polynomgleichung EN 
ı :  AdSssWtHhp)=0 .x 2) 


den positiven, reellen Parameter 2, ein, Ben wir h( ») 
mit A, multiplizieren. 


Wir erhalten 
4,0) =9(P) + A hip) = 
und behaupten: i 
| Hat Alp) : =0 
N, Wurzeln mit positivem Realteil (R >0) 
N, Wurzeln mit negativem Realteil (R<0) 
N, Wurzeln auf der imaginären Achse (R = 
so hat auch A FRrDE ebensoviele Wurzeln in jeder 
beiden Halbebenen und auf der Bi a As 3 


und 


{ we 
£ , 
& Er ra a * N 
Fa BR iz .- ak 
A Ins 2 il Pr‘ we m PR 
Ir Ve 
KArıh R & aa! BR nah, 
aR " se. 
"N et = 
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j Durch die Einführung von A, kann keine Wurzel die 


imaginäre Achse verlassen oder überschreiten oder auf 
die imaginäre Achse rücken. Dies folgt am ein- 
| fachsten aus der Tatsache, daß die Lage der Schnitt- 
punkte der Wurzelkurven a 
| P (a, y) 
9 (x, y) 2 ee) 
von "A(p) = P(&, y) + iQ(&,y) = 0, = CH iy 
mit der imaginären Achse von /,> 0 unabhängig ist. 
Diese Scehnittpunkte y ergeben sich als Wurzeln der 


Gleichungen (3) für © = 0; diese Gl. lauten: 
P,Yy=% — Gy ray — 9 =0 
N) - an - array —=0 .. (A). 


EN 
2 DW HM DW 08:05:04 02 0 


ıq 
\ 


10:08 -056 04 02 0 


—r 
70 -08 -06 -04 -02 0 


„20 18-16 14 12 


10 -08 -06_-04 -02 0 


2 

LEE EN ERS Q 

? 
10 
os 
s 
Fr’ Ye 

N ER EEE en ER RE = 

4ıp) ArziP) 
’ Teil b 


Die Gleichungen der Wurzelkurven für A 4,(p) unter 
scheiden sich nur dadurch von denen für A(p), daß 
Q,(y) mit A, multipliziert ist. Die Wurzelkurven von 
A(p) und A,(p) schneiden also die imaginäre Achse 
indenselben Punkten, wie dies in der Darstellung 
- der Wurzelkurven zum 1. Beispiel mehrfach an- 


und (2) denselben Wert. 
Da ferner die Wurzeln von (1) stetige Funktionen 


_ Realteil einer Wurzel von A(p) infolge der Einführung 
von /, sein Vorzeichen wechseln. Könnte nämlich eine 
Wurzel die imaginäre Achse überschreiten, so müßte 
es auch einen Wert A, geben, für den R=0 würde. 

Das ist aber nieht möglich, da A, die Schnittpunkte 


A; 


gedeutet ist. N, hat also für die Gleichungen (1) 


der Gleichungskoeffizienten sind, so kann auch kein 


» 
der Wurzelkurven mit der imaginären Achse un- 
verändert läßt. 
Für die Beantwortung der Stabilitätsfrage kann 
also die Gleichung (1) durch (2) ersetzt werden. Setzt 
man nun in (2) 


hı= wm) __atapft capı Bir=" 
h(pı) G1Pıtesp+*-- 
worin 9, negativ, reell ist, z. B. 
9,=—]1, 
so wird 


A) = ap) hip)— g(pı) A(p) =O 


use 
30:08 -06 :04 02 0 70 08 -06 4 0 0 


30 


20 


ost 
0% Es existiert noch eine 
Ü zweite reelle Wurzel 
02 beix--346 
(nicht dargestellt) 
Re nr Er ET 
d;(p) 2r,(P) 


Teil c 


10 -08 -06 -04 02 0 10 08 06 04 -02 0 10 -08 -06 04 02 0_, x 
— 


q 
30 N SL 
52 ? 
20 
Q u 
10 
az 
067 Es existiert noch dieselbe 
reelle Wurzel 
ya beix.-3.46 p 
wie bei Ans 
! ge $ (nicht dargestellt) io 
en a I er De = 
* Z,(p) 7,(P) 4,1P} % 
Teil d uud e 


ZuBeispieli1. Wurzelkurven der zu: 
A(p)=8p°+19 9° + 1209! -+137 9° +450 9° + 2509 +322— 0 
gehörigen reduzierten Gleichungen 4, (p) =0 und der Gleiehun- 
gen d,(p)=0. 
4 S$ = Wurzelschwerpunkt 


Dieser Ausdruck verschwindet für p — p,; er ist also 
durch p — p, teilbar, so daß wir die Gleichung 


_ 20) Ko) — oo) hp) | 
ii | . 
=nPp)+h(WM)=0 « 


vom Graden— 1 erhalten, die in der negativen Halb- 
ebene eine Wurzel weniger hat als die Ausgangs- 
gleichung; ihr haben wir durch die Wahl des Para- 


A,(p) 


(6) 


*- 
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meters A, zunächst die Wurzel p, aufgezwungen, die 
wir sodann abgespalten haben, 

Weist Gl. (6) einen Vorzeichenwechsel auf, so ist 
sicher N, 720, die Stabilitätsfrage also zu verneinen; 
hat (6) keinen Vorzeichenwechsel, so muß man mit (6) 
wie mit (1) verfahren, wobei man aber auch eine be- 
liebige andere (p} > — 1) negativ reelle Wurzel p, ab- 
np: alten kann. Man erhält so eine Gleichung vom 


Grade n— 2 mit den Anzahlen N,, Nn_s, N, der 
Wurzeln. Man hat also zu bilden: 
AP =nP)+Rhılp) ....(M 
y 91(Pı) 
mit Mo = — —— 
s h(Pı) 
und aus (7) ‘ 
A, ,(P) 
4.p = —— =: +h,lpP)=0. (8). 
P=Pı 


Die Vorzeichen von A,(p) entscheiden über die Not- 
wendigkeit der Bildung von A,(p). 

Dieses Verfahren der Reduktion des Grades der vor- 
gelegten Gleichung setzt man fort, solange in der 
reduzierten Gleichung ER Vorzeichenwechsel auf- 
tritt, längstens bis zur Gl. Grades, die man nach 
n — 2 Schritten erreicht. Sina alle reduzierten Glei- 
chungen von Zeichenwechseln frei, so ist N,— 0; da 
N, durch das Verfahren nicht verändert wird, so 
findet dann die Folge der Gleichungen A,(p) mit 
derjenigen Gleichung ihr Ende, welche nur noch die 
rein imaginären Wurzeln von A(p) hat, falls A(p) solche 
Wurzeln besitzt. Im allgemeinen ist es zweckmäßig, 
für alle Reduktionsschnitte p, = — 1zu setzen, so daß 

_ ot gt + ES} 
RE 
wird; die notwendigen Divisionen durch (p-+]1) 
führt man mit dem Hornerschen Schema aus. 

Schur [1] ordnet dem vorgelegten Polynom A(p) 
vom Grade n ein zweites /*(p) von demselben Grade zu, 
dessen Wurzeln in bezug auf die imaginäre Achse 
Spiegelbilder der Wurzeln von A(p) sind; er benutzt 
den vom Hermite [2] bereits bei ähnlichen Unter- 
suchungen verwendeten Ausdruck 


F(p, pı) = 4*pı) AP) - Alp) Ap)=0 ... 9) 
für das Reduktionsverfahren; F(p, p,) ist durch p — pı 
teilbar. 
Setzt man bei reellen Koeffizienten 
Alp) =og(p) + Ph(p) 
A*p) =a,g(p) + Pıh(p), 


so wird 
0 = |: le@ sn swo ro) 
oder 
ER ap, 
pP, P1) = h(pı) aß w-% un | 


Der Ausdruck (9) stimmt also bis auf einen ‚belang- 
losen Faktor mit (2) bzw. mit (5) überein, 


1. Beispiel;stabiler Fall, 
Vorgelegt sei die Gleichung 


Alp) = 8p° + 1395 + 120p* + 137p° + 450p° + 


+ 250p + 322 = g(p) + h(p) = 0. 
Mit p, =— 1 wird 
a 322 +450+120+8 _ 2,25 
250 + 137 +13 


und 
A,(p) = 8p° + 29,25 p5 + 120 pi + 308,25, p° + 

-+ 450 p? + 562,5p + 322 = 0. 
Wie auch das Bild der Wurzelkurven zeigt, hat A, (p) 
— und ebenfalls A, (p), 4,,(p) und A, (pP) — dieer- 
zwungene Wurzel p, = — 1; ein ursprünglich kom- 
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plexes Wurzelpaar ist durch die Kinführung y von en 
in zwei reelle Wurzeln zerfallen, deren eine (p,) mit 
Hilfe des Hornerschen Schemas abgespalten wird, Rev 


Das Schema liefert: 


8 29,25 120 308,25 450 562,5 | 322 
—1] —8 -—21,25 —98,75 — 209,5 — 240,5 — 322 
8 21,25 98,75 209,5 240,5 322 


Die 1. reduzierte Gleichung ist also 
A,(p) =8p° + 21,25 9? + 98,75 p® + 209,5 9° + 
+ 240,5p + 323 = 9,(p) * hp) = 0. 
Mit: A, = 1,59 ; 
wird Wer (p) = 5896 p? + 9838,75 pt + 72778 ‚75 p° En 
- + 96998,65p? + 177248,5p + 149086 = 0 
und 


As(p) = 5896 p! + 3942,75, p° + 68836,p2 + 28162,5 p 
+ 149086 — 9,(p) + h,(p) = 0. 


Mit 7, = 6,97 findet: man in derselben Weise 


A,(p) = p° + 3,662 p? + 8,013 p + 25,286 — 9, (p) 
+ hp) = 0 


14 = 3,212 
4,(p) = p?+ 0,1409 + 7,873 = g,(p) + h,(p) = 0. 


Da in der Reihe der A,,(p) keine der Gleichungen einen 
Vorzeichenwechsel aufweist, hat: A(p) keine Wurzeln 


und:mit 


-mit: R>0; da ferner A N) keine rein imaginären 


Wurzeln hat, liegt auch der Grenzfall der Stabilität :g 
nicht vor. A(p) hat also nur Wurzeln mit R<0. j 


2.Beispiel; GrenzfallderStabilität. 
Ap=p’+4p+sp +l2p + 13pP +2 
+6p+4=0 
führt mit b r 


9 


8 7 
2ı re, hen joe S 


‘zu den reduzierten Gleichungen 
A,(p) =8p° + 2 + 44p! + 40p? +'64p? + 20p 
3 


=,P)+hp)=0 AIR 
a u Be ar ee + 40p + 140 f 
= 9:(p) + he(p) = b 


4d,(p) = 28 pt + 56? +28 = Sn +h(p)=0 


mit dem en De —+i,dasauch Wurdel« 
paar- von A(p) ist. 


3. Beispiel; instabiler Fall. 

Das von Schur [1] gegebene Beispiel ist 
Ap=8pP+4p+7p+5®R Hp +1=0. 
Wir erhalten mit 2, -5 Baader: "z 


4; (p) = 40p° + 32p' + 35p° + or +5p+8= 

Das Hornersche Schema 

40, 32 35 40 bi 38 
—40 +8 —43 +3 

4 —8 +43 — 3 Bu; 


liefert in der letzten Zeile die Koeffizienten ° u Alp) 
mit Zeichenwechiel; es liegt also petaiBEe von 
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Das Mehrstellenverfah ren bei Plattenaufgaben. 
“ Beiden Randwertaufgaben der Potentialtheorie hat 


sich unter den verschiedenen Arten von Differenzen- 


. verfahren das Mehrstellenverfahren als brauchbar ge- 


zeigt. Es liegt daher nahe, es auch bei der Platten- 
gleichung anzuwenden. Das Mehrstellenverfahren ist 
eine der verschiedenen Verschärfungen des gewöhn- 
lichen Differenzenverfahrens, welche den Vorzug hat, 
bei Differentialgleichungen einfacher Bauart die Re- 
chenarbeit gegenüber dem gewöhnlichen Differenzen- 
verfahren nur wenig zu erhöhen, während die erhalte- 
nen Resultate Fehler aufweisen, die mit einer höheren 
Potenz der Maschenweite h gegen Null gehen und da- 
her im. allgemeinen wesentlich kleiner sind, als beim 
gewöhnlichen Differenzenverfahren. 


1. Formeln des. Mehrstellenverfah- 
rens im quadratischen Gitter für 
die Plattengleichung (2) 

Hier soll nicht das Mehrstellenverfahren allgemein 
entwickelt !), sondern nur seine Anwendung auf die 


"Plattengleichung beschrieben werden, Wie bei den 


anderen Arten von Differenzenverfahren wird ein 
etwa quadratisches Gitter von der Maschenweite A 
mit der Gitterpunkten ach Bild1 

n=mrih, „eytkkh...() 
verwendet. Zur Anwendung N die Plattengleichung 
u Our, du 9 S 


AAu— na dy! Saaayaı N. . (2) 


Bild ı 


(mit den üblichen Bezeiehnungen: Durchbiegung 
u (2, y) der Mittelebene der als dünn vorausgesetzten 
homogenen Platte an der Stelle x, y, Belastungsdichte 
p (x, y) an der Stelle x, y, konstante Plattensteifigkeit, 
N) wird eine Summe gebildet _ 


Pal Jh. Gr A (AA) X an (3) & 


wo die nn über eine gewisse Anzahl von Gitter- 


_ punkten zu erstrecken ist, bei u und AA u die ange-_ 


 hängten Indizes j, k die Werte an der Stelle %,, Y, be 
zeichnen und die Konstanten @,4, A;,s0 zu bestimmen 


Pe sind, daß bei Taylorentwicklung aller Glieder der 
Summe nach oo und den partiellen Ableitungen an 


& &+ ; 
telle 20 Yo die Koeffizienten von a) R j 
0,0: 


0x% Oyf 
ae hohen Be von & und ß ver- 


Er Nr war 


RER des Mehrntellekerlährens BA: 


468 u, — 144 3 up — us 


BIT EP 
22 Zn SE En Mo = ° & "5 


enaeben mit a Anw en 
if., 361 ff. uf 2 


eine Gleichung für die Näherungswerte U, (Näherung 
für u,,) sein. (In dieser Gleichung wird die Differen- 


. tialgleichung (3) an ‚‚mehreren Stellen‘“ %;, Y, Ver- 


wendet, während sie bei (6) und (8) nur an der Stelle 0 
benutzt wird.) 

Bei der Plattengleichung (2) wird man unter Aus- 
nutzung der Symmetrien als Ansatz (3) etwa ver- 
wenden 


B=a,yUugt ap üp +08 Ug + An | 
En a8, us + Ay AAuy 
a... | 

wobei Up die Summe der u-Werte in den vier Punk- 

ten P, also nach Bild 1 


Düp —UyotUu,ıt %-,0ot+ &u-ı 
bedeutet, ebenso sind die anderen Bezeichnungen 

RE 3 AAup, ... wohl ohne weiteres ver- 
ständlich. 

Bei der Durchführung des eben beschriebenen Ab- 
gleichs der Koeffizienten in der Taylorentwicklung 
erhält man als einfachste Möglichkeit den Fall des 
gewöhnlichen Differenzenverfahrens 


0ug—- 8yup +23 ug +3un — htAA ol 
— Restglied 6. Ordnung] 


und als genauere Formel 


(5), 


36 U, — 10. 3) up — ER ug + Up + PAR 


EN 5 ae HT): 


— Restglied 8. Ordnung 


Gegenüber anderen Arten der Verbesserung des 
Verfahrens, wie z. B. der Formel 


184 u, — 77 3 up + 20,3 ug + 14 3) un 
mans SD üm— 6ht AAu, . (8) 
— Restglied 8. Ordnung 

hat das Mehrstellenverfahren den Vorzug, bei gleichem 

Annäherungsgrade nicht so weit herausreichende 

Punkte zu verwenden, was besonders in Randnähe 
wichtig ist. 

Man kann sogar, ohne das Quadrat von der Seiten- 

länge 4h zu verlassen, auch einen Ausdruck des Mehr- 


stellenverfahrens aufstellen, bei dem das Restglied von 
10. Ordnung wird:, 


18 un+8 dus 
Eu 332 AAu,g+ 72 9 AAup 


+ 26 3/44 ug — D!AA un) 
— Restglied 10, Ordnung. 
Dabei muß man allerdings unangenehme Koeffi- 
zienten mit in Kauf nehmen. 


2; Beispiele mit quadratischem 
Gitter 


"Einige ganz einfache-Beispiele mögen die Anwen- 
dung der Formeln erläutern. Es wurde darauf ver- 


zichtet, durch Verwendung sehr feiner Gitter eine hohe 
Genauigkeit. zu erzielen; für technische Zwecke wird 
es oft genügen, mit grober Maschenweite und geringer 


 Rechenarbeit einen Überblick über die Lösungsfunk- 


tion. zu erhalten. 


% Beispiel I: Eigenschwingungen einer quadratischen 


_ homogenen ringsum frei aufliegenden Platte von der 


x  Seitenlänge A. 


Sau der Differentialgleichung 
k% RANDE SHN. u a) 
ER 28 


= UP 
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% Baer » 
treten dann die Randbedingungen 
u=Au=0 


am Rande des Quadrates, 


) ist dem Quadrat der Eigenfrequenz & proportio- 
nal. Die Eigenwerte sind hier exakt bekannt, so daß 
man für die Näherungswerte die Fehler angeben kann. 
Die Tafel1l gibt die für verschiedene Maschenweiten 


nach (6) und nach (7) erhaltenen Näherungswerte AA! 


für die Schwingungsformen, die bezüglich der Diago- 
nalen D,, D, und Seitenmittellinien S,, S;symmetrisch 
sind (Gitter nach Bild 2). 


Y 


Tafell 


Gewöhnl. Difi.- Verfahren Mehrstellen- 


verfahren Gl. (7) 


n=A | AA = 324 (-17%) | 405 (+ 4%) 


exakter erster Eigenwert A, mit 
}, A! = Ant 389,64 


Beispiel IL: Eine homogene Platte mit dem Sechs- 
eck BCDEFGB als Grundriß (Bild 3) sei längs BCD 
eingespannt, längs der anderen Sechseckseiten frei 


| Bild 3 RR 


aufgelagert. und gleichmäß' g belastet“ Ni const- 


Nach (7) hat man bei der Maschenweite = 4/2 die 


drei Gleichungen 


ige m rmen 


mit der Lösung 


N i Bild 2 


‘erhält dabei aus 


, 351,4(-10%)|394,01 (+ 1,1%) 
h=7 |44:= 14096 10 240 
I 937 37258 


ag. eier. 
"= 9940700180 At, 2 
616 
0,0 24 
Ne ze 
425 v 
ge le PER: 


Die Rechnung nach dem hen Differenzen- 
verfahren mit Gleichung (6) re ebensovie] 


Rechenarbeit bei weniger guten Ergebnissen. Man 


18a — 1b + 2:=0= 0 
— 84+21b—- 8c=e % 
a— 8b TEN 3 
die Werte 
709. „P 
= au e = 000 A02 £ 
612° p 
ER 4 & 
381205 9,0100 A » > 
/ i “ dal f ! 
ee 38207 0,0069 A E- PIDER % 


"Zur or chkun von Spannungen und Momenten 
kann man die Differenzenausdrücke wie beim gewöhn- 
lichen D.fferenzenverfahren, aber gebildet mit den 
nach dem hier beschriebenen genaueren ee be- 
rechneten Funktionswerten, benutzen. 


\ 


ir 3. Dreiecksnetz r 
Da Dreiecksgitter (der Einfachheit halber oh nur 
Gitter mit gleichseitigen Dreiecken betrachtet werden) 
in der Literatur masgterholt aufgetreten : sind le Ei Be 


ser 


BEE. © 
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x 
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nen wir uns auch hier kurz fassen. Es werden wie in 
° Nr. 1 Ausdrücke der Form (3) gebildet und unter Aus- 
nutzung der Symmetrien, Bild 4, der Ausdruck B in 
der Form (5) angesetzt, wobei die Summen wieder 
über die symmetrisch gelegenen Punkte zu erstrecken 
sind, also z. B. Su, geht über die sechs Punkte 


Sup=ulag + h,yo) + ul — h, Yo) 
h Ba 
+u (" SEIT 23) 


h h 
ru (= ER Yo 3 /3) 


h 3) 


2 


RG Re 
+u (= 3.9 373) 


Wieder sind alle Glieder in (5) nach dem Taylorschen 
Satz an der Stelle x,, y, zu entwickeln und die noch 
zur Verfügung stehenden Koeffizienten Ag Apı ..: 


sind so zu ermitteln, daß in der Taylorentwicklung das 


h 
+u ( + ER ze 
En 


geben sich so die Formeln ®) 
. Int 
2u9 3 up+ Zug 75 AA 
| = Restglied 6. Ordnung 
42u, — 10 3 up +23 ug + PX 
3 

BT n4(6 AA ug + >; AA up)g (U). 

j = Restglied 8. Ordnung 


(10) 


4.Andere Gleichungen 


Für andere Gleichungen hat man die Formeln des 
Mehrstellenverfahrens jeweils durch Taylorentwick- 
‚lung neu aufzustellen ; insbesondere, wenn die Platten- 


steifigkeit und damit die Koeffizienten der Ableitungen _ 


von win der Differentialgleichung Veränderliche sind, 
kann die Durchführung des Mehrstellenverfahrens 
sehr mühsam werden. Wenn die Koeffizienten der Ab- 
leitungen von win der Differentialgleichung aber noch 
konstant sind, lohnt es sich, durch einmalige Aufstel- 
lung der Taylorentwicklung sich die zugehörigen For- 
meln des Mehrstellenverfahrens herzuleiten, eine Ar- 
"beit, die dann fürandere Randbedingungen und andere 
Belastungen: übernommen werden kann. 


Restglied von möglichst hoher Ordnung wird. Es er- 


Als Beispiel greifen wir die Differentialgleichung 
heraus 


e DT 
Ju Beh u +J Ser 
MER I 7 a uU TEE 
x Y b 
(kreuzweise bewehrte Kisenbetonplatte; einzelne 


Kisenstäbe bilden unter Mitwirkung des Betons eine 
fachwerkartige Scheibe. Näheres bei Mareus, a.a. 0. 
S. 105). 

Für J_ = J, entsteht 


U 


D 
u ae de 
zzyy yuuy 


EB, J 


% 


Liu] Su ggen 

Wieder wird unter Benutzung der Symmetrien ein 
Ansatz B nach (5) gemacht, wobei nur AAw durch 
L[u] zu ersetzen ist. Man erhält als einfachste Möglich- 
keit die Formeln des gewöhnlichen Differenzenver- 
fahrens 


16u,— 6 NY up + PT, + up -- hL[u]) 


(12) 
—= Restglied 6. Ordnung J 
und als genauere Formel 
40 u, — 10.3 Up— 4 N ug + 2 Sup 4 PAR 
2 ht 
Fans 2 Lu, +3 L[u]p} (13). 


—= Restglied 8. Ordnung 


x AAUu=%k | Urer + Urzryy + uyyyy =—E 
Gewöhnl. Diff.-Verfahren Mehrstellenverfahren Gewöhnl. Diff.-Verfahren Mehrstellenverfahren 
Gl. (6) Gl. (7) Gl. (12) Gl. (13) 
. a 1971 3 - 5407 | ‘1221 309593 
= ne nn 23 ee ——— I 265 
E kA% 625 - 1522 D:0p 207 625 11398 0,00228 | 625 - 802 Re 625187172 No 
. & ee = u z phmanähien 
f b 1053 32893 636 4 160965 
; er a0 u 192. 137 
E e kA? 625 - 1522 2001 62511398 2 625 - 802 u 625187172 ODER 
| 7: 13874 Aul% 3.3587 795 201699 
-. = a m a nn a ——dy 2 
“ kA* 625 - 1522 Kamin: 62511398 0127 625 - 802 O:0222 625 -187172 wur 
2 |. 4698. 28= „ | 3.1942 421 107748 
r — u —— {DU — 092 
.kA* 2 625:1522 leg '625- 11398 a? 625 - 802 ame 625: 187172 
Hin \ h 


*, Bei Göttlicher,a.a. O., is. die Formel aufgestellt (um- 
Ry ges: 'hrieben auf die hier benutzten Bezeichnungen: 


u Fi Le ” ; 1 
Beau Op +2 ZugtZun-4arg=0 un) 


Als einfaches Beispiel sei eine quadratische, homo- 
gene, gleichmäßig belastete Platte der Seitenlänge A 
gewählt, die an zwei gegenüberliegenden Seiten einge- 
‚spannt, an den anderen beiden Seiten frei aufgelagert 
‚ist. Bei der Maschenweite A = 4/5 und vier unbe- 
kannten Funktionswerten a, b, ec, d nach Bild 5 wurde 
für die beiden Gleichungen 


Adv =% (14) 

und 
n N ER.) 5 
} RBIRR ıv De, = Ki. (18) 
28* 
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sowohl nach dem gewöhnlichen Differenzenverfahren 
(6) (12), als auch nach dem hier dieselbe Rechenarbeit 
verursachenden Mehrstellenverfahren (7) (13) gerech- 
net. Die Ergebnisse der Tabelle zeigen, daß das ge- 
wöhnliche Differenzenverfahren Werte liefert, dienoch 
um etwa 10% unterhalb der vertrauenswürdigeren 
Werte des Mehrstellenverfahrens liegen. 


Hannover. L. Collatz. 


Über die Koppelkurve des Schubkurbelge- 
triebes. 

Es sollen einige für die Getriebesynthese wichtige 
Eigenschaften der Koppelkurve eines allgemeinen 
Schubkurbelgetriebes abgeleitet werden, die nicht ein- 
fach durch Grenzübergang aus bekannten analogen 
Sätzen für das endliche Dreistabgetriebe gewonnen 
werden können, weil die Koppelkurve des Schubkur- 
belgetriebes auch von den imaginären Kreispunkten 
verschiedene uneigentliche Punkte besitzt. Zur be- 
kannten Gleichung der Koppelkurve wird zunächst 
eine Bemerkung über die Orientierung des Koppel- 
dreiecks und des Koordinatensystems hinzugefügt. 
Hierauf werden die Asymptoten der Koppelkurve be- 
stimmt und es wird nachgewiesen, daß die reellen 
Sehnittpunkte der beiden Paare konjugiert imaginärer 
Asymptoten auf der Verbindungsgeraden der beiden 
Doppelpunkte der Kurve liegen. Dabei besitzen die 
Paare der Doppelpunkte und der Asymptoten-Sehnitt- 
punkte denselben Mittelpunkt. Zu einer Koppelkurve 
eines Schubkurbelgetriebes gibt es stets noch eine 
2. Erzeugung dieser Art. Die Gerade, die den Schnitt- 
punkt der Kreuzkopfbahnen dieser beiden Getriebe 
mit dem ihnen gemeinsamen festen Kurbeldrehpunkt 
verbindet, steht senkrecht zur Verbindungsgeraden der 
beiden Doppelpunkte der Kurve. 


1.Die Gleiehung der Koppelkurve. 
F,ABF,. sei (Bild 1) ein allgemeines Schubkurbel- 
getriebe mit der Exzentrizität e, dem Arm FR, A=a 
und der Koppel AB=c. Der Punkt X der Koppel- 
ebene sei durch das Koppeldreieck AXB mit der Kop- 


Bild 1 


pel fest verbunden und so gelegen, daß das Dreieck 
AXBin diesem Sinne positiv umlaufenist. Wir führen 
ein kartesisches Rechts-Koordinatensystem mit dem 
Ursprung F, und der x-Achse F, F,„ ein, wobei die 
Richtung von F, gegen F,. bei e +0 mit der gegen 
(B) übereinstimmen soll. Ist dann AX =», BX = 
und sind die Winkel des Koppeldreiecks bei A, B, X 
gleich &, ß, y, so lautet die Gleichung der von X be- 
scehriebenen Koppelkurve!) 


U? + V2— W? 
mit 
U=ucoy(a + y?+ v2 — a) — r(z — e), 
V=usiny(e+ y*+ v%— a?) — 2uy (2 — e), 
W = 2uv (2sin y — y008 y). 
*) Vgl. F. Ebner, Leitfaden d, techn. wichtigen Kurven 


, Leipzig 1906, 8.18 (Gleichung der Koppelkurve für ein negativ 
id, umlaufenes Koppeldreieck AXB nd) ein rien): AN 


Die Kurve (X) gehört also zu er Arkılaren Kurven 
4. Ordnung. 

Bezüglich der obigen Orientierungsfestsetzung | sei 
bemerkt: Bei jedem Schubkurbelgetriebe ist zu jeder 
gegebenen Stellung F, ABF, „ (Bild 2) auch die be- 
züglich der Steggeraden F,F,% symmetrische Stel- 


jung F, ABF,„ möglich. Es sei X ein Punkt der 


Koppelebene, der so gelegen sei, daß das Koppeldrei- 
eck AX Bin diesem Sinn positiv umlaufen ist. Ist X*, g 


X” der zu X bezüglich der Koppel; bzw. Steggeraden 


symmetrische Punkt, so ist X* die zu der Stellung AB 
der Koppel gehörige Lage des Punktes X*. Die Bahn- 
kurven der Punkte X und X* sind also bezüglich der 4 
Steggeraden symmetrische Kurven. Beziehen wir 
demnach die Koppelkurve (X) auf das Rechtssystem \ 
F,(#, y), die Koppelkurve (X*) auf das hierzu bezüg- - 
Ber der Steggeraden symmetrische Linkssystem 
F,(x, Y), so stimmen die Gleichungen der beiden Kop- 
pelkurven miteinander überein. Wir legen also unseren 
Überlegungen ein im + Sinn umlaufenes Koppaldrei- 
eck AXB und das Rechtssystem F,(z, y) zugrunde und 
können alle Formeln auch für ein — umlaufenes 
Koppeldreieck AXB verwenden, wenn wir sie in die- 
sem Falle in dem Linkssystem r, (2, y) deuten, 


| 
e 


2. Die Doppelpunkte der Koppel- 
kurve U=0, V=0 sind Kegelschnitte, W = 0 
eine Gerade g, die gegen die +x-Achse unter dem 
Winßely geneigt ist. Gemeinsame Schnittpunkte 
dieser drei Linien sind Doppelpunkte der Ko duo \ 
(X). Die Gerade g schneidet den Kegelschnitt U — 0 k 
in 2 Punkten mit den Abszissen E 


a 2 2 + } Vale 7 ; 3 
RT v2 + (a Ver tee) ‘v) SR ‘ j 
OP a 
. Diese Schnittpunkte RN auch auf dem Kegelsehitt 
V = 0 und sind daher Doppelpunkte D,, D, der Kop- 
pelkurve. Die auf der Geraden g En rel Eure 
nungen F\, D,, F},D, sind (Bild Chr j PERZATNE 


ya = 


ee MEN LIESCSNEN. Be: 
A ad - - = ae } 


aan aon | 


Die allgemeine Koppelkurve a: Schuber 
bes ist demnach eine zirkulare Kurve 4. Orc 
ae en auf der Geraden g durch F,, 
F, F, „ unter dem Winkel y geneigt ist 2. Die: 
schaften genügen jedoch noch Ei Ke 
nung dieser Koppelkurven unter Ku i 


nung. je Di ng 
3.Die Apymptoten 
a) Die . 


man die Gliedergruppe 4. 
der ope)kur sp gleis ) N erg 
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Seite zeigt, daß die Koppelkurve einfach zirkular ist. 
Die restliche quadratische Gleichung lautet 


u— Auvceosy+t 4: 


= 0 


u: vsin 
m?’—4 en Ah 


und liefert die beiden konjugiert imaginären Wurzeln 
2vsinyHi(u — 20c08Yy) _w 
u Ira, 

Im Gegensatz zu der Koppelkurve des endlichen Drei- 
stabgetriebes besitzt also die des Schubkurbelgetriebes 
unendlich ferne Punkte auch außerhalb derimaginären 
Kreispunkte. 

b)Die Asymptoten-Schnittpunkte., 
Die Asymptoten, die in den imaginären Kreispunkten 
berühren, schneiden sich im festen Drehpunkt F, der 
Kurbel; denn setzt man in der Gleichung der Koppel- 
kurve y= -+ia, 2°-+ y®=0, so verschwinden die 
Koeffizienten von x* und 2°. F, ist demnach singu- 
lärer Brennpunkt der Koppelkurve und er liegt auf der 
Verbindungsgeraden der Doppelpunktte D,, D,. 

Aber auch die Asymptoten mit den Steigungen 


u 


F Man sind konjugiert imaginär und schneiden sich 


j in einem reellen Punkt, der ebenfalls auf der Geraden 
j gliegt. Ist nämlich G, ein Punkt von gin der Entfer- 
nung s von F, und sind g,,g, die Geraden mit den 

’ - Steigungen w/u durch @,, so sind die Parameterdar- 

| - stellungen von g,, 95 


@=8c0o8y+ ul, y=ssiny+ wi 


mit dem Parameter t. Setzt man diese Werte in die 
Gleichung der Koppelkurve (X) ein, so verschwindet 
bei beliebigem s der Koeffizient von i?. Für 


$ 2 ev 
5 20cosy— u 
| wird auch der Koeffizient von i? gleich 0. Die konju- 
giert imaginären Asymptoten g,, 9, schneiden sich also 
in der Tat in einem Punkt G, von g, der diese Entfer- 
= nung s von F\, besitzt, 
Er 4. Eine Beziehung zwischen den 
. Doppelpunkten unddenreellen Asym- 
ptoten-Schnittpunkten. Auf der Verbin- 
dungsgeraden g der Doppelpunkte D,, D, liegen nach 
‚obigem der reelle Schnittpunkt F, der Asymptoten mit 
fi den Steigungen -+i und der reelle Schnittpunkt @, der 
= Asymptoten g,, 9, mit den Steigungen w/u. Die vier 
A Punkte D,, D,, F,, 6, sind voneinander nicht unab- 
hängig. Wird vielmehr auf der Geraden g ein Abszis- 
.... sensystem beliebig festgelegt und sind d,, d,, fı, gı die 
h Abszissen der genannten 4 Punkte, so ist 


d,+%,=fıt 9 


d.h.: „Die Summen der von einem beliebigen An- 
fangspunkt auf der Geraden g gemessenen Abszissen 
der beiden Doppelpunkte D,, D, und der rellen Asym- 
_  ptoten-Schnittpunkte F,,@, sind einander gleich“, 
Wählt man nämlich als Anfangspunkt der Abszissen 

auf gden Punkt 7, und als + Richtung jene, die sich 
aus der Richtung der +z-Achse durch Drehung um 
+-y ergibt, so ist d,,, = 81,3, "ı = 0, g, = s und nach 
Art. 2,3 in der Tat 


9-8 


D 


zii 2 

‚2vcosy—u' 

kann diese Beziehung auch so formulieren: „Die 

i iden Punktepaare D,, D, und F,, @, haben denselben 

2 Mit ? unkt.‘“ u 

g: Ki Punkt @, ist also leicht zu konstruieren, wenn 
“ri ‚und F, bekannt sind. 

Dieke nnzeichnendenKEige nschaf- 

ı ne Bir emeinen Koppelkurve 


Y rbelgetriebes. Die allge- 


a < r 


Kleine Mitteilungen 


. 


389 


lare Kurve 4. Ordnung mit 2 Doppelpunkten D,, D,, 
einem singulären Brennpunkt F, auf der Verbindungs- 
geraden g der beiden Doppelpunkte und einem reellen 
Schnittpunkt @, der beiden nicht isotropen Asympto- 
ten, der ebenfalls auf der Geraden g und zwar so liegt, 
daß die Punktepaare D,, D,; F,, G, denselben Mittel- 
punkt: besitzen. ‘‘ 


6. Die Anzahl der zur Bestimmung 
der Koppelkurve verfügbaren Para- 
meter. Als algebraische Kurve 4. Ordnung ist 
die allgemeine Koppelkurve eines Schubkurbelge- 
triebes durch 14 einfache Bedingungen bestimmt. Daß 
sie zirkular ist, 2 Doppelpunkte besitzt und die reellen 
Schnittpunkte der beiden Paare konjugiert imaginärer 
Asymptoten auf der Verbindungsgeraden g der beiden 
Doppelpunkte liegen, bedeutet je 2 Bedingungen; dazu 
kommt die Entfernungsbedingung unter den Doppel- 
punkten und den reellen Asymptoten-Schnittpunkten 
(Art. 4). Es scheint demnach, daß zur Bestimmung 
der Koppelkurve nur 14—2—2—2—1 = 7 Bedingun- 
gen frei seien. Wählt man jedoch die Doppelpunkte 
D,, D, und die Asymptoten-Schnittpunkte 7,, @, ent- 
sprechend den im Begriff ‚‚Koppelkurve eines Schub- 
kurbelgetriebes‘‘ enthaltenen 7 Bedingungen und dazu 
die konjugiert imaginären uneigentlichen Punkte U,, 
U, der nichtisotropen Asymptoten, so sind wohl für 
die Koppelkurve 2+6-+4-+2 —= 14 Bedingungen vor- 
geschrieben; diese werden aber nicht bloß von einer 
Koppelkurve erfüllt, sondern von einem ganzen Bü- 
schel von Koppelkurven, das man als durch 2 zerfal- 
lende seiner Kurven bestimmt ansehen kann: Eine 
davon besteht aus der Geraden g und der uneigent- 


‚lichen Geraden, beide doppelt gezählt; die zweite zer- 


fällt in 2 konjugiert imaginäre Kegelschnitte k,, k,; 
k, ist bestimmt durch die Punkte D,, D,, /,, U, und 
die Tangenten g, = F,/,, kı =@1Uı,; k, durch die 
Punkte D,, D,, I,, U, und die Tangenten g, = F\,, Is, 
h,=G,U,. Daß die je 6 Bestimmungsstücke für 
jeden der Kegelschnitte %k,, k, miteinander verträglich 
sind, folgt aus der Entfernungsbedingung für die 
Punkte D,, D,, F}, @. Wird nämlich ein Kegelschnitt 
von 2 parallelen Geraden g, win den Punkten D,, D,; 
I,, U, geschnitten und sind F,, @, die Schnittpunkte 
von g mit den Tangenten in /,, U,, so haben die beiden 
Punktepaare D,, D,; F,, @, offenbar denselben Mittel- 
punkt. Dies bleibt auch richtig, wenn w mit der un- 
eigentlichen Geraden der Ebene identisch wird. — Die 
aus den beiden konjugiert imaginären Kegelschnitten 
kı,k, bestehende Kurve 4. Ordnung ist reell, wenn 
auch nullteilig. Es gibt also ein ganzes Büschel von 
reellen Koppelkurven, die die angeführten Bedingun- 
gen erfüllen; von ihnen geht eine und nur eine durch 
jeden Punkt allgemeiner Lage der Ebene. Zur Be- 
stimmung einer Koppelkurve eines Schubkurbelgetrie- 
bes können also außer den in dem Begriff einer solchen 
Kurve enthaltenen noch 8 einfache Bedingungen frei 
gewählt werden, oder es gibt oo® Kurven 4. Ordnung 
mit den für solehe Koppelkurven kennzeichnenden 
Eigenschaften. 

7.Eine Beziehungfür dieeine Kop- 
pelkurveerzeugenden beidenSchub- 
kurbelgetriebe. Die Bahnkurve eines beliebi- 
gen Punktes der Koppelebene eines Schubkurbelge- 
triebes #| ABFa% läßt noch eine 2. Erzeugung dieser 
Art vermittels eines Schubkurbelgetriebes F, A, PB, 
F,% zu ®). Die Geraden (B), (B,) mit dem Schnitt- 
punkt H seien die Kreuzkopfbahnen dieser beiden Ge- 
triebe (Bild 3), h die Gerade F, H. Dann besteht der 
folgende für die Getriebesynthese wichtige 

Satz: „Bei einer beliebigen Koppelkurve eines 
‚Schubkurbelgetriebes steht die Verbindungsgerade h 
des singulären Brennpunktes 7, mit dem Schnittpunkt 
H der beiden Kreuzkopfbahnen (B), (B,) normal auf 
der Verbindungsgeraden g der beiden Doppelpunkte 


j D 1, D 2°. 1 
es, Schubkurbelgetriebes läßt sich 


») W,Schmid,,,Über die mehrfache Erzeugung von Koppel- 


kennzeichnen: „Sie ist eine zirku- _ kurven“, Z. angew. Math. Mech. 30 (1950), 8. 330—333. Art, 2. 
> et Pe, L 2 1 


1 
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Sg 3,2, —7,8, für 6, s— — 4,6. 
Dresden. _. W. Schmid. 


fen 
Bild 3 : 


die Steigung der Kreuzkopfbahn (B,) gegen die »- Grenzschichtströmung kurz nach dem plötz- 
Achse F} F, „ gleich cotga, daher ihre Gleichung _ eher An bzw. Abstoppen eines rotieren- 
en Bodens. 


a a RT Die Feststellung eines Rechenfehlers in der früheren 


Die Gleichung von (B) ist Arbeit des Verfassers über dasselbe Problem [1], auf 
den freundlicherweise Herr Prof. H. Görtler (vgl. 
ze) auch [4]) aufmerksam machte, gab Veranlassung} eine 
der Schnittpunkt H von (B) und (B,) hat demnach Neuberechnung vorzunehmen. - 
die Koordinaten Um einen unmittelbaren Vergleich mit den Mes- 
sungen im rotierenden Laboratorium des Kaiser- 
ee are 3 Wilhelm-Instituts für Strömungsforschung durch- a 
Le aim a Per führen zu können, wurde die Flüssigkeitsbewegung als = 
Be . i Relativbewegung gegen den ruhend gedachten Boden 
Die Steigung der Geraden h = F}H ist also zurückgeführt. Mit den Zylinderkoordinaten s,gundz 
er und den Geschwindigkeitskomponenten « in tangen- 
SUPER tialer, vin radialer und win vertikaler Riehtung zu der j 


der inkompressiblen reibenden Flüssigkeit über dem 
% | Boden angegeben werden für den Fall 


+’ 


zu Null werden: 


und dies ist (Bild4) gleich — cotgy, während g die 
Steigung tg y besitzt. Es ist also in der Tat; } normal 
zu g. SE ER u 


har *) Vgl. W. Schmid, a.a.0,, Art.l,2. 


nf % 
er’. 


4 


ee 


Die Randbedingungen lauten für >20: 


BORN =. :=6,0/= 6,0) = 0), 
lo) =—1, lm. 
P, (oo) = 2 .6,(00) = G, (00) = 0) 


Zur Berechnung von w machen wir den Ansatz: 
w=— 40 ri (wie, + NM +: an. 


Die in (1) und (2) eingesetzten Ansätze (4) und (5) lie- 
fern nach durchgeführtem Koeffizientenvergleich be- 


züglich steigender wt-Potenzen die Differential- 
gleichungen: 

F+27B—=0 28% ER EN TR AR (9), 
a er ZB Bm de ao), 


S Y : 


3 RR, ho 
wo H,=/f @,dn gesetzt ist. 
0 


Die Beiwertfunktionen F,, G,, H, und F, wurden | 


geschlossen integriert und ergaben unter Berück- 
Sack) der BERNER „9 und (7), wenn 


fer dn= a, fer and, ei dn=e und, 


2 =.a RAR wird, die RN 


a\ 3/ Ur 
U2 N Pr 10,4 827,3 
er 145)’ mn na? 
3 sr | Er 3 5 (14) 
ERS E BR ET En A Te 
SIT EEE AT ET 2 2 
? lan 3)" Falls, ( 2) 
, 25 ö Up 1er j 2 Pezie 
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a hm r Na? — — Ma — 0° 
3 27 
u ur gl ANETTR 2 
Bi 2 16 7 aa 
r } 3 Ast 
+ an n°a 3 a En na 3_— zab+ gone 
243 243 / 3 
BTW + + a I, a 
A 
It zum mn gn „>aa’ 
15 15 
! — gmranaa + er = nta 
B12 39 9 51 
Bu SNEEK 2 
Es =) Pe ke Fjg n) ? 
Low — Mair eye + Be = n? a’ 
2n 2n31l2 128" 647 
ONdT 1 
F3 (2: a 2) nat nel 2 
64 ; 
+ 572 ( +37 
ul. 078. 68 (15) 
Ya \36 2m 5n 157? 


f 3 
x (Mar swa+ ger gut za) 


Bi en nal ——— n?a® 
a? — EEE 
as! 2Ya 
ge en) 2 
322 480 


+) re 


nad 


32 160 


3, Tn— 128\, 
128"  640n 
@ En) Ss 


64 F 960 7 


+ 
N. 1, 
+ (mt BETT ) 2 : 


I wi 
*( ann 3 u | 


30Y 8 20 10Yr 


2. Die: tionen sind aufgetragen in Zahlentafel 1 und 
" dargestellt i in Bild I, wobei 3 Stellen gesichert, die 4. 


nicht ganz sichergestellt ist. Die Berechnung von F, 


5 Filste so, daß das inhomogene Glied der Differential- 
Be an) nach Einsetzen der Beiwertfunktionen 


@, H,in den Bestandteil des Auslaufs und eine 


2 Restsumme aufgespalten wurde, Zudem wurde die 
Auslaufdifferentialgleichung nochmals in 2 Teilen inte- 
. griert, was in den eckigen Klammern zum Ausdruck 
ee, Die Probe durch Rückdifferenzieren zeigte, 

BR die Lösung jeweilsrichtig war. Da die Differential- 
 gleichung linear ist, können die Teillösungen zur 


esamtlösung zusammengesetzt werden. 
‚us diesen Ergebnissen erhält man die Neigung m 


den logarithmischen Spiralen r = e”"? der Strom- 
Fa Boden i in zweiter Ordnung bezüglich ot: 
4 86 (0) 

0) + (et)? F; (0) 
0,7268 @1 ®) 

I + 0,167 (oi)% 


(16). 


Zahlentafell. BE a Elton en 
des Anlaufproblems 


N F,(n) | 6, (N) H, (N) | F,(n) 
0 0 0 0 0 
0,1 —0,1125. | 0,0639 0,0035 — 0,0187 
0,2 — 0,2227 0,0975 0,0118 — 0,0331 
0,3 — 0,3286 0,1099 0,0219 —0,0422 
0,4 — 0,4284 0,1084 0,0333 — 0,0474, 
0,5 — 0,5205 0,0984 0,0436 — 0,0489 
0,6 — 0,6039 0,0846 0,0529 — 0,0466. 
0,7 — 0,6778 0,0692 0,0606 — 0,0428. 
0,8 — 0,7421 0,0547 0,0667 — 0,0379, 
0,9 — 0,7969 0,0418 0,0716 —0,0326. 
- 1,0 — 0,8427 0,0311 0,0752 — 0,0272 
sl. — ‚08802 0,0225 0,0778 — 0,0222 
1,2 —0,9103 0,0158 0,0797 | —0,0178. 
13 — 0,9340 0,0110 0,0811 — 0,0138 
1,4 |. —0,9523 0,0077 0,0820 — 0,0106, 
1,5 — 0,9661, 0,0049 0,0825 — 0,0079 
” 1,6 —0,9763 0,0028 0,0827 —0,0060. 
; 17 — 0,9838 0,0022 0,0831 — 0,0044 
1,8 — 0,9891 0,0015 0,0834 — 0,0030 
1,9 —0,9928 0,0009 0,0835 — 0,0020 
2,0 — 0,9953 0,0004 0,0836 — 0,0013. 
2,1 — 0,9970 0,0002 0,0836 — 0,0007 
2,2 — 0,9981 0,0000 0,0836 — 0,0004 
2,3 —0,9989 0,0000 0,0836 —0,0002 
2,4 — 0,9993 0,0000 0,0836 —0,0001. 
2.5 —.0,9996 0,0000 0,0836 —.0,0000 
2,6 — 0,9998 0,0000 0,0836 —0,0000 \ 
Der von Cochran [2] angegebene Wert für große 
Zeiten nach Bewegungsbeginn ist m = — 0,828. Die 
erste Näherung m = — 0,7268 ot und die zweite 
Näherung (16) sind in Bild 2 aufgetragen. h 


0,9 
-m De er 
08 vH 

07 

» 

06 R 
an —— —— I.Näherung 

‘ Z.Näherung 


nz vermuteter wahrer Verlauf 


0,5 20 15 20 235 


wr 


II, Näherung, verglichen mit dem wahrscheinlichen Verlauf auf 
Grund des asymptotischen Wertes nach Arien 


b) Auslauf. Die Ditterentialgleichungen lauten 


E jetzt: 
a ET (17) 
N II ARLEREN. Re (19) 


unter den Randbedingungen für £ > 0, daß fürz = 0: 
vev=w=0;füre = 0; u=or,» = 0. Weiter ist 
Die Beiwerte Re Auslaufs werden 


IE 2 a 
r 


2) Es ergab Blche in dieser Rechnung 40107 statt 0,23 in der 
üheren. A ‘ 


Parse F, siehe geschweifte Re (15) und die fol. 


Bild 2. Stromliniensteigung am Boden für den Anlauf in I. und 


Mniler Benutzung des Ansatzes @ und: (5% zur Unter- 
dung vom Anlauf mit einem Querstrich versehen: Ex 


F,(0) = F (0), = rg) AUBTES 0,0) =::-=0 (20); - j 
f F,(x) =1, yl . n 

Zn m a HR (21). 

Flo) =: = 6,(00) = G,(&) = 1 


Ansatz (8) gilt wieder, jedoch für die überstrichenen 


Beiwerte. Die Differentialgleichungen werden dann: 
EP 24, = 0 A (22) 
+ 276 — 46, Rn ei Er 


ta BF SR N 


wo H,= 2 fe G, dn gesetzt ist. 


Die Beiwertfunktionen wurden WiodbR. geschlossen 
integriert und ergaben unter Berücksichtigung der 
Randbedingungen (20) und (21) mit den gleichen Ab- 


kürzungen a, b und c wie oben: 3 
De L LAN Ne (25): 
Ya E 
LT 8 RE > 
ee et am 35 aa | h 
e (27) 
2 =) 2 { N 
re, meter) a 
e 8. 8 4 i F 
u a RR TERN zer : “ } ; 
8 a ae 
taten ae 1 
a7 de ver pre u 2 
2 AuTE EAN B 
Fa ger 372 a 


Yes, r 


Zahlentafel2, Beiwertfunktionen 
des ERLERNT 


genden Ausführungen Dr 


LEERE ERDE 
oo EUER | 
0,1. | 0,1125. | 0,1240 | 0,0065 | 0,0096 

0,2 18 | —0,0236 | 0,0184 
03 | 0,0478 | 0,0263. 
04 —0,0761 | 0,0531 
706% | 0,1061 | 0,084. 
0,6 0.1357 | 00008. FIR 
07 0, ‚1636 0,0405 % 
{ 0,8 j e h al 

03 

a 

LTE 

1,2 
Ws 

1 

155." 


a) 


% 7r Diese. unktionen and ee in Zahlentafel 2 


und dargestellt i in Bild 3, wo 3 Stellen gesichert, die 4. 


nicht, es grobergestellt ist. Im Gegensatz zur früheren 


Bild 3. Beiwertfunktionen des Auslaufs 


f 
’ 


Arbeit wurde auch hier mit dem richtiggestellten in- 
homogenen Glied der Gleichung (24) statt numerisch 
in geschlossenen Ausdrücken integriert. 

Die Stromlinien am Boden ergeben sich jetzt in 
zweiter Ordnung zu: 5 


RR ERBETEN 1,2732 ot, 
[pr 1 + 0,0794 (ot)? 


LZI) 
m (29) 


Mit der ersten Näherung m = — 1,2732 wt ist die 


zweite Näherung (29) in Bild 4 aufgetragen. Der Wert 
-+0,0794 (ot)? gibt einen besseren Verlauf gegen den 
von Bö de wa dt [3] kurz nach der Arbeit [1] ge- 


Bodewadt 


} u 
.— 
.— 


-—---/ Naherung 
IT Näherung 
—:—— vermufefer wahrer Verlauf 


X 


terung, Feen mit dem wahrscheinlichen Verlauf auf 
Bruns des asymptotischen ee usoh Bödewadt 
na ö 


= 2 1218 als der 5 


0.1 40), 8. 1 bis 


Son. Bd, 1:30 01934), 


N Kleine Mitteilungen 


4P,,, (2) = — 


Avomlisiknsfsisuig am Boden für den Buslauf in I. und Rn : (3) 
Apyula) = hy! ("= ») (3 BES ) 


I Für die Abweichung y (x) — p,, 
ch <r<hgilt‘!) 


bang vo == Pg1,0, x 


flow Be oa ta : 
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[3] U.T. Bödewadt: Die Drehströmung über 
festem Grunde. Z. angew. Math. Mech. Bd. 20 
(1940). S. 241—- 253. 

[4] H. Görtler: Zähe Flüssigkeiten. Naturfor- 
schung und Medizin in Deutschland 1939— 1946, 
Teil ILL. angew. Math. Mech,, 5, S. 69. 


Friedberg i./H. K. -H. Thiriot. 


Notiz zur Stabbiegung. 


Ein gerader, horizontaler, an seinen beiden Enden 


vw =— + und > gelenkig gelagerter Stabmitder _ 


konstanten Biegesteifigkeit EJ sei vertikal belastet 
(siehe Bild 1). Die Belastung g(x) sei vorläufig stetig, 
aber sonst beliebig vorausgesetzt. Sie werde wie üblich 


Bild 1 


positiv gerechnet, wenn sie in Richtung der lrd- 
schwere wirkt, desgleichen die Durchbiegung y (x). 
Letztere bestimmt sich bei den Randbedingungen 


y & a = Y (2)=% y" B =) (3) =0 (1) 


als viermal stetig differenzierbare Funktion von & 
eindeutig aus der Differentialgleichung 
EIyHla)= g(2). er to 
Es sei p, ,„,(2) (mit ganzzahligem &,,a, > 1) die ein- 
deutig bestimmteParabel (höchstens) (&, + &,— 1)-ter 
Ordnung in z, welche die Biegelinie y (x) in 
h.bzw. x = Tg h &,-punktig bzw. 
&,-punktig berührt. p, „, (2) hat also mit y (x) bei- 
spielsweise in x = — h sämtliche Ableitungen bis ein- 
schließlich der (&, — I)-ten gemeinsam. 
Mit den Abkürzungen 


l ’ ’ U 4 
v(5)=# v( 2)= 4 2h=l 


} erhält man für die Paare (a,,&,) = (2,2), (3,1) (1,3) 
bei Berücksichtigung von (1) die Parabelgleichungen 


w” j 
\ ( = ») (n a 


® } 
try ) 


At) = (1) ” 


‚(a4 Ma m Yard (M, 


ti Lan)i 2, 


- 4) A. Markofi: Differenzenrechnung (deutsche Über- 
Bunt ap Th, Friesendorff und E. Prümm). Leipzig 1896, 
ns 


4‘ 


a, (2) im Intervall a 


uk ‚ ” N 


- 


‚schen Physik während der letzten Jahrzehnte „eine 
zunehmende Entfremdung zwischen diesen beiden 


wo z eine gewisse, uns hier näher nicht interessierende 
Zwischenstelle mit -h<z<h ist. 


Wir machennun weiter diefürs Folgende wes°ntliche 


Voraussetzung, daß die stetige Flächenbelastung 
durchweg gleichgerichtet ist. Es sei z. B. q(x) 
durchweg niehtnegativ, also für -h<x<h 


4(@)20. . (5) 
und natürlich 
1) #0. (5a). 
Dann folgt aus (4), (2), (5) für -h<a<h 
Poa(2) <y(z) <P31(®), Punsl®) - . (6). 


Die Parabeln p,,.(2) bzw. P3,1(2), P1,3(2) schließen 
also (siehe Bilal) im Intervalle —h<xr<h die 
Biegelinie y(x) ein. Nach Voraussetzung schneiden 
sich 91,5(2), P3,1(2) in = — h und @ = h. Die Ab- 
szisse ©* ihres dritten Schnittpunktes liegt im Innern 
des Intervalles (— h, h). Dies folgt aus der Anschau- 
ung oder aus den GI. (4), (5), (52). Unter Verwendung 
von (3) erhält man 


w* 
1 < rs z v2 
h Yan Yy_ 


(N, 


(Bei unseren Voraussetzungen ist P1,3(2) = P3,1(%) 
nicht möglich, so daß es immer genau eine solche Stelle 
x* gibt.) 

Nunmehr folgt aus (6) 


P3,1(8) < P1,3 (8); 
lc RS Br n.| 


(8) 
Pı»3 (2) < Par1(®), | 
( h<a*<a<h) 
Aus (7) und (8) lassen sich in einfacher Weise einige 
Folgerungen ziehen, 

a) Der Schnittpunkt S der Tangenten an die Biege- 
(siehe Bild 2)! hat die 


Paa(a) Syla) S | 


linie in Ah und a=7 
Abszisse a 


Bild 2 


Dr. phil, R, Becker (ord. Prof. a. d, Univ. Göttingen) 
Vorstufe zur theoretischen Physik. VII u. 
172 S. mit 94 Abb. Berlin/Göttingen/Heidelberg 
1950. Springer-Verlag. Preis kart. 7,50 DM. 

Das Buch ist aus dem Empfinden heraus entstan- 
den, daß trotz der innigen Verknüpfung von Mathe- 
matik und Physik in der Entwicklung der theoreti- 


Te: 


u. 7a x / 
Buchbesprechungen 
= E - 


Wegen (7) gilt also 


El. 9. 
b) Aus (8) erhält man, wenn man noch die spitzen 

Winkel &_, + (siehe Bild 2) einführt, also 
y=tgo-, y =— tgo+ und weiter y(0) = yo 

setzt, für y, die Ungleichung 3 


e 3 
z wo trigo.)<y<z 


Bei zu x = O symmetrischer Belastung folgt aus (10) 


ER; 


_ 


I Min (tg ©-, tg ©+) (10). 


UA) 


N 


I 3 
4 Wwo-<sy<sz itgo- 


— 


sowie die Näherung 


5 Be 
Hnrrgitgo-. A 


Sind die Eindtangenten zeichnerisch gegeben, so führt 
die Ungleichung (10) auf die folgende Konstruktion 
zur Eingrenzung von HE (siehe Bild 2): n 


Durch die Punkte » = + 4 der »-Achse ziehe man zu 


den gegebenen Tangenten die Parallelen. Man erhält 
so auf der Geraden x = 04 Schnittpunkte O, U, 0’, U’, 
Dabei gehöre zur steileren der beiden Richtungen der 
Buchstabe U, zur anderen der Buchstabe O0. Es sei 
M gry, die Mitte der Strecke 0’ U’, M yo die Mitte von 
0'0. Dann schneidet die Biegelinie die Gerade x = 0 
in einem Punkte von Myy Moro- — Falls 0; = o_ 


ist, fallen O und U, desgleichen 0’, U’, M y,y, zusam- 
men, Man erhält so die Eingrenzung entsprechend 
(10a). f > 
Aus (10) folgt weiter 
tgo- <tgo; <2tgu ; 
= ss Ss (11). 
nt <tg w_ <2tg + ef; 


Man bestätigt leicht, daß (11) gleichbedeutend ist 
mit (9). | 7. IRRE 


c) Die für stetiges q (#2)> 0 aufgestellten Unglei- 
chungen (9) und (10) bleiben auch bei unsteti- 
gem qg(x) sowie im-Falle des Auftretens von konzen- 
trierten Einzellasten gültig, falls sämtliche 
Lasten vertikal und gleichgerichtet wirken und de 
Diskontinuitäten nur in endlicher Zahl auftreten. ’ 
Wegen der Gültigkeit des Superpositionsprinzips ge- 
nügt es, die Behauptung für den Fall einer unstetigen 
Belastung des Stabes durch eine Einzellast und eine 
unstetige Flächenlast nachzuweisen. Man kann nun 
letztere Lasten beliebig genau durch stetige nicht. 
negative Flächenlasten annähern. Daher treffen di 
für’ solche Näherungen gültigen leichungen 
und (10) auch in der Grenze zit 


von einem kurzen Exkurs in die Elektrostatik werden 
die behandelten Beispiele der Mechanik und der Wär- 
melehre entnommen, den beiden Teilgebieten also, 
die sich durch eine besondere verschiedenartige Be- 
handlung der Probleme auszeichnen, so daß dem Leser 
auf knappem Raum eine große Auswahl physikali- 
scher und mathematischer Fragestellungen nahe- 
gebracht wird. In der Mechanik wird zunächst die 
Mechanik eines und vieler Massenpunkte behandelt, 
und dann werden die Längsschwingungen einer line- 
aren, elastischen Kette mit dem Grenzfall des elasti- 
schen Stabes ausführlich diskutiert. In der Wärme- 
lehre nehmen Beispiele aus der kinetischen Gastheorie 
den breitesten Raum ein. Es schließt sich ein weiterer 
Teil ‚Mathematische Erinnerungen und Beispiele‘ 
an, der in seiner Unvollständigkeit und Knappheit 
nicht ganz so überzeugend wirkt wie die Hauptteile 
des Buches. Hier wäre eher ein Hinweis auf die ein- 
schlägige mathematische Literatur zu empfehlen. 
Die mathematischen und physikalischen Voraus- 
setzungen entsprechen im allgemeinen dem Bildungs- 
stand eines 3. oder 4. Semesters, doch werden an 
einigen Stellen in dem Bestreben, den Leser an die 
moderne Betrachtungsweise heranzuführen, auch Tat- 
sachen gebracht, die vorläufig einfach hingenommen 
werden müssen. Besonders hervorzuheben ist die 
klare, einprägsame Sprache, die den Leser oft persön- 
2 lich anspricht und ihn zum eigenen Mit- und Weiter- 
denken veranlaßt; und das ist zum rechten Verständ- 
nis dieses trotz des bescheidenen Titels recht an- 
spruchsvollen Buches unentbehrlich. 


j Dresden. G. Opitz. 
j 
Prof. Dr. Richard Gans, Vektoranalysis 
| mit Anwendungen aus Physik und 
| Technik. (Teubners Mathematische Leitfäden 
Bd. 16), siebente Aufl. 1208. m. 44 Abb. Leipzig 
| 1950. B. G. Teubners Verlagsgesellschaft. Preis kart. 
5,90 DM. ’ | 
Die sechste Auflage dieses 1905 zuerst erschienenen 
Bändchens wurde in Bd.10 (1920) Seite 308/309 
I dieser Zeitschrift angezeigt. Der Inhalt dieser sieben- 
ten, von Dr. Stein durchgesehenen Auflage ist im 
wesentlichen unverändert geblieben. Das Bändchen 
‚ unterscheidet sich von manchen anderen einführenden 
Büchern dadurch, daß die Beweisführung in größerem 
Umfang rein vektoriell ist. Die Anwendungen sind 
in der Hauptsache dem Gebiet der Elektrodynamik 
- entnommen. . 


Dresden. Willers. 


Fr. A, Willers, Methoden der praktischen 
Analysis. (Göschens Lehrbücherei Band 12); 
Br 2. Aufl. 410 S. m. 93 Abb. Berlin 1950. Walter de 

Gruyter & Co. Preis geb. 24,— DM. 

Es gibt wohl kaum einen Mathematiker, dem beim 
Studium der numerischen, graphischen und instru- 
mentellen Verfahren das Lehrbuch von Willers unbe- 

kannt geblieben ist und der es nicht wegen seiner 
 Reichhaltigkeit und vorbildlichen Klarheit schätzen 

gelernt hat. Darüber hinaus ist eseinem großen Kreise 
von Ingenieuren und Physikern ein zuverlässiger 
Berater in der Praxis ihres Berufs geworden, Wenn 
nun das Buch über zehn Jahre vergriffen und bereits 
sehr gesucht war, so ist das lange Ausbleiben einer 
Neuauflage allein den widrigen Verhältnissen der 


-  fasser hat inzwischen sein Werk einer gründlichen 
_  Neubearbeitung unterzogen. Einige veraltete Ver- 
fahren und solche mit beschränktem Anwendungs- 
bereich wurden fortgelassen. Dafür ist aber jedes 
- Kapitel durch Zusätze und Literaturhinweise be- 
rt worden, 
die neueste Zeit berücksichtigen. Trotz dieser 
Br nhalts ist 
uch i n Sinne des Wortes geblieben 
seiner Frische eingebüßt. Wenn 
A Yen 29 r ey 


4 : Buchbesprechungen 


Kriegs- und Nachkriegszeit zuzuschreiben. Der Ver- 


‚die alle wesentlichen Fortschritte 


das Buch gut lesbar und 
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uns der Verlag jetzt die zweite Auflage auf gutem 
Papier gedruckt zu einem angemessenen Preise vor- 
legt, soist esein Wiedersehen mit einem guten Freund, 
auf dessen erfahrenen Rat wir nicht mehr ‚verzichten 
mögen. 

Inhalt: 1. Das Zahlenrechnen und seine Hilfsmittel. 
2. Interpolation im engeren Sinne. 3. Angenäherte 
Integration und Differentiation. 4. Interpolation im 
weiteren Sinne; trigonometrische Interpolation. 
5. Praktische Gleichungslehre. 6. Genäherte Integra- 
tion von Differentialgleichungen. 1 


Aachen. Günther Schulz. 
Dr, phil. Hans Hornich (ord. Prof. a. d. Techn. Hoch- 
schule Graz). Lehrbuch der Funktionentheo- 
rie. VII u. 2168. m. 34 Abb. Wien 1950. Springer- 

Verlag. Preis 19,50 DM, geb. 21,60 DM. 
In wohl ausgewogener Auswahl aus dem Gesamt- 
gebiet der Funktionentheorie hat der Verf. in seinen 
Vorlesungen die klassischen Bestandteile und Grund- 
pfeiler jeder Einführung in das Gebiet knapp und 
treffend dargestellt. Er ist dabei mehrfach Wege ge- 
gangen, die von dem Vorgehen früherer Vorlagen deut- 
lich abweichen; wir erwähnen etwa die Einführung der 
analytischen Fortsetzung auf Grund eines Äquivalenz- 
prinzips, die naheliegend, aber wohl-bisher in Lehr- 
büchern nicht mit solcher Folgerichtigkeit entwickelt 
wurde. Gegenüber der kurz zuvor erschienenen Dar- 
stellung von Heinhold betont Hornich stärker die 
den Universitäten eigentümliche Seite der Funktionen- 
theorie; wenn dabei das anschauliche Moment zu- 
rücktritt, wie es in vielen Figuren über konforme Ab- 
bildungen dem Ingenieurmathematiker näher liegt; 
so muß doch betont werden, daß das Buch auch für 
diesen sehr wertvoll ist, weil es die unentbehrlichen 
Grundgedanken ohne zu viel Belastung mit Theorie so 
darstellt, daß sie rasch zu handwerklichem Gebrauch 
eingesetzt werden können. Die einzelnen Kapitel be- 
treffen: Komplexe Zahlen, Differenzierbare Funktio- 
nen, Potenzreihen, Integrale im Komplexen, Satz von 
Cauchy, Isolierte Singularitäten, Reihen von Funk- 
tionen, Analytische Fortsetzung, Untersuchung spe- 
zieller Funktionen, Algebraische Funktionen und ihre 
Integrale. — Das letzte Kapitel gibt eine kurze, aber 
vorzügliche Darstellung der Grundgedanken der Theo- 
rie der algebraischen Funktionen und der Abelschen 
Integrale (ohne aber auf volle Beweise mit arithme- 
tischen Mitteln einzugehen); und es umfaßt einen Aus- 
blick auf die elliptischen Funktionen und Integrale. 
Große Theorien sind im vorletzten Kapitel etwas iso- 
liert vertreten durch den Riemannschen Abbildungs- 
satz und den Picardschen Satz; der Ingenieurmathe- 
matiker wird das für die konforme Abbildung etwas 
"bedauern; die Beweistechnik erscheint dort — not- 
wendig ohne größeren Zusammenhang wiedergegeben 
— etwas isoliert. — Die Druckerei war offenbar ihren 
“ mathematischen Aufgaben noch nicht ganz gewachsen, 
‚worunter die Formelbilder öfter leiden. — Eine An- 
zahl trefflich gewählter Übungsaufgaben ist anerken- 
nend zu nennen. In allem: Ein sehr gutes Buch. 


Gießen. E. Ullrich. 


Dr.-Ing. Hans Umstätter (Kaiser-Wilhelm-Institut 
für physikalische Chemie und Elektrochemie). Struk- 
'turmechanik. Ein Beitrag zur Physik der Kolloide. 
XVI u. 237 S. mit 1 Titelbild, 143 Abb, u. 34. Tab. 
Dresden und Leipzig 1948. Verlag ‘Theodor Stein- 
kopff. Preis: br. 17,— DM., geb. 19,— DM, 

Der Verfasser hat sich mit der vorliegenden Mono- 
graphie die Aufgabe gestellt, die mannigfaltigen Er- 
scheinungen und Gesetzmäßigkeiten der Kolloid- 
physik unter einem einheitlichen Gesichtspunkt nicht 
nur darzustellen, sondern darüber hinaus ihnen eine 
‚theoretische Grundlage zu geben. Da die Kolloide 


_ infolge ihrer „‚Struktur‘“ den Bedingungen der klas- 


 sischen Kontinuumsmechanik nicht mehr gehorchen, ' 
a es erforderlich, diese entsprechend zu modifizieren, 
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bzw. zu ergänzen. In diesem Sinne bezeichnet der 


Verfasser und stellt als die drei wichtigsten Postu- 


late der Strukturmechanik an die Spitze: die Stetigkeit 

der Zeit, die Quantelung des Raumes und die Endlich- 

keit der Fortpflanzungsgeschwindigkeit der Reibungs- 

störungen. Die durch die Struktur bedingte Quante- 

lung des Raums führt rechnerisch zur teilweisen Ein- 

führung von ‚‚Extremalwerten“ Az, Ay, .... in die 

Grundgleichungen an Stelle der Differentialwerte d, 

dy,...., eben wegen der notwendigen Forderung, daß 

die Strukturelemente unveränderlich sein sollen, wenn 

sich die Qualität der Substanz nicht ändert. Dies gilt 

z. B. für den Gradienten einer Zustandsgröße (Ge- 

schwindigkeit, Konzentration, Temperatur) Azx/Ay. 

Es bleibe dahingestellt, wie weit man im Einzelnen 

| der mathematischen Behandlung und den gezogenen 

Folgerungen zu folgen geneigt ist. Zu würdigen ist 

jedenfalls der Versuch, auf solche Weise ein einheit- 

liches Gebäude zu schaffen. Zu betonen ist ferner 

noch die Reichhaltigkeit des Inhalts, der stets das 

En Bestreben zeigt, die gezogenen Folgerungen an den 

| zahlreichen bereits vorliegenden Beobachtungsergeb- 

' nissen der Rheologie u. a. zu prüfen. So darf man im 

{, ‚. Ganzen das Buch als einen zwar vielfach kühnen, oft 

nicht ganz durchsichtigen, aber trotzdem durchaus 

anzuerkennenden Vorstoß in dieses schwierige Gebiet 
bezeichnen, 


Weilburg (Lahn). L. Schiller. 


3 Veröffentlichungen des Bamberger Instituts für 
Be) -  Erdmessung, Bamberger Verlagshaus Meisenbach & Co. 
I Bamberg 1949. 

i Band 1: Das Zentraleuropäische Dreiecksnetz, 
Grundlagen. 122 S. und 5 Tafeln. Preis 
> 4,— DM. 

! Band 3: Professor Hans Boltz und sein Werk. 
- 77 8., Preis 3,— DM. 
n Band 4: Beiträge zur Ausgleichung astronomisch- 
$ geodätischer Netze unter ‚besonderer Be- 
rücksichtigung des Zentraleuropäischen 
Netzes. 143 S. mit 6 Abb. Preis 4,— DM. 
Beiträge zur Lotabweichungsausgleichung 
und Geoidbestimmung. Allgemeinere Be- 
trachtungen und Durchführung im Bereich 
des Zentraleuropäischen Netzes. I. Teil: 
106 S. mit 6 Abb., 2. Teil: 25 Tafeln. 
Preis 5,— DM. 

Beiträge zur Zentrierung von Richtungen. 
47 8. mit 9 Abb. Preis 3,— DM. 
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Band 9: 
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AM } einer Schriftenreihe, mit der das von Direktor Gigas 
hiesk > geleitete Bamberger Institut für Erdmessung über seine 
RN - im Jahr 1945 begonnene Tätigkeit berichtet. In außer- 
ordentlich kurzer Zeit ist u. a, bereits eine umfang- 
reiche Arbeit von großer internationaler Bedeutung ge- 
leistet worden: die Ausgleichung des Zentraleuro- 
päischen Dreiecksnetzes (ZEN). 

Die in den Bänden 1, 3, 4, 6 enthaltenen Veröffent- 
lichungen bilden insofern eine in sich zusammen- 
hängende Gruppe, als sie sich mittelbar auf diese Arbeit 
bzw. die bei derselben benutzten Boltzschen Aus- 
gleichungsverfahren beziehen, 

Der Band 1 bringt nach einem Geleitwort von 
Gigas eine von Levasseur verfaßte allgemeine 
Übersicht über die Ausdehnung, Anlage und Aus- 
gleichung des ZEN. Daran schließt sich eine grund- 
sätzliche Studie von Wolf über die Ausgleichung 
astronomisch-geodätischer Dreiecksnetze mit beson- 
derer Berücksichtigung des ZEN. Gewisse Einzel- 
probleme, die sich bei der Ausgleichung des ZEN aus 
dessen Eigenart ergaben, werden in den Bänden 4 und 6 
von Wolf und Berroth in einer ganzen Reihe 
besonderer Aufsätze ausführlicher behandelt. Im 
2. Teil des Bandes 6 ist in 25 Tafeln das auf die ein- 
zelnen Arbeitsgänge und die Ergebnisse der Ausglei- 
“ chung des ZEN bezügliche zahlenmäßige Material in 
übersichtlicher Weise zusammengestellt. Der Band 3 
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Diese Bände bilden den verheißungsvollen Anfang - 


bringt nach einem Nachruf von Gigas auf denim. 
Jahre 1947 verstorbenen Prof. Dr.Hans Boltz eine 

Studie von Wolf, in der die Stellung des Boltz-_ 
schen Entwicklungsverfahrens in der Ausgleichungs- 

rechnung nach bedingten Beobachtungen behandelt 

wird. t 

Das ZEN ist ein gitterförmig gestaltetes Rahmen- 
netz, das im Norden durch die Nord- und Ostsee, im 
Osten durch die Ostgrenzen Lettlands und Polens, im 
Süden durch den Parallelkreis in 48° Breite und im 
Westen durch den Rhein begrenzt wird. Das geodä- 
tische und astronomische Beobachtungsmaterial ist aus 
selbständigen Haupttriangulationen der vom ZEN er- 
faßten Länder entnommen. Als Bezugsfläche ist das 
Hayford-Ellipsoid gewählt. 

Die bereits im Jahr 1947 zum Abschluß gebrachte 
Ausgleichung I des ZEN stellt einen Mittelweg zwischen 
der Bowie- Methode und einem neuen von Wolf 
entwickelten Verfahren vor, das der strengen Aus- 
gleichung nach Eggert gleichkommt,, 


Die Bo wie - Methode und in der Sowjetunion von 
Krassowskij und Urmajew entwickelte Ver- 
fahren haben als Vorbilder gedient. Es sind jedoch 
verschiedene grundsätzliche Ergänzungen und Abän- 
derungen eingeführt worden. Durch eine besondere 
Lotabweichungsausgleichung haben die Ausgangsdaten 
im alten Zentralpunkt Potsdam (Helmert- Turm) v 
Verbesserungen erhalten. Die Laplaceschen Azi- E 
mute und dıe Grundlinien mit ihren Vergrößerungs- j 
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. netzen wurden nicht als fehlerfreie Elemente behandelt, 


sondern durch besondere Azimut- und Basisaus- : 
gleichungen mit Verbesserungen versehen, die zur Til- : 
gung der auf die Azimut- und Grundlinienbedingungen 
bezüglichen Widersprüche beitragen. Die Gewichte für j 
die bei der Schleifenausgleichung benutzten Breiten- 
und Längenunterschiede wurden als Funktionsgewichte 4 
aus den Verbindungsketten entnommen. Die Schleifen- Fe 
ausgleichung erfolgte nach bedingten Beobachtungen. | 
Nach Abschluß der Ausgleichung I ist eine noch- 
malige Ausgleichung II des (durch Einbeziehung von 
zwei neuen Verbindungsketten etwas veränderten) 
ZEN nach dem neuen von Wolf entwickelten Ver- 
fahren ‚durchgeführt worden. 
Ebenso wie beim Verfahren von Eggert werden 
auch bei dem von Wolf die Bo wieschen Knoten- 
netze durch einfache Seiten, Knotenstrahlen, ersetzt, 
so daß die Verbindungsketten alle Richtungsmessungen 
enthalten. Beim Heranführen der Verbindungsketten 
an einen solchen Knotenstrahl ergeben sich je ein bis 
zwei den übergreifenden Ketten gemeinsame Dreiecke, 
deren Identität in den betreffenden Ketten durch Ein- 
führung von je zwei Winkelidentitätsbedingungen ge- 
‚sichert wird. Diese Winkelidentitätsbedingungen wer- 
den mit den Azimut-, Grundlinien-, Breiten- und 


gemessenen La pla ceschen Azimute und die G 
linien werden aber, ebenso wie bei der Ausgleichung 
rfreie, sondern als verh 
rungsbedürftige Elemente I a 
Mit der Ausgleichung. 
Wolf entwickelten Ve 
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. Math, Mech, 
Ba, 30 Nr. 11/12 Nov./Dez. 1950 
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Dank seiner Lage, dem umfangreichen Beobachtungs- 
material und seiner ho:hentwickelten Ausgleichung er- 
scheint das ZEN sehr wohl als Anschlußblock für wei- 
tere großräumige Ausgleichungen geeignet. U. a. hat 
das ZEN mit dem Östseering der Baltischen Geodä- 
tischen Kommission, dessen Ausgleichung durch 
Ölander im Jahr 1949 zum Abschluß gelangt ist, 
vier aufeinanderfolgende Stationen, also drei Polygon- 
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‘seiten, gemeinsam. Es ist also Hoffnung vorhanden, 

daß demnächst durch Zusammenschluß des Ostsee- 
rings und des ZEN ein weiterer Schritt zur Schaffung 
e eines gesamteuropäischen Dreiecksnetzes erfolgen 
f werde. 


Die bei der Ausgleichung des ZEN gewonnenen Er- 
fahrungen und insbesondere das von Wolf ent- 
wickelte neue Verfahren werden gewiß für künftige 
Arbeiten ähnlicher Art von sehr wesentlicher, neue 
Wege weisender Bedeutung sein. 

Der Band 9 enthält zwei Abhandlungen, die sich auf 
die Zentrierung von und nach trigonometrischen Sig- 
nalen gemessener Horizontalwinkel auf das unter dem 
Signal vermarkte Zentrum beziehen. 

In der ersten Abhandlung von Wolf wird der Vor- 
gang bei der Ablotungsmessung beschrieben und die 
Aufgabenstellung präzisiert. Der zur Bestimmung der 
Zentrierungselemente führende Rechnungsgang wird 
aufgezeigt, wobei die Lösung für zwei Fälle — ohne und 
mit Berücksichtigung überschüssiger Ablotestände — 
angegeben wird, und es wird eine Diskussion der Fehler- 
einflüsse durchgeführt. Es ergibt sich, daß größere, 
etwa 100—150 m überschreitende, Entfernungen der 
Ablotestände vom Signal tunlichst zu vermeiden sind, 
und daß es bei vorgegebenem Quantum der für die ge- 
samte Ablotung anzusetzenden Arbeit vorteilhafter ist, 
bei geringerer Satzanzahl der Winkelmessungen die 
Anzahl der Ablotestände zu erhöhen, als in umge- 
kehrter Weise zu verfahren. ( 

Im zweiten Aufsatz beschäftigt sich Mulert mit 
dem (auch bei der, Ausgleichung des ZEN häufig vor- 
gekommenen) Fall, daß auf nicht koordinierte Punkte 
bezügliche Beobachtungen auf koordinierte Punkte zu 
reduzieren sind. Es werden Formeln und Tafeln zur 
Berechnung der Zentrierungsbeträge aus den betreffen- 

den Gauß-Krüger- Koordinanten angegeben. 


Dresden. A: Buchholtz. 


Hugo Sirk (o. Prof. a.d. Univ. Wien), Mathema - 
tik für Naturwissenschaftler und Chemiker. 
6. verbesserte Aufl., XII u. 3028. m. 126 Abb. u. 
1 Ausschlagtafel. Dresden u. Leipzig 1950. Th. Stein- 

’ kopff. Preis 12,— DM. ; 

2 - Das Werk hat erhebliche Verbreitung gefunden; 
und es mag daher am Platze sein, es nochmals vom 
Standpunkt der Mathematik zu prüfen. Wir müssen 
uns dabei vor Augen halten, wozu der junge Chemiker 
die Mathematik braucht: Er muß soviel Mathematik 

‘ verstehen, daß ihm die Grundtatsachen der Thermo- 
dymamik, der theoretischen Chemie und schließlich 
auch der Vorstellungen von Atomhülle und Atom- 
kern im Laufe seiner Studienzeit überhaupt erschlos- 
sen werden können; wir sind der Überzeugung, laß 


ersten Studienjahr, mit Mathematik bekannt zu ma- 
chen. Statt dessen sehen wir allenthalben Bücher er- 
scheinen, die dem Chemiker das einfache Rüstzeug 
vorenthalten, oder in einer Form darbieten, welche 
mit Umstand und Ballast überladen ist, welche davor 
scheut die kleinen Hürden zu Ama „eher Khnllei 
sie aus in Sümpfe und Dorngestrüpp und begünstigen _ 
ee EEE Haltung. Man 
' Entgleisungen bei Asmus (de Gruyter) 
ren Werken, die sogar mehrere unverän- _ 
zeitigt haben. — Angesichts solcher 
Banaich die Erstauflage von Sirk 
en Versuchen dieser Art‘ gestellt. 
izise Einwände geschehen. Die 
uf solche Einwände keine 
mer noch wird der absolute 
a er 
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es richtig wäre, die Chemiker heute frühzeitig, im ni 
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Betrag vermieden und schon dadurch eine schleppende, 
stellenweise unexakte Darstellung herbeigeführt. Die 
Kunst, umständliche Formeln zu drucken, ist ge- 
rade dort besonders verbreitet, wo die Kunst fehlt, 
auch nur einfache Formeln zu lesen. Diese Erschei- 
nung teilt das Sirksche Buch mit anderen „einfachen‘‘ 
Einführungen in die Mathematik für Nicht-Mathe- 
matiker. Wir nennen nur zwei Leseproben: S. 155/56, 
wo erst die Differentiation von e”* verderbt darge- 
stellt ist (solcher Entstellungen der logischen Ordnung 
gibt es mehr), und dann die Formel i 


RE 1 © 
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mit typischem Mißbrauch anfänglicher Formel- 


schreiberei entwickelt ist: z* und «&° werden mit- 
geschleppt um.xt zu erhalten. S. 158 wird die Reihen- 
entwicklung für log (a + £) mit so umfänglichen For- 
meln gewonnen, daß kein junger Chemiker die Ge- 
duld ‚haben dürfte durchzuhalten. — Immer noch 
scheuen die Chemischen Rennpferde vor der Hürde 
des Grenzprozesses, dessen „Erwähnung in der täg- 
lichen Praxis zu umständlich‘‘ wäre (S. 24): darum 
verwende man den Begriff des Unendlich kleinen. 
In der Volkskunde kennt man die Erscheinung des 
„gesunkenen Kulturguts‘‘ (Hans Naumann): Was 
kulturell schöpferische Schichten geschaffen, das brei- 
tet sich langsam, unter sinkender Intensität in 
breitere, nachahmende Schichten aus. Es scheint 
uns, daß die Chemiker-Mathematik ein besonders ein- 
drucksvolles Beispiel ist. Hoffen wir, daß der einheit- 
liche, und nach vielen Seiten schmiegsame Grenz- 
begriff innerhalb der nächsten 50 Jahre dorthin ein- 
sinkt. — Es kann nicht ünsere Aufgabe sein, die vielen 
kleinen und größeren Mängel der Darstellung aufzu- 
zählen. Nach 6 Auflagen können wir das Buch zwar 
noch als besser als manche andre — doch darum 
keineswegs als gut ansehen. Es ist so wie es ist: 
irreparabel. 


Gießen. Egon Ullrich. 


Dr. 6. Hamel (em. o. Prof. a.d. Techn. Univ. 
Berlin, Theoretische Mechanik. Eine 
einheitliche Einführung in die gesamte Mechanik. 
(Die Grundlehren der mathematischen Wissenschaften 
in Einzeldarstellungen. Band L VIL) XVL-+ 796 S. 
mit 161 Abb. Berlin-Göttingen-Heidelberg 1949. 
Springer-Verlag. Preis brosch. 63,— DM, geb. 
66,— DM. 

Der Verf., langjähriger Professor der Mechanik an 
der Technischen Hochschule Charlottenburg, der 
jetzigen Technischen Universität, und einer der 
führenden Vertreter der theoretischen Mechanik in 
Deutschland gibt nunmehr den Hauptinhalt seiner 
Vorlesungen in Buchform heraus. 

Die Bezeichnung theoretische Mechanik anstelle 


der üblichen Bezeichnung ‚analytische Mechanik‘ 


wurde gewählt, um den deduktiven Charakter der 
entwickelten Lösungsmethoden zu betonen und es zu 
vermeiden, daß etwa die Analysis, z.B. die Inte- 
rationstheorie der Differentialgleichungen, in den 
Sekrond gerückt würde. Andererseits sollte der 
Buchtitel klar den Gegensatz zu den nicht-deduk- 
tiven, bzw. syntestischen und mehr anschaulichen 


* Methoden heraustellen, welche auf dem Gebiete der 


technischen Mechanik meist Anwendung finden. Der 
Verf. hat sich mit Erfolg bemüht, die gesamte theo- 
retische Mechanik in möglichst einheitlicher Dar- 
stellung erscheinen zu lassen, und zwar nicht getrennt 
in Punktmechanik und Mechanik der Kontinua, son- 
dern eingeteilt nach Gesichtspunkten, die sich einer- 


seits aus den Begriffsdefinitionen, anderseits aus der 
"Art des Bewegungszustandes der betrachteten Sy- 


steme und der zugehörigen mathematischen Behand- 
lung ergeben. Die einzelnen Kapitel tragen folgende 
Überschriften: Der Begriff der Kraft und das New- 


to nsche Grundgesetz; Statik gebundener Systeme 
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’ 
von endlichem Freiheitsgrad; Statik der Systeme von 
unendlich vielen Freiheitsgraden; die ersten allge- 
meinen Prinzipien der Kinetik; holonome Systeme 
mit endlichem Freiheitsgrad und die Lagrange- 
schen Gleichungen; mathematische Durcharbeitung; 
die Minimal-Prinzipien; der starıe Körper im Raum; 
holonome Systeme von endlichem Freiheitsgrad ; Auf- 
gaben und Probleme der theoretischen Mechanik. Es 
war also notwendig, nicht die eigentliche Stoffein- 
teilung der Mechanik, die der Verf. in seiner Fin- 
leitung besonders herausstellt (Kinematik, Massen- 
kinematik, Dynamik, Statik, Kinetik, Kinetostatik) 
für die Gliederung des Buches als maßgeblich zu ver- 
wenden; sondern es wurde durch dıe einheitliche 
Behandlung des Gesamtgebietes eine abweichende 
Aufteilung erforderlich, wobei gemäß dem Vorbild 
vonLagrangedas Prinzip der virtuellen Arbeiten, 
dasD’Alembertsche PrınzipunddasLagran- 
3esche Prinzip der Befreiung konsequent durch- 
geführt wurden. Hierin unterscheidet sıch das Buch 
wesentlich von dem Enzyklopädieartikel von Prange 
und der analytischen Dynamik von Whittaker, 
wenn sich auch der Verf. in seiner Darstellung zu- 
weilen an diese beiden Standardwerke anlehnt. Um 
die gestellte Aufgabe erfüllen zu können, ohne das 
Buch umfangsmäßig in ein Handbuch ausarten zu 
lassen, mußte der Verf. sich auf das beschränken, 
was für ein Lehrbuch über die Gesamtmechanik als 
wesentlich anzusehen ist. Dennoch zeigt sich eine 
große Vielseitigkeit des Inhaltes; u.a. finden sich 
Minimal-Prinzipien von Castigliano, sowie die 
Theorie dünner Häute, Schalen und Platten in den 
Abschnitten über Elastizitätstheorie, die Behandlung 
spezieller nicht-holonomer Systeme von endlichem 
Freiheitsgrad, Bemerkungen zur Integrationstheorie 
der kanonischen Gleichungen, die Anwendung der 
Lagrangeschen Methoden in der Kontinuums- 
mechanik. Auch eine kurze elegante Darstellung des 
Gravitationsfeldes in der allgemeinen Relativitäts- 
theorie wird gegeben. Vor allem enthält das Buch in 
seinem letzten Teil eine Sammlung interessanter Auf- 
gaben und Probleme‘ des Gesamtgebietes, welche 
absichtlich aus den vorangegangenen theoretischen 
Ableitungen herausgelöst und mit in das Schluß- 
kapitel aufgenommen wurden, so z. B. die Unter- 
suchungen über Reibung von Painleve& u.a, 
Behandlung des Mehrkörperproblems in einer erwei- 
terten Formulierung der elementaren Lösung von 
Lagrange, Aufstellung der allgemeinen Kreisel- 
beziehungen mit Hilfe der Lagrangeschen Glei- 
chungen unter Verwendung komplexer Koordinaten, 
mathematische Probleme zur Trennung der Variablen. 
Das interessante Buch, das offenbar einen Teil der 
Lebensarbeit des bekannten Verf. darstellt, stellt sich 
somit nicht nur würdig an die Seite der beiden erwähn- 
ten Standardwerke, sondern liefert in mancherlei Hin- 
sicht wertvolle Ergänzungen, die von der Fachwelt 
sicherlich begrüßt werden. Die Darstellung ist über- 
sichtlich und trotz der mathematischen Strenge klar 
und verständlich. Im Gegensatz zu den Werken von 
Whittaker und Prange, welche fast ohne 
Abbildungen erschienen waren, enthält das Ha- 
melsche Werk auch eine größere Anzahl erläutern- 
der. Skizzen, welche dem Leser das Verständnis er- 
leichtern. 


Dresden. 


H. Neuber. 


Lothar Hefiter, Grundlagen u. analy- 
tischer Aufbauder Projektiven, Eu- 
klidischen, Niehteuklidischen Geo- 
metrie. 2. Aufl., 1928. Leipzig 1950 B. G. Teuh- 
ner Verlagsgesellschaft. Preis geb. 12,60 DM. 

Das Buch ist im wesentlichen des dem dreibändigen 
Heffterschen ‚‚Lehrbuch der analytischen Geometrie‘* 
eigentümlichen systematischen Aufbaues; es bringt 
darüber hinaus aber mancherlei Verbesserungen u. 
Zusätze, die sich weder im ‚Lehrbuch‘ noch in 


ee, Buehbespreehungen 


‚allgemeinen ohne grundlegende Ableitungen. Es legt 


Stichworte geben auch die Hauptteile 


v 


anderen Veröffentlichungen finden. Es vereinigt in 
sich in hervorragendem Maße alle Vorzüge zur Eımög- 
lichung einer gründlichen Orientierung über die Pro- 
bleme der modernen Geometrie auf kürzestem Wege. 
In knapper und klarer Form werden, stilistisch stets 
‚mit größter Eindringlichkeit, die Grundlagen für die 
Klassifizierung der Geometrien herausgearbeitet, so 
daß der Leser ein plastisches Bild von den Unter- 
schieden und inneren Zusammenhängen gewinnt. . 

An die Darstellung der Grundlagen der Geometrie, 
in der alle Textstellen, die sich auf den Anschauungs- 
raum beziehen, für den axiomatischen Weg also nicht 
benötigt werden, auch äußerlich durch anderen Druck 
kenntlich gemacht werden, schließt sich die Behand- 
lung der Parallel- und Orthogonalgeometrie an. Im 
letzten Abschnitt wird, ebenfalls von der projektiven 
Geometrie aus, die Cayley-Kleinsche Geometrie ent- 
wickelt. 

Wollte die Besprechung Einzelheiten des reichen 
Inhaltes herausheben, so müßte sie bei der Gleich- 
wertigkeit fast auf alles eingehen. Dasist unmöglich. 
Darum sei mit einer warmen Empfehlung des Werkes 
geschlossen. 3 


Dresden. Draeger., 
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Otto Lorenz, DerIngenieurberuf, Vor- 
aussetzung— Ausbildung— Laufbah- 
nen. VI -+ 1538. m. 4Abb. Düsseldorf 1950. 
Deutscher Ingenieur-Verlag G.m.b. H. Preis brosch. 
4,80 DM für VDI Mitglieder 4,30 DM. 4 ! | 

Die vorliegende Schrift gibt einen guten Überblick 2 
über die verschiedenen Arten des Ingenieur-Berufes- 
für solche Kreise, die der Technik an sich fern stehen. ö 
Wer über die eigentliche Problematik des Studiums \ 
unterrichtet sein will, wird nicht befriedigt sein; ü 
ebenso sind die Angaben über die derzeitigen Ein- 
richtungen in der Ostzone unvollständig. 


Dresden. Heidebroek. 


Dr, phil. V. Happach, Ausgleichsrechnung. 
Ein Lehrbuch der Fehlerausgleichung nach der Me- 
thode der kleinsten Quadrate in Wissenschaft und _ 
Technik. (Teubners Mathematische Leitfäden Bd.18) 
2. Aufl, 104S. m. 26 Abb. Leipzig 1950. B.G. 
Teubner Verlagsgesellschaft. Preis kart. 6,20 DM. 

Das vorliegende Buch gibt die für die einzelnen 
Rechnungsgänge nötigen Formeln und Verfahren im 


u 
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das Hauptgewicht auf die Erläuterung ihrer An- 
wendung an zahlreichen aus verschiedenen Gebieten, 
insbesondere aus der geodätischen Vermessungstech- 
nik entnommenen Beispielen, und gibt so eine gute 
Einführung in die Technik der Ausgleichsrechnung.. 


‚Dresden. \ Wi llers. 


Bodo von Borries, Die Übermikroskopie. 
Einführung, Untersuchung ihrer Grenzen und Abriß 
ihrer Ergebnisse. 416 Seiten mit 225 Abbildungen. 
Berlin 1949. Verlag Dr. Werner Saenger. Preisgeb. 
485— DM. \ Br 3 IR 


Der Verfasser des vorliegenden Werkes, bekannt 


Ergebnisse der Übermikroskopie briı 
an. Die Einführung behandelt im wesent 
rechnerische Beweise die Grun 
optik, sie erläutert das Prinzip 
und die technische Ausführung un: 
gleich die verschiedenen 
präparation. Besonderer 
ausführliche Dars 
wicklung gelegt. 2 
der > 
jenigen Eigens« 
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E E r. Res er a 1950; Eingegangene Bücher 8 Nachrichten 399 
Bildfehler und Beugung, sowie technische Unvoll- B. Eck, Technische Strömungslehre. 
kommenheiten der Linsen. Ferner sind diejenigen Dritte, 2 Bahn und erweiterte Auflage. 
Leistungsgrenzen behandelt, die aus der Wechsel- X + 398 Seiten mit 372 Abb. Berlin- Göttingen- 


wirkung zwischen Elektronenstrahlen und Objekt 

entstehen: die beschränkte Dicke der als durchstrahlt 
erscheinenden Objekte; die untere Grenze des klein- 
sten mit hinreichendem Kontrast abgebildeten Details, 
die durch Ionisierung oder Erwärmung des Objekts 
bestimmt ist. Der letzte Hauptteil schließlich enthält 
Ergebnisse aus Chemie und Technik, sowie aus Bio- 
logıe und Medizin mit zahlreichen Abblidungen. 

Als leichtverständliche Darstellung wird das Buch 
allen denen willkommen sein, die mit dem Übermi- 
kroskop arbeiten müssen und "sich in die Theorie des 
Gerätes nicht im einzelnen einarbeiten können. Der 
zweite Hauptteil über die Grenzen der Übermikro- 
skopie mit seiner ausführlichen Wiedergabe der Ar- 
beiten des Verfassers wird auch denen viel zu bieten 
haben, die sich für die konstruktive Ausgestaltung 
des Gerätes interessieren. Tabellen wie auf Seite 240 

- bis 241 oder die Abbildung auf Seite 369 könnten nach 
Meinung des Referenten fortbleiben, ohne daß der 
Wert des Buches beeinträchtigt würde. Im Ganzen 
darf aber die Arbeit als gut gelungen bezeichnet 
werden. 


Dresden. Recknagel. 


Heidelberg 1949. Springer-Verlag. Ladenpreis: 24,— 
DM, Ganzleinen 27,— DM. 

Das Buch von Eck hat das Ziel, in einer dem 
Ingenieur leicht zugänglichen Form diejenigen Lehren 
und Versuchsergebnisse der: modernen Strömungs- 
forschung darzustellen, die für die Praxis wichtig 
erscheinen. Die jetzige dritte Auflage des beifällig 
aufgenommenen Werkes enthält wesentliche Erwei- 
terungen, so Ausführungen über die Aerodynamik der 
Verbrennung, über die Bewegung von Schwebe- 
teilchen, die z. B. für die pneumatische Förderung 
und die Entstaubungstechnik wichtigist, und manches 
andere mehr. Die steten Bemühungen des Verfassers, 
in seinem Buche mit den Fortschritten der Strömungs- 
wissenschaft Schritt zu halten, sind sehr anzuerkennen. 
Die Sichtung und Verarbeitung der in der umfäng- 
lichen Originalliteratur verstreuten neueren Ergeb- 
nisse und ihre eingängige Darstellung ist eine be- 
deutende Leistung, für die der in der Praxis stehende 
Ingenieur dem Verfasser Dank wissen wird. 


Göttingen. W. Tollmien. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen, 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


‚6. Zwikker (Technical Director Light Division 
‚Philips, Eindhoven), Advanced Plane Geometry 
VII + 299 S. mit 273 Abb. Amsterdam 1950. North- 
Holland Publishing Company. 


Dr.-Ing. Franz Ollendorff (Professor am Hebrew 
Technical College, Haifa). Die Welt der Vektoren. 
Einführung in Theorie und Anwendung der Vektoren, 
Tensoren und Operatoren. VIII -+ 470 S. mit 68 Abb. 
Wien 1950. Springer-Verlag. Preis: 37,50 DM., geb. 
40,— DM. 


Dr. phil. Ad. Duschek (o. Prof. a. d. Techn. Hoch- 
schule Wien) und Dr. techn. A. Hochrainer. Grund- 
züge der Tensorrechnung in analytischer 
Darstellung. Teil II: Tensoranalysis. VII + 338 S. 
mit 64 Abb. Wien 1950. Springer-Verlag. Preis: 
24,80 DM. 


Technische Bücher 1945 bis März 1950. Berlin- 
Göttingen-Heidelberg: Springer-Verlag. 


F, Strecker, Praktische Stabilitäts- 
prüfung mittelsOrtskurven und nu- 
merischer Verfahren. XI+189S. mit 101 
Abb. Berlin- Göttingen-Heidelberg 1950. Springer- 
Verlag, Preis brosch. 15,— DM, geb. 18,— DM. 


Herbert Beckert. 


8 
U 


Hristeng- und Eindeutigkeits- 


fie Di u u Da , 25 


“ RR von Prof, Sommerfeld an den inter- 
nationalen Mathematiker-Kongreß in Cambridge Mass. 
Die Ans 
tion für angewandte Mathematik von Herrn J. von 
Neumann i in engl. Übersetzung verlesen. 


Re = q Meine Damen undHerren! 


Ich bin sehr. traurig, daß ich heute nicht unter 
se und persönlich zu Ihnen sprechen kann, 
ir das ehrenvollerweise ‚von der TasunE dieses 


'beweise für das ee zur Lösung des » 


prache wurde in der ersten Sitzung der Sek- 


für meine Stellung zur an 


Anfangwertproblems, des gemischten Anfangs-, Rand- 
wert- und des charakteristischen Problems einer 
hyperbolischen Ditferentialgleichung zweiter Ordnung 
mit zwei unabhängigen Variablen. (Bericht über die 
Verhandlungen der Sächs. Akademie der Wissen- 
schaften zu Leipzig. Mathematisch-naturwissen- 
schaftl. Klasse, Bd. 97, Heft 4. 42 S. mit 9 Abb. 
Berlin: Akademie- Verlag 1950. Preis: 9,— DM. 


. Mathematisch-Physikalische- Semesterberichte zur 
Pflege des Zusammenhangs von Schule und Uni- 
versität. Herausgeg. von H. Behnke und W. Lietz- 

mann in Verbindung mit dem mathematischen For- 
schungsinstitut in Oberwolfach. Bd.1, 318 S. Göt- 
tingen 1949/1950. Vandenhoek & Ruprecht. Preis: 
19,60 DM. (Subskriptionspreis 17,— DM.) 


Mitteilungen aus dem Max-Planck-Institut 
für Strömungsforschung, 

Unter Mitwirkung von L. Prandtl, W. Tollmien, 
G.Vogelpohl herausgegeben von A. Betz. Göttingen 
1950. Selbstverlag Max-Planck-Institut für Strö- 
mungsforschung. 

Heft 1: J. Rotta, Über die Theorie der turbu- 
lenten Grenzschicht. 54 S. mit 10 Abb. Preis: 
4,— DM. 

Heft2: H, Himmelskamp, Profiluntersuchun - 

en an einem umlaufenden Propeller. 38 8. 
mit 59 Abb. Preis: 8,— DM 


NACHRICHTEN 


Kongresses vorgeschlagen war. Ich bin dreimal in 
Ihrem Lande gewesen, 1922/23 als Carl-Schurz-Pro- 
fessor in Madison Wis, 1929 als Gastprofessor in 
Pasadena, 1931 bei der Sommersession in Ann. Arbor. 
und habe mich jedesmal unter meinen vielen ameri- 
kanischen Freunden und Schülern herzlich wohl ge- 
fühlt. 

- Ich hätte auch gern persönlich Zeugnis abgelegt 
nandten Mathematik. 


Die Zeiten sind längst vorüber, in Ihrem Lande wie 
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in meinem, wo ein bedeutender englischer Mathe- 
matiker sagte: Gott bewahre die Mathematik vor ihren 
Anwendungen! Wir wissen alle, daß die abstrakte 
Mathematik ihre eigenen Wege gehen muß, ohne 
Rücksicht auf mögliche Anwendungen, die oft lange 
auf sich warten lassen. Wir wissen aber auch, mit 
welcher Liebe sich die größten Mathematiker, Gauss, 
Riemann, Cauchy, Poincare den Anwendungen ge- 
widmet haben. 

Als Planck kurz vor meiner letzten Reise nach 
Amerika mir die ZPlanck-Medaille überreichte, 
apostrophierte er mich, sich gewissermaßen entschul- 
digend, als angewandten Mathematiker. Ich ant- 
wortete ihm, daß ich diese Benennung gern acceptiere 
und daß ich nicht so denke, wie mein Kollege Alfred 
Pringsheim, der Weierstrass-Schüler, der einmal bei 
einem‘ studentischen Feste sagte: Sie wissen, daß 
ich immer einen gewandten Mathematiker mehr 
geschätzt habe als einen angewandten Mathe- 
matiker. Hoffentlich konnten Sie die Antithese ge- 
wandt, angewandt, ungewandt, die ganz im Sinne 
des witzigen br ingsheim lag, auch in der Übersetzung 
verstehen. Planck gegenüber betonte ich, mein 
Streben sei immer gewesen, zugleich ein gewandter 
und angewandter Mathematiker zu sein. 

Mein Freund Constantin Carathedory, mit dem 
ich gehofft hatte gemeinsam zu diesem Kongreß zu 
reisen, war in den abstrakten Regionen der Mengen- 
lehre zu Hause und gleichzeitig erfolgreich in der 
Thermodynamik, der geometrischen ‘Optik und der 
allgemeinen Mechanik. Wir alle bedauern tief seinen 
plötzlichen Tod. Er wäre eine Zierde beider Sektionen 
dieses Kongresses, der reinen und der angewandten, 
gewesen. 

Für den umfassenden Geist Einsteins ist Mathe- 
matik, Physik, Kosmologie und Philosophie eine 
Einheit. Auch sonst beschränkt sich in Ihrem Lande 
die angewandte Mathematik nicht wie früher auf das 
Gebiet des numerischen Rechnens, der Astronomie 
und Geodäsie. Sie richtet sich auf die Ordnung aller 
physikalischen Erscheinungen und technischen Er- 
fahrungen, ja sie wagt es mathemathische Methoden 
in die Psychologie und die Technik des Denkens zu 
übertragen. 

Mein Wunsch ist, daß dieser Kongreß für die 
Mathematik im: weitesten und breitesten Sinne 
fruchtbar sein möge. 


Internationaler Mathematikerkongreß 


In der Zeit vom 30 August bis 6. September 1950 
fand an der Harvard-Universität in Cambridge, Mas- 
sachusetts USA, ein Internationaler Mathematiker- 
kongreß statt, der von über 2000 Personen besucht war. 
Es wurden über 500 Vorträge gehalten, darunter einige 
zusammenfassende Vorträge von zwanzig bis sechzig 
Minuten Dauer, alle anderen Vorträge waren Kurz- 
vorträge von zehn Minuten Dauer. Die Vorträge wur- 
den in parallel verlaufenden Sektionen abgehalten; es 
waren 7 Sektionen aufgestellt: I. Algebra und Zahlen- 
theorie, II. Analysis, III. Geometrie und Topologie, 
IV. Wahrscheinlichkeitsrechnung ünd Statistik, V. 
Mathematische Physik and angewandte Mathematik, 
VI. Logik und Philosophie, VII. Geschichte der Mathe- 
matik und Erziehung. Hier seien von den Vorträgen 
nur die zusammenfassenden Vorträge in der Sektion V 
genannt. Sie wurden gehalten von C.D. Darwin 
(Refractive Index ionisierter Gase), F., Rellich 
(Störungstheorie der Spektralzerlegung), J. v. Neu - 


. mann(Shock Interaction), R.Courant ren: - 


aufgaben ara modernen Hydrodynamik), 8. Gold- 


Verantwortlich für den Inhalt: Prof. Dr. Fr. A. Willers, Dresden; für 
GmbH., Berlin NW 7, Schiffbauerdamm 19, Fernsprecher: Sammelnummmeı 
VerlagsnummerdiesesHeftes: 1009/30/11/12. Die Zeitschrift für angew 
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erkannt werden, anderenfalls wird das folgende Quartal noch geliefert 
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zu 64 Zahlen) erfolgt in einem mechanisch arbeitenden 


‚gegeben, die Resultate ebenfalls Überrgun 


a ine (Aus der Gasdynamik), W. "Hei p en b ne 
(Stabilitätsfragen der Flüssigkeitsdynamik), W.Pra-. 
ger (Randwertprobleme der Plastizität), J. I. Sto- 
ker (Wasserwellen), H. Lewy (Analytische Rand- 


bedingungen). Neben den Vorträgen fanden Be- 
sichtigungen der Institute statt, so konnte man z.B. 
im MIT (Massachusetts Institute of Technology) einen 
Differential Analyser, ein Zyklotron usw. und im 
Computation Laboratory einen I.B.M. Automatic 
Sequence Controlled Calculator sehen. 


Hannover. L. Collatz. 


Sektion für Kristallkunde der Deutschen 
Mineralogischen Gesellschaft 


Auf der diesjährigen Hauptversammlung der Deut- 
schen Mineralogischen Gesellschaft, die Ende August 
in Göttingen stattfand, wurde auf ‚Anregung einer 
Reihe interessierter Wissenschaftler eine „Sektion für 
Kristallkunde‘“ eingerichtet und diese Sektion sinn- 
gemäß in den Satzungen der Deutschen Mineralogi-, 
schen Gesellschaft verankert. Zum Leiter der Sektion, 
der gleichzeitig einer der beiden stellvertretenden Vor- 
sitzenden der Deutschen Mineralogischen Gesellschaft _ 
ist, wurde Professor Dr. H.O’Daniel, a N 
Frankfurt, gewählt. 

Die Gesellschaft wünscht, in der „Sektion für 
Kristallkunde‘“ auch diejenigen Fachrichtungen zu 
pflegen, die — über den Rahmen speziell mineralogi- 
scher Fragestellungen hinausgehend — sich mit dem 
Aufbau und den Eigenschaften des festen ‘/kristalli- 
sierten Körpers befassen. Die Gesellschaft hofft, auf 
diesem Wege einen engen Kontakt mit den Nachbar- 
disziplinen der Mathematik, der Physik und der Che- 
mie herzustellen. Eine jährliche 'Vortragstagung soll 
vorzüglich hierzu dienen. 


GaMM-Tagung 2 
Die wissenschaftliche Jahrestagung der Gesellschaft 
für angswandte Mathematik und Mechanik findet in 
der Zeit vom 28. bis 31. März 1951 in Freiburg 
i. Br. statt. Anmeldungen zur Teilnahme sowie An- 
meldungen von Vorträgen (mit kurzer Inhaltsangabe 
bis zu 15 Zeilen) sind an den örtlichen Tagungsleiter, 
Prof. Dr. Görtler (17b) Freiburg i. Br., Stadtstr. 57 
zu richten, und zwar beides bis spätestens 28. Be- 
bruar 1951. Als normale Vortragsdauer sind 10 bis 
15 Minuten vorgesehen. Die Mr N im 
Hauptgebäude der Universität statt, . 


Z a Im Juli 1950 wurde am Institut für ange- 
wandte Mathematik der ETH. Zürich eine von der 
Zuse KG. in Neukirch (Kreis Hünfeld) gebaute pro- 
grammgesteuerte Relais-Rechenmaschine „Modell Z4“ 
in Betrieb genommen. Das elektro-mechanisch arbei-r 
tende Gerät hat 2200 Telephon-Relais und 21 Schritt- 
schalter. Die Speicherung der Zwischenresultate (bis 


Speicherwerk. Die Rechnung erfolgt im Dualsystem 
mit gleitendem Komma. Die Daten werden mittels 
Tasten oder Lochkarten in ee 


gelocht oder gedruckt. Die Üb 
die Dualschreibweise 7a ‚autom 
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